Coupled Fields
The definition of each coupled field can be found at the CCCma TWiki.
For those outside of the CCCma, the raw TWiki code is included below:
---+ <center> !CanAM5 Coupled with NEMO/LIM2 using the conservative options </center>
<center>
%TABLE{caption="<font size=+1><b>Representative Areas</b></font>" headerrows="1" }%
| *Index* | *Name* | *Area Description* |
| 1 | C | the entire grid cell |
| 2 | W=OI | the non-land (water) portion of the grid cell |
| 3 | O | the open water portion of the non-land portion |
| 4 | I | the ice covered portion of the non-land portion |
</center>
<button class="twistyExpandAll twikiButton">Expand all</button> <button class="twistyCollapseAll twikiButton">Collapse all</button>
%TWISTY{mode="div" showlink="About the *Field Sent* and *Field Received* columns in these tables." hidelink="<font color=#cf00ff>Hide About fields sent/received</font>" remember="on"}%
The *Field Sent* and *Field Received* columns in the following tables each contain 2 pieces of information on separate lines.
The first is the name of the Fortran variable in the relevant component model (either atmosphere or ocean) containing the data that is sent or received.
The second is the name associated with this data in the coupler.
If the data is sent to or received from NEMO then the coupler name will be followed by the name of the constant used in NEMO to index the =ssnd= or =srcv= array associated with the field that is sent or received, e.g. for sea ice fraction this would be !OIceFrc=(jps_fice).
The *Source* and *Destination* columns contain similar information to the *Field Sent* and *Field Received* columns. These provide continuity between complimentary tables. Tables 1 and 2 are complimentary as are tables 3 and 4. A table with a *Field Sent* column will also contain a *Destination* column and a table with a *Field Received* will also contain a *Source* column.
In some cases a source field contributes to more than one destination field. In these cases there will be more that one Fortran variable name on the first line of the corresponding *Field Received* entry.
Conversely, more than one source field can contribute to a single destination field. In these cases there will be more than one Fortran variable name on the first line of the corresponding *Source* column entry.
%ENDTWISTY%
%TWISTY{mode="div" showlink="Table 1. Fields Sent by the NEMO Ocean" hidelink="<font color=#cf00ff>Hide Table 1.</font>" remember="on"}%
%TABLE{caption="<font size=+1><b>Table 1. Fields Sent by the NEMO Ocean</b></font>" headerrows="1" }%
| *Row* | *Config* | *Namelist Param* | *Field Sent* | *Field Description* | *Rep Area* | *Sign* | *Units* | *Grid* | *Time Avg* | *Pre-proc* | *Remap Options* | *Post-proc* | *Destination* | *Comments* |
| %CALC{"$ROW(-1)"}% | 2 | | fr_i<br />OIceFrc=(jps_fice) | sea ice fraction | OI | | 1 | T | Inst | | cv-fracarea<br />mask orca land | | SICNPAK<br />SICN_atm | |
| %CALC{"$ROW(-1)"}% | 2 | sn_snd_thick%cldes=<br />'weighted iwe and swe' | ht_i<br />OIceTck=(jps_hice) | ice water equivalent of sea ice<br />(avg over ocn and ice) | OI | | kg m-2 | T | Inst | | cv-fracarea<br />mask orca land | | SICPAK<br />SIC_atm | constant ice density<br />ρ<sub>i</sub> = 900 kg m-3 |
| %CALC{"$ROW(-1)"}% | 2 | sn_snd_thick%cldes=<br />'weighted iwe and swe' | ht_s<br />OSnwTck=(jps_hsnw) | snow water equivalent over sea ice<br />(avg over ocn and ice) | OI | | kg m-2 | T | Inst | ht_s/fr_i | cv-fracarea<br />mask orca land | | SNOPAKO<br />SNO_atm | constant snow density<br />ρ<sub>s</sub> = 330 kg m-3 |
| %CALC{"$ROW(-1)"}% | 2 | sn_snd_temp%cldes=<br />'oce and ice' | tsn+273.15<br />O_SSTSST=(jps_toce) | sea surface temperature | O | | K | T | Inst | | cv-fracarea<br />mask orca land | | GTPAT<br />GT_atm | |
| %CALC{"$ROW(-1)"}% | 2 | sn_snd_temp%cldes=<br />'oce and ice' | tn_ice<br />O_TepIce=(jps_tice) | sea ice surface temperature | I | | K | T | Inst | | no remap | | coupler | |
| %CALC{"$ROW(-1)"}% | 2 | sn_snd_co2%cldes=<br />'none' | oce_co2<br />O_CO2FLX=(jps_co2) | CO2 flux | OI | | mol C m-2 s-1 | T | Inst | | cv-fracarea<br />mask orca land | *0.044011*-1.0 | XSFXPAK(1,ICO2)<br />CO2flx_atm | #ifdef key_cpl_carbon_cycle |
%TWISTY{mode="div" showlink="<font color=#cf00ff>Table 1 Notes</font>" hidelink="<font color=#cf00ff>Hide Notes</font>" remember="on"}%
The values of !OIceTck and !OSnwTck depend on the namelist parameter sn_snd_thick%cldes as follows
* If sn_snd_thick%cldes='ice and snow' then
* !OIceTck = ht_i
* !OSnwTck = ht_s
* If sn_snd_thick%cldes='weighted ice and snow' then
* !OIceTck = ht_i*fr_i
* !OSnwTck = ht_s*fr_i
* If sn_snd_thick%cldes='weighted iwe and swe' then
* !OIceTck = rhoic*ht_i*fr_i ...constant rhoic=900 kg m-3
* !OSnwTck = rhosn*ht_s*fr_i ...constant rhosn=330 kg m-3
The 'weighted iwe and swe' option was added so that ice and snow water equivalent values may be evaluated on the ORCA grid and passed directly to the AGCM. No pre or post processing is required (except interpolation onto the AGCM grid).
The ice and snow densities required for the conversion are the constant values that are used in NEMO and defined in the =phycst= module.
The value for ocean surface temperature sent from NEMO to the coupler depends on the namelist parameter sn_snd_temp%cldes as follows
* If sn_snd_temp%cldes='oce only' then
* O_SSTSST = tsn+273.15
* If sn_snd_temp%cldes='weighted oce and ice' then
* O_SSTSST = (tsn+273.15)*(1-fr_i)
* O_TepIce = tn_ice*fr_i
* If sn_snd_temp%cldes='oce and ice' then
* O_SSTSST = tsn+273.15
* O_TepIce = tn_ice
* If sn_snd_temp%cldes='mixed oce-ice' then
* O_TepMix = (tsn+273.15)*(1-fr_i) + tn_ice*fr_i
The 'oce and ice' option was added so that the value of ice surface temperature (tn_ice) that is used in NEMO can be passed to the coupler along with sea surface temperature.
The ice surface temperature will be required when a bulk model is added to the coupler and was useful for debugging during coupler development.
Note that when the namelist parameter sn_snd_temp%cldes='oce only' or 'oce and ice' the sea surface temperature received by the coupler, O_SSTSST, is defined over the entire ocean including under sea ice.
In this sense the representative area is not strictly ice free or partially ice covered ocean. This choice is made because the AGCM requires SST values everywhere.
In the case sn_snd_temp%cldes='weighted oce and ice', O_SSTSST is the ocean cell mean value and will be zero under grid cells that are entirely covered with ice.
%ENDTWISTY%
%ENDTWISTY%
%TWISTY{mode="div" showlink="Table 2. Fields Received by the !CanAM5 Atmosphere" hidelink="<font color=#cf00ff>Hide Table 2.</font>" remember="on"}%
%TABLE{caption="<font size=+1><b>Table 2. Fields Received by the !CanAM5 Atmosphere</b></font>" headerrows="1" }%
| *Row* | *Config* | *Field Received* | *Field Description* | *Rep Area* | *Sign* | *Units* | *Grid* | *Time Avg* | *Post-proc* | *Remap Options* | *Pre-proc* | *Source* | *Comments* |
| %CALC{"$ROW(-1)"}% | 2 | SICNPAK<br />SICN_atm | sea ice fraction | OI | | 1 | G | Inst | | cv-fracarea<br />mask orca land | | fr_i<br />OIceFrc=(jps_fice) | |
| %CALC{"$ROW(-1)"}% | 2 | SICPAK<br />SIC_atm | ice water equivalent of sea ice<br />(avg over ocn and ice) | OI | | kg m-2 | G | Inst | | cv-fracarea<br />mask orca land | | ht_i<br />OIceTck=(jps_hice) | |
| %CALC{"$ROW(-1)"}% | 2 | SNOPAT(1,IOSIC)<br />SNO_atm | snow water equivalent over sea ice<br />(NOT cell avg) | I | | kg m-2 | G | Inst | | cv-fracarea<br />mask orca land | ht_s/fr_i | ht_s<br />OSnwTck=(jps_hsnw) | |
| %CALC{"$ROW(-1)"}% | 2 | GTPAT<br />GT_atm | sea surface temperature | O | | K | G | Inst | | cv-fracarea<br />mask orca land | | tsn+273.15<br />O_SSTSST=(jps_toce) | |
| %CALC{"$ROW(-1)"}% | 2 | XSFXPAK(1,ICO2)<br />CO2flx_atm | CO2 flux | OI | | kg CO2 m-2 s-1 | G | Inst | | cv-fracarea<br />mask orca land | *0.044011*-1.0 | oce_co2<br />O_CO2FLX=(jps_co2) | #ifdef biogeochem |
%ENDTWISTY%
%TWISTY{mode="div" showlink="Table 3. Fields Sent by the !CanAM5 Atmosphere" hidelink="<font color=#cf00ff>Hide Table 3.</font>" remember="on"}%
%TABLE{caption="<font size=+1><b>Table 3. Fields Sent by the !CanAM5 Atmosphere</b></font>" headerrows="1" }%
| *Row* | *Config* | *Field Sent* | *Field Description* | *Rep Area* | *Sign* | *Units* | *Grid* | *Time Avg* | *Pre-proc* | *Remap Options* | *Post-proc* | *Destination* | *Comments* |
| %CALC{"$ROW(-1)"}% | 2 | UFSOPAL<br />UFSO_atm | Atm-ocn wind stress (x) | O | pos up | N m-2 | G | avg | -UFSO_atm | cv-fracarea<br />mask orca land | | utau<br />O_OTaux1=(jpr_otx1) | |
| %CALC{"$ROW(-1)"}% | 2 | VFSOPAL<br />VFSO_atm | Atm-ocn wind stress (y) | O | pos up | N m-2 | G | avg | -VFSO_atm | cv-fracarea<br />mask orca land | | vtau<br />O_OTauy1=(jpr_oty1) | |
| %CALC{"$ROW(-1)"}% | 2 | UFSIPAL<br />UFSI_atm | Atm-ice wind stress (x) | I | pos up | N m-2 | G | avg | -UFSI_atm | cv-fracarea<br />mask orca land | | utau_ice<br />O_ITaux1=(jpr_itx1) | |
| %CALC{"$ROW(-1)"}% | 2 | VFSIPAL<br />VFSI_atm | Atm-ice wind stress (y) | I | pos up | N m-2 | G | avg | -VFSI_atm | cv-fracarea<br />mask orca land | | vtau_ice<br />O_ITauy1=(jpr_ity1) | |
| %CALC{"$ROW(-1)"}% | 2 | FSGOPAL<br />FSGO_atm | Solar heat flux over ocean | O | pos down | W m-2 | G | avg | | cv-fracarea<br />mask orca land | | qsr_tot<br />O_QsrOce=(jpr_qsroce) | |
| %CALC{"$ROW(-1)"}% | 2 | FSGIPAL<br />FSGI_atm | Solar heat flux over ice | I | pos down | W m-2 | G | avg | | cv-fracarea<br />mask orca land | | qsr_tot, qsr_ice<br />O_QsrOce=(jpr_qsroce) | |
| %CALC{"$ROW(-1)"}% | 2 | BEGOPAL<br />BEGO_atm | Total heat flux over ocean | O | pos down | W m-2 | G | avg | BEGO_atm - FSGO_atm | cv-fracarea<br />mask orca land | | qns_tot<br />O_QnsOce=(jpr_qnsoce) | |
| %CALC{"$ROW(-1)"}% | 2 | BEGIPAL<br />BEGI_atm | Total heat flux over sea ice | I | pos down | W m-2 | G | avg | BEGI_atm - FSGI_atm | cv-fracarea<br />mask orca land | | qns_tot, qns_ice<br />O_QnsIce=(jpr_qnsice) | |
| %CALC{"$ROW(-1)"}% | 2 | RAINSPAL<br />RAIN_atm | Total liquid precipitation | C | pos down | kg m-2 s-1 | G | avg | | cv-fracarea<br />mask orca land | | emp_tot<br />OTotRain=(jpr_rain) | |
| %CALC{"$ROW(-1)"}% | 2 | SNOWSPAL<br />SNOW_atm | Total solid precipitation | C | pos down | kg m-2 s-1 | G | avg | | cv-fracarea<br />mask orca land | | emp_tot, emp_ice, sprecip<br />OTotSnow=(jpr_snow) | |
| %CALC{"$ROW(-1)"}% | 2 | BWGOPAL<br />BWGO_atm | ocean freshwater budget (P-E) | O | pos down | kg m-2 s-1 | G | avg | -BWGO_atm | cv-fracarea<br />mask orca land | | emp_tot<br />OOEvaMPr=(jpr_oemp) | |
| %CALC{"$ROW(-1)"}% | 2 | BWGIPAL<br />BWGI_atm | sea ice fresh water budget | I | pos down | kg m-2 s-1 | G | avg | -BWGI_atm | cv-fracarea<br />mask orca land | | emp_tot<br />OSubMPre=(jpr_sbpr) | |
| %CALC{"$ROW(-1)"}% | 2 | SLIMPAL<br />SLIM_atm | non-Solar sensitivity to temperature over sea ice | I | ?? | W m-2 K-1 | G | avg | | cv-fracarea<br />mask orca land | | dqns_ice<br />O_dQnsdT=(jpr_dqnsdt) | |
| %CALC{"$ROW(-1)"}% | 2 | RIVOGRD<br />RIVO_atm | River discharge | OI | pos down | kg m-2 s-1 | G | avg | | cv-fracarea<br />mask orca land | | subtract from emp_tot<br />O_Runoff=(jpr_rnf) | #ifndef agcm_river_routing<br />RIVOGRD = 0<br />#endif |
| %CALC{"$ROW(-1)"}% | 2 | SWAPAL<br />SWMX_atm | Mixed 10 meter wind | C | N/A | m s-1 | G | avg | | cv-fracarea<br />mask orca land | | wndm<br />O_Wind10=(jpr_w10m) | |
| %CALC{"$ROW(-1)"}% | 2 | PMSLPAL<br />PMSL_atm | Mean sea level pressure | C | N/A | hPa | G | avg | PMSL_atm/1013.25 | cv-fracarea<br />mask orca land | | atm_mslp<br />O_MSLP=(jpr_mslp) | #ifdef biogeochem |
| %CALC{"$ROW(-1)"}% | 2 | XSRFPAL(1,ICO2)<br />CO2_atm | Atm CO2 concentration | C | N/A | mass mixing ratio | G | avg | 1e6 * 28.97/44.0 * CO2_atm | cv-fracarea<br />mask orca land | | atm_co2<br />O_AtmCO2=(jpr_co2) | #ifdef biogeochem |
%ENDTWISTY%
%TWISTY{mode="div" showlink="Table 4. Fields Received by the NEMO Ocean" hidelink="<font color=#cf00ff>Hide Table 4.</font>" remember="on"}%
%TABLE{caption="<font size=+1><b>Table 4. Fields Received by the NEMO Ocean</b></font>" headerrows="1" }%
| *Row* | *Config* | *Namelist Param* | *Field Received* | *Field Description* | *Rep Area* | *Sign* | *Units* | *Grid* | *Time Avg* | *Post-proc* | *Remap Options* | *Pre-proc* | *Source* | *Comments* |
| %CALC{"$ROW(-1)"}% | 2 | sn_rcv_tau%cldes=<br />'oce and ice' | utau<br />O_OTaux1=(jpr_otx1) | Atm-ocn wind stress (x) | O | pos down | N m-2 | T | avg | | cv-fracarea<br />mask orca land | -UFSO_atm | UFSOPAL<br />UFSO_atm | |
| %CALC{"$ROW(-1)"}% | 2 | sn_rcv_tau%cldes=<br />'oce and ice' | vtau<br />O_OTauy1=(jpr_oty1) | Atm-ocn wind stress (y) | O | pos down | N m-2 | T | avg | | cv-fracarea<br />mask orca land | -VFSO_atm | VFSOPAL<br />VFSO_atm | |
| %CALC{"$ROW(-1)"}% | 2 | sn_rcv_tau%cldes=<br />'oce and ice' | utau_ice<br />O_ITaux1=(jpr_itx1) | Atm-ice wind stress (x) | I | pos down | N m-2 | T | avg | | cv-fracarea<br />mask orca land | -UFSI_atm | UFSIPAL<br />UFSI_atm | |
| %CALC{"$ROW(-1)"}% | 2 | sn_rcv_tau%cldes=<br />'oce and ice' | vtau_ice<br />O_ITauy1=(jpr_ity1) | Atm-ice wind stress (y) | I | pos down | N m-2 | T | avg | | cv-fracarea<br />mask orca land | -VFSI_atm | VFSIPAL<br />VFSI_atm | |
| %CALC{"$ROW(-1)"}% | 2 | sn_rcv_qsr%cldes=<br />'conservative' | <font color="#0000a0">qsr_tot</font><br />O_QsrMix=(jpr_qsrmix) | solar heat flux mixed over ocean-ice | OI | pos down | W m-2 | T | avg | | cv-fracarea<br />mask orca land | SICN_atm*FSGI_atm + <br/> (1-SICN_atm)*FSGO_atm| FSGOPAL, FSGO_atm <br/> FSGIPAL, FSGI_atm| |
| %CALC{"$ROW(-1)"}% | 2 | sn_rcv_qsr%cldes=<br />'conservative' | qsr_ice, <font color="#0000a0">qsr_tot</font><br />O_QsrIce=(jpr_qsrice) | solar heat flux over sea ice | I | pos down | W m-2 | T | avg | | cv-fracarea<br />mask orca land | | FSGIPAL<br />FSGI_atm | |
| %CALC{"$ROW(-1)"}% | 2 | sn_rcv_qns%cldes=<br />'conservative' | <font color="#0000a0">qns_mix</font><br />O_QnsMix=(jpr_qnsmix) | non-solar heat flux mixed over ocean-ice | OI | pos down | W m-2 | T | avg | | cv-fracarea<br />mask orca land | (BEGO_atm - FSGO_atm)*<br/>(1-SICN_atm) + <br/> (BEGI_atm - FSGI_atm)*SICN_atm | BEGOPAL,FSGOPAL<br />BEGO_atm,FSGO_atm <br/> BEGIPAL,FSGIPAL<br />BEGI_atm, FSGI_atm| |
| %CALC{"$ROW(-1)"}% | 2 | sn_rcv_qns%cldes=<br />'conservative' | qns_ice, <font color="#0000a0">qns_tot</font><br />O_QnsIce=(jpr_qnsice) | non-solar heat flux over sea ice | I | pos down | W m-2 | T | avg | | cv-fracarea<br />mask orca land | BEGI_atm -<br />FSGI_atm | BEGIPAL,FSGIPAL<br />BEGI_atm, FSGI_atm | |
| %CALC{"$ROW(-1)"}% | 2 | sn_rcv_emp%cldes=<br />'conservative' | <font color="#0000a0">emp_tot</font><br />OTotEvap=(jpr_tevp) | Total evaporation (ocean evap + ice sublimation ) | OI | pos up | kg m-2 s-1 | T | avg | | cv-fracarea<br />mask orca land | RAIN_atm + SNOW_atm - <br/>( BWGI_atm*SICN_atm + BWGO_atm*(1-SICN_atm)) | RAINSPAL, SNOWSPAL, BWGIPAL, BWGOPAL<br />RAIN_atm, SNOW_atm, BWGI_atm<br />BWGO_atm | |
| %CALC{"$ROW(-1)"}% | 2 | sn_rcv_emp%cldes=<br />'conservative' | <font color="#0000a0">sprecip</font><br />OIceEvap=(jpr_ievp) | sublimation over sea ice | I | pos up | kg m-2 s-1 | T | avg | | cv-fracarea<br />mask orca land | RAIN_atm + SNOW_atm - BWGI_atm | RAINSPAL, SNOWSPAL, BWGIPAL<br />RAIN_atm, SNOW_atm, BWGI_atm | |
| %CALC{"$ROW(-1)"}% | 2 | sn_rcv_emp%cldes=<br />'conservative' | <font color="#0000a0">emp_tot</font><br />OTotSnow=(jpr_snow) | Snow | C | pos down | kg m-2 s-1 | T | avg | | cv-fracarea<br />mask orca land | | SNOWSPAL<br/>SNOW_atm | |
| %CALC{"$ROW(-1)"}% | 2 | sn_rcv_emp%cldes=<br />'conservative' | emp_tot<br />OTotRain=(jpr_rain) | Rain | C | pos down | kg m-2 s-1 | T | avg | | cv-fracarea<br />mask orca land | | RAINSPAL <br/>RAIN_atm | |
| %CALC{"$ROW(-1)"}% | 2 | sn_rcv_dqnsdt%cldes=<br />'coupled' | dqns_ice<br />O_dQnsdT=(jpr_dqnsdt) | non-Solar sensitivity to temperature over sea ice | I | ?? | W m-2 K-1 | T | avg | | cv-fracarea<br />mask orca land | | SLIMPAL<br />SLIM_atm | |
| %CALC{"$ROW(-1)"}% | 2 | sn_rcv_rnf%cldes=<br />'coupled' | <font color="#0000a0">emp_tot</font><br />O_Runoff=(jpr_rnf) | runoff | N/A | N/A | kg m-2 s-1 | T | avg | | cv-fracarea<br />mask orca land | | RIVOGRD<br />RIVO_atm | #ifndef agcm_river_routing<br />RIVOGRD = 0<br />#endif |
| %CALC{"$ROW(-1)"}% | 2 | sn_rcv_w10m%cldes=<br />'coupled' | wndm<br />O_Wind10=(jpr_w10m) | 10 meter wind | C | N/A | m s-1 | T | avg | | cv-fracarea<br />mask orca land | | SWAPAL<br />SWMX_atm | |
| %CALC{"$ROW(-1)"}% | 2 | sn_rcv_taumod%cldes=<br />'coupled' | taum<br />O_TauMod=(jpr_taum) | ocean wind stress modulus | O | N/A | N m-2 | T | avg | | cv-fracarea<br />mask orca land | sqrt(UFSO_atm<sup>2</sup> + VFSO_atm<sup>2</sup>) | UFSPAL,VFSPAL<br />UFSO_atm, VFSO_atm | |
| %CALC{"$ROW(-1)"}% | 2 | | atm_mslp<br />O_MSLP=(jpr_mslp) | Mean sea level pressure | C | N/A | atm | T | avg | | cv-fracarea<br />mask orca land | PMSL_atm/1013.25 | PMSLPAL<br />PMSL_atm | #ifdef key_cpl_carbon_cycle |
| %CALC{"$ROW(-1)"}% | 2 | sn_rcv_co2%cldes=<br />'none' | atm_co2<br />O_AtmCO2=(jpr_co2) | atm CO2 concentration | C | N/A | ppm | T | avg | | cv-fracarea<br />mask orca land | 1e6 * 28.97/44.0 * CO2_atm | XSRFPAL(1,ICO2)<br />CO2_atm | #ifdef key_cpl_carbon_cycle |
%TWISTY{mode="div" showlink="<font color=#cf00ff>Table 4 Notes</font>" hidelink="<font color=#cf00ff>Hide Notes</font>" remember="on"}%
In what follows the variable fr_i is sea ice fraction on the ORCA grid. Also the namelist element "clcat" which is used to indicate multiple ice categories is defined as 'no' in all cases for the NEMO/LIM2 configuration.
Other than that, if a namelist element is not defined here then it is assumed to contain the default value as defined in the subroutine sbc_cpl_init.
The following combination of namelist parameters will result in 4 stress fields associated with indices of the srcv array (jpr_otx1, jpr_oty1, jpr_itx1, jpr_ity1) being passed to the ocean.
These are the x and y components of atm-ocn and atm-ice wind stress on the T grid which will be rotated onto the ORCA grid and assigned to the NEMO variables utau, vtau, utau_ice and vtau_ice respectively.
This is what is used in the NEMO/LIM2 configuration.
* sn_rcv_tau%cldes='oce and ice'
* sn_rcv_tau%clvref='spherical'
* sn_rcv_tau%clvor='eastward-northward'
* sn_rcv_tau%clvgrd='T'
An alternative is to send only 2 stress components (x and y for for atm-ocn stress associated with jpr_otx1 and jpr_oty1). This is accomplished using the following combination of namelist parameters.
In this cased utau and utau_ice both get the x components and vtau and vtau_ice both get the y components of the stress fields that are passed to the ocean.
Note that the above set of namelist values and this one will result in identical wind stress assignments in NEMO when the AGCM only supplies a single x component and a single y component for wind stress (as is what happens currently in the AGCM).
* sn_rcv_tau%cldes='mixed oce-ice'
* sn_rcv_tau%clvref='spherical'
* sn_rcv_tau%clvor='eastward-northward'
* sn_rcv_tau%clvgrd='T'
The values assigned to the NEMO variables qsr_tot and qsr_ice in the subroutine sbc_cpl_ice_flx depend on the namelist parameter sn_rcv_qsr%cldes as follows
* if sn_rcv_qsr%cldes = 'oce only' then
* qsr_tot = O_QsrOce
* qsr_ice (is not used ??)
* if sn_rcv_qsr%cldes = 'conservative' then
* qsr_tot = O_QsrMix
* qsr_ice = O_QsrIce
* if sn_rcv_qsr%cldes = 'oce and ice' then
* qsr_tot = (1-fr_i)*O_QsrOce + fr_i*O_QsrIce
* qsr_ice = O_QsrIce
* if sn_rcv_qsr%cldes = 'mixed oce-ice' then
* qsr_tot = O_QsrMix
* qsr_ice = O_QsrMix*(1-palbi)/( 1-( albedo_oce_mix*(1-fr_i) + palbi*fr_i))
* palbi is ice albedo and albedo_oce_mix is ocean albedo (mixed clear/overcast)
The values assigned to the NEMO variables qns_tot and qns_ice in the subroutine sbc_cpl_ice_flx depend on the namelist parameter sn_rcv_qns%cldes as follows
* if sn_rcv_qns%cldes = 'oce only' then
* qns_tot = O_QnsOce
* qns_ice (is not used ??)
* if sn_rcv_qns%cldes = 'conservative' then
* qns_tot = O_QnsMix
* qns_ice = O_QnsIce
* if sn_rcv_qns%cldes = 'oce and ice' then
* qns_tot = (1-fr_i)*O_QnsOce + fr_i*O_QnsIce
* qns_ice = O_QnsIce
* if sn_rcv_qns%cldes = 'mixed oce-ice' then
* qns_tot = O_QnsMix
* qns_ice = O_dQnsdT*( pist - ((273.15 + psst)*(1-fr_i)+pist*fr_i) )
* psst is sea surface temperature (in Celcius) and pist is ice surface temperature (in Kelvin)
* In all cases the latent heat of melting snow is added so that qns_tot = qns_tot - (1-fr_i) * sprecip * lfus
The values assigned to the NEMO variables emp_tot (total evaporation - precipitation(liquid and solid) -runoff), emp_ice (ice sublimation - solid precipitation over the ice) and sprecip (solid precipitation over the ocean) in the subroutine sbc_cpl_ice_flx depend on the namelist parameter sn_rcv_emp%cldes as follows
* if sn_rcv_emp%cldes = 'conservative' then
* emp_tot = !OTotEvap - !OTotRain - !OTotSnow
* emp_ice = !OIceEvap - !OTotSnow
* sprecip = !OTotSnow
* if sn_rcv_emp%cldes = 'oce and ice' then
* emp_tot = (1-fr_i) * !OOEvaMPr + fr_i * !OSubMPre
* emp_ice = !OISubMSn
* sprecip = !OIceEvap - !OISubMSn
* in both cases runoff is removed so that emp_tot = emp_tot - O_Runoff
In terms of AGCM variables we have the following definitions. Here ℓ<sub>s</sub> = 2.834e6 J kg-1 is the latent heat of sublimation for fresh ice or snow.
%NAVY%There are 2 definitions for !OTotEvap shown. The red version is what appears in the old tables and is what is currently used in the coupler, but when this definition is used to calculate emp_tot below there appears to be an inconsistency between the 'conservative' and 'oce and ice' options. If the green definition is used then this inconsistency goes away. Why does this inconsistency exist? Does it impact our understanding of what NEMO thinks !OTotEvap actually is?%ENDCOLOR%
| !OTotRain | = | RAIN_atm | total liquid precipitation |
| !OTotSnow | = | SNOW_atm | total solid precipitation |
| !OTotEvap | = | %RED%HFLI_atm/ℓ<sub>s</sub> + RAIN_atm + SNOW_atm - BWGO_atm%ENDCOLOR% | total evaporation (over oce + ice sublimation) |
| !OTotEvap | = | %GREEN%HFLI_atm/ℓ<sub>s</sub> - BWGO_atm%ENDCOLOR% | total evaporation (over oce + ice sublimation) |
| !OIceEvap | = | HFLI_atm/ℓ<sub>s</sub> | evaporation over ice = sublimation |
| !OOEvaMPr | = | -BWGO_atm | ocean water budget = ocean Evap - ocean precip |
| !OSubMPre | = | HFLI_atm/ℓ<sub>s</sub> - RAIN_atm - SNOW_atm | sublimation - liquid precipitation - solid precipitation |
| !OISubMSn | = | HFLI_atm/ℓ<sub>s</sub> - SNOW_atm | ice solid water budget = sublimation - solid precipitation |
Rewriting the previous definitions for emp_tot, emp_ice and sprecip in terms of AGCM variables we have
* if sn_rcv_emp%cldes = 'conservative' then
* %RED%emp_tot = - BWGO_atm + HFLI_atm/ℓ<sub>s</sub>%ENDCOLOR%
* %GREEN%emp_tot = - BWGO_atm + HFLI_atm/ℓ<sub>s</sub> - RAIN_atm - SNOW_atm%ENDCOLOR%
* emp_ice = HFLI_atm/ℓ<sub>s</sub> - SNOW_atm
* sprecip = SNOW_atm
* if sn_rcv_emp%cldes = 'oce and ice' then
* emp_tot = -(1-fr_i) * BWGO_atm + fr_i * (HFLI_atm/ℓ<sub>s</sub> - RAIN_atm - SNOW_atm)
* emp_ice = HFLI_atm/ℓ<sub>s</sub> - SNOW_atm
* sprecip = SNOW_atm
* in both cases runoff is removed so that emp_tot = emp_tot - O_Runoff
%ENDTWISTY%
%ENDTWISTY%
%TWISTY{mode="div" showlink="Outstanding Questions" hidelink="<font color=#cf00ff>Hide Questions</font>" remember="on"}%
---+++ <center> Outstanding Questions </center>
1 Sea ice height and snow height over sea ice are delivered by LIM2. These heights, in meters, need to be converted to ice water equivalent (IWE) and snow water equivalent (SWE), in kg m-2,
when passed to the AGCM by multiplying by ice and snow densities. Can a 2D field for snow (and possibly ice) density be extracted from LIM2 to be used in the coupler to make this conversion?
If not, what is the proper value for a constant density to be used in this calculation?
1 The AGCM uses 273.16 for the conversion from Celcius to Kelvin but NEMO uses 273.15 for this conversion. Is this a concern?
1 Should fields sent from the ocean contain instantaneous values or values averaged over the coupling interval?
1 Should the run also include the WLOST updates (both lsmod and lssub) and the two radiation bugfixes, ie in gcm18_127c?
%ENDTWISTY%
-- Main.LarrySolheim - 2016-04-12