Changeset 7766 for branches/NERC
- Timestamp:
- 2017-03-07T19:30:03+01:00 (8 years ago)
- Location:
- branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM
- Files:
-
- 2 deleted
- 10 edited
-
CONFIG/SHARED/iodef_medusa_et_al.xml (modified) (8 diffs)
-
NEMO/TOP_SRC/MEDUSA/sms_medusa.F90 (modified) (1 diff)
-
NEMO/TOP_SRC/MEDUSA/trcbio_medusa.F90 (modified) (3 diffs)
-
NEMO/TOP_SRC/MEDUSA/trcdms_medusa.F90 (modified) (5 diffs)
-
NEMO/TOP_SRC/MEDUSA/trcdms_medusa.F90_COLIN4 (deleted)
-
NEMO/TOP_SRC/MEDUSA/trcnam_medusa.F90 (modified) (1 diff)
-
NEMO/TOP_SRC/MEDUSA/trcnam_medusa.F90_save (deleted)
-
NEMO/TOP_SRC/MEDUSA/trcopt_medusa.F90 (modified) (1 diff)
-
NEMO/TOP_SRC/MEDUSA/trcsed_medusa.F90 (modified) (5 diffs)
-
NEMO/TOP_SRC/MEDUSA/trcsms_medusa.F90 (modified) (3 diffs)
-
NEMO/TOP_SRC/trcrst.F90 (modified) (5 diffs)
-
NEMO/TOP_SRC/trcstp.F90 (modified) (1 diff)
Legend:
- Unmodified
- Added
- Removed
-
branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/CONFIG/SHARED/iodef_medusa_et_al.xml
r7516 r7766 1 1 <!-- GO6.0 version of the XIOS namelist --> 2 <!-- XIOS2 version namelist --> 3 2 <!-- XIOS version 2 --> 4 3 <?xml version="1.0"?> 4 5 5 <simulation> 6 6 7 <context id="nemo" time_origin="1950-01-01 00:00:00">7 <context id="nemo"> 8 8 9 9 <!-- $id$ --> … … 15 15 ============================================================================================================ 16 16 --> 17 <field_definition src="./field_def_bgc.xml"/> 17 18 <field_definition src="./field_def.xml"/> 18 <field_definition src="./field_def_bgc.xml"/>19 19 <!-- 20 20 ============================================================================================================ … … 24 24 ============================================================================================================ 25 25 --> 26 27 <file_definition type="multiple_file" name="@expname@_@freq@_@startdate@_@enddate@" sync_freq="1 0d" min_digits="4">26 27 <file_definition type="multiple_file" name="@expname@_@freq@_@startdate@_@enddate@" sync_freq="1mo" min_digits="4"> 28 28 29 29 <file_group id="1ts" output_freq="1ts" output_level="10" enabled=".TRUE."/> <!-- 1 time step files --> … … 33 33 <file_group id="3h" output_freq="3h" output_level="10" enabled=".TRUE."/> <!-- 3h files --> 34 34 <file_group id="4h" output_freq="4h" output_level="10" enabled=".TRUE."/> <!-- 4h files --> 35 <file_group id="6h" output_freq="6h" output_level="10" enabled=".FALSE."> <!-- 6h files --> 36 <file id="file8" name_suffix="_UK_shelf" description="Forcings for the UK shelf seas model" > 37 <field field_ref="toce" name="votemper" grid_ref="grid_T_3D_zoom" operation="instant" standard_name="sea_water_potential_temperature"/> 38 <field field_ref="soce" name="vosaline" grid_ref="grid_T_3D_zoom" operation="instant" standard_name="sea_water_salinity"/> 39 <field field_ref="ssh" name="sossheig" grid_ref="grid_T_2D_zoom" operation="instant" standard_name="sea_surface_height_above_geoid"/> 40 <field field_ref="uoce" name="vozocrtx" grid_ref="grid_U_3D_zoom" operation="instant" standard_name="sea_water_x_velocity"/> 41 <field field_ref="voce" name="vomecrty" grid_ref="grid_V_3D_zoom" operation="instant" standard_name="sea_water_y_velocity"/> 42 <field field_ref="ubar" name="ubar" grid_ref="grid_U_3D_zoom" operation="instant" standard_name="ocean_barotropic_current_along_i-axis"/> 43 <field field_ref="vbar" name="vbar" grid_ref="grid_V_3D_zoom" operation="instant" standard_name="ocean_barotropic_current_along_j-axis"/> 44 </file> 45 </file_group> 35 <file_group id="6h" output_freq="6h" output_level="10" enabled=".TRUE."/> <!-- 6h files --> 46 36 47 37 <file_group id="1d" output_freq="1d" output_level="10" enabled=".TRUE."/> <!-- 1d files --> 48 38 <file_group id="3d" output_freq="3d" output_level="10" enabled=".TRUE."/> <!-- 3d files --> 49 39 50 <file_group id="10d" output_freq="10d" output_level="10" enabled=".TRUE."> <!-- 10d files -->51 40 <file_group id="10d" output_freq="10d" output_level="10" enabled=".TRUE."> <!-- 10d files --> 41 52 42 <file id="file1" name_suffix="_grid_T" description="ocean T grid variables" > 53 <field_group group_ref="SPIN-UP_T" /> 43 <field field_ref="e3t" name="e3t" standard_name="cell_thickness" /> 44 <field field_ref="toce" name="votemper" operation="average" freq_op="1mo" > @toce_e3t / @e3t </field> 45 <field field_ref="toce" name="votemper2" operation="average" freq_op="1mo" > @toce2_e3t / @e3t </field> 46 <field field_ref="soce" name="vosaline" operation="average" freq_op="1mo" > @soce_e3t / @e3t </field> 47 <field field_ref="soce" name="vosaline2" operation="average" freq_op="1mo" > @soce2_e3t / @e3t </field> 48 <field field_ref="ttrd_totad" name="opottempadvect" operation="average" freq_op="1mo" > @ttrd_totad_e3t * 1026.0 * 3991.86795711963 </field> 49 <field field_ref="ttrd_iso" name="opottemppmdiff" operation="average" freq_op="1mo" > @ttrd_iso_e3t * 1026.0 * 3991.86795711963 </field> 50 <field field_ref="ttrd_zdfp" name="opottempdiff" operation="average" freq_op="1mo" > @ttrd_zdfp_e3t * 1026.0 * 3991.86795711963 </field> 51 <field field_ref="ttrd_evd" name="ttrd_evd" operation="average" freq_op="1mo" > @ttrd_evd_e3t / @e3t </field> 52 <field field_ref="ttrd_qns" name="ttrd_qns" operation="average" freq_op="1mo" > @ttrd_qns_e3t / @e3t_surf </field> 53 <field field_ref="ttrd_qsr" name="rsdoabsorb" operation="average" freq_op="1mo" > @ttrd_qsr_e3t * 1026.0 * 3991.86795711963 </field> 54 <field field_ref="ttrd_bbl" name="ttrd_bbl" operation="average" freq_op="1mo" > @ttrd_bbl_e3t / @e3t </field> 55 <field field_ref="ttrd_tot" name="opottemptend" operation="average" freq_op="1mo" > @ttrd_tot_e3t * 1026.0 * 3991.86795711963 </field> 56 <field field_ref="strd_totad" name="osaltadvect" operation="average" freq_op="1mo" > @strd_totad_e3t * 1026.0 * 0.001 </field> 57 <field field_ref="strd_iso" name="osaltpmdiff" operation="average" freq_op="1mo" > @strd_iso_e3t * 1026.0 * 0.001 </field> 58 <field field_ref="strd_zdfp" name="osaltdiff" operation="average" freq_op="1mo" > @strd_zdfp_e3t * 1026.0 * 0.001 </field> 59 <field field_ref="strd_evd" name="strd_evd" operation="average" freq_op="1mo" > @strd_evd_e3t / @e3t </field> 60 <field field_ref="strd_tot" name="osalttend" operation="average" freq_op="1mo" > @strd_tot_e3t * 1026.0 * 0.001 </field> 61 <field field_ref="strd_bbl" name="strd_bbl" operation="average" freq_op="1mo" > @strd_bbl_e3t / @e3t </field> 62 <field field_ref="ssh" name="sossheig" standard_name="sea_surface_height_above_geoid" /> 63 <field field_ref="ssh2" name="zossq" standard_name="square_of_sea_surface_height_above_geoid" /> 64 <field field_ref="empmr" name="sowaflup" standard_name="water_flux_out_of_sea_ice_and_sea_water" /> 65 <field field_ref="saltflx" name="sosafldo" standard_name="salt_flux_into_sea_water" /> 66 <field field_ref="qt" name="sohefldo" standard_name="surface_downward_heat_flux_in_sea_water" /> 67 <field field_ref="qsr" name="soshfldo" standard_name="downwelling_photosynthetic_radiative_flux_in_sea_water" /> 68 <field field_ref="mldkz5" name="somixhgt" standard_name="ocean_mixed_layer_thickness_defined_by_vertical_tracer_diffusivity" /> 69 <field field_ref="mldr10_1" name="somxl010" standard_name="ocean_mixed_layer_thickness_defined_by_sigma_theta" /> 70 <field field_ref="mldzint_1" name="sokaraml" standard_name="ocean_mixed_layer_thickness" /> 71 <field field_ref="mldzint_2" name="somlddbm" standard_name="ocean_mixed_layer_thickness" /> 72 <field field_ref="ice_cover" name="soicecov" standard_name="sea_ice_area_fraction" /> 73 <field field_ref="wspd" name="sowindsp" standard_name="wind_speed" /> 74 <field field_ref="qhcisf" name="sohflisf" standard_name="ice_shelf_heat_flux" /> 75 <field field_ref="qlatisf" name="soqlatisf" standard_name="ice_shelf_latent_heat_flux" /> 76 <field field_ref="fwfisf" name="sowflisf" standard_name="ice_shelf_melting" /> 77 <field field_ref="berg_total_melt" name="berg_total_melt" standard_name="" /> 78 <field field_ref="berg_total_heat_flux" name="berg_total_heat_flux" standard_name="" /> 79 <field field_ref="runoffs" name="sorunoff" standard_name="water_flux_into_sea_water_from_rivers" /> 80 <field field_ref="hflx_rain_cea" name="hfrainds" standard_name="temperature_flux_due_to_rainfall_expressed_as_heat_flux_into_sea_water"/> 81 <field field_ref="hflx_evap_cea" name="hfevapds" standard_name="temperature_flux_due_to_evaporation_expressed_as_heat_flux_out_of_sea_water"/> 82 <field field_ref="evap_ao_cea" name="evap_ao_cea" standard_name="water_evaporation_flux" /> 83 <field field_ref="rain" name="pr" standard_name="rainfall_flux" /> 84 <field field_ref="snowpre" name="snowpre" standard_name="snowfall_flux" /> 85 <field field_ref="snow_ao_cea" name="snow_ao_cea" standard_name="snowfall_flux_over_ocean" /> 86 <field field_ref="snow_ai_cea" name="snow_ai_cea" standard_name="snowfall_flux_over_sea_ice" /> 54 87 </file> 55 88 56 89 <file id="file2" name_suffix="_grid_U" description="ocean U grid variables" > 57 <field_group group_ref="SPIN-UP_U" /> 90 <field field_ref="e3u" name="e3u" standard_name="cell_thickness" /> 91 <field field_ref="uoce" name="vozocrtx" operation="average" freq_op="1mo" > @uoce_e3u / @e3u </field> 92 <field field_ref="uoce" name="vozocrtx2" operation="average" freq_op="1mo" > @uoce2_e3u / @e3u </field> 93 <field field_ref="ut" name="ut" operation="average" freq_op="1mo" > @ut_e3u / @e3u </field> 94 <field field_ref="us" name="us" operation="average" freq_op="1mo" > @us_e3u / @e3u </field> 95 <field field_ref="utau" name="sozotaux" standard_name="surface_downward_x_stress" /> 58 96 </file> 59 97 60 98 <file id="file3" name_suffix="_grid_V" description="ocean V grid variables" > 61 <field_group group_ref="SPIN-UP_V" /> 62 </file> 63 64 <file id="file4" name_suffix="_grid_W" description="ocean W grid variables" > 65 <field_group group_ref="SPIN-UP_W" /> 99 <field field_ref="e3v" name="e3v" standard_name="cell_thickness" /> 100 <field field_ref="voce" name="vomecrty" operation="average" freq_op="1mo" > @voce_e3v / @e3v </field> 101 <field field_ref="voce" name="vomecrty2" operation="average" freq_op="1mo" > @voce2_e3v / @e3v </field> 102 <field field_ref="vt" name="vt" operation="average" freq_op="1mo" > @vt_e3v / @e3v </field> 103 <field field_ref="vs" name="vs" operation="average" freq_op="1mo" > @vs_e3v / @e3v </field> 104 <field field_ref="vtau" name="sometauy" standard_name="surface_downward_y_stress" /> 66 105 </file> 67 106 68 <file id="file5" name_suffix="_diaptr" description="process diagnostics" > 69 <field_group group_ref="SPIN-UP_diaptr" /> 70 </file> 71 107 <file id="file4" name_suffix="_grid_W" description="ocean W grid variables" > 108 <field field_ref="woce" name="vovecrtz" standard_name="upward_sea_water_velocity" /> 109 <field field_ref="wt" name="wt" standard_name="" /> 110 <field field_ref="ws" name="ws" standard_name="" /> 111 <field field_ref="avt" name="votkeavt" standard_name="ocean_vertical_heat_diffusivity" /> 112 <field field_ref="avm" name="votkeavm" standard_name="ocean_vertical_momentum_diffusivity" /> 113 <field field_ref="aht2d" name="soleahtw" standard_name="vertical_component_of_ocean_xy_tracer_diffusivity" /> 114 </file> 115 116 <file id="file5" name_suffix="_diaptr" description="processed diagnostics" > 117 <field field_ref="zomsfglo" name="zomsfglo" standard_name="meridional_streamfunction_global" grid_ref="gznl_W_3D" /> 118 <field field_ref="zotemglo" name="zotemglo" standard_name="zonal_mean_temperature_global" grid_ref="gznl_T_3D" /> 119 <field field_ref="zosrfglo" name="zosrfglo" standard_name="zonal_mean_surface_global" grid_ref="gznl_T_3D" /> 120 <field field_ref="zosalglo" name="zosalglo" standard_name="zonal_mean_salinity_global" grid_ref="gznl_T_3D" /> 121 <field field_ref="zomsfatl" name="zomsfatl" standard_name="meridional_streamfunction_atlantic" grid_ref="gznl_W_3D" /> 122 <field field_ref="zotematl" name="zotematl" standard_name="zonal_mean_temperature_atlantic" grid_ref="gznl_T_3D" /> 123 <field field_ref="zosrfatl" name="zosrfatl" standard_name="zonal_mean_surface_atlantic" grid_ref="gznl_T_3D" /> 124 <field field_ref="zosalatl" name="zosalatl" standard_name="zonal_mean_salinity_atlantic" grid_ref="gznl_T_3D" /> 125 <field field_ref="zomsfpac" name="zomsfpac" standard_name="meridional_streamfunction_pacific" grid_ref="gznl_W_3D" /> 126 <field field_ref="zotempac" name="zotempac" standard_name="zonal_mean_temperature_pacific" grid_ref="gznl_T_3D" /> 127 <field field_ref="zosrfpac" name="zosrfpac" standard_name="zonal_mean_surface_pacific" grid_ref="gznl_T_3D" /> 128 <field field_ref="zosalpac" name="zosalpac" standard_name="zonal_mean_salinity_pacific" grid_ref="gznl_T_3D" /> 129 <field field_ref="zomsfind" name="zomsfind" standard_name="meridional_streamfunction_indian" grid_ref="gznl_W_3D" /> 130 <field field_ref="zotemind" name="zotemind" standard_name="zonal_mean_temperature_indian" grid_ref="gznl_T_3D" /> 131 <field field_ref="zosrfind" name="zosrfind" standard_name="zonal_mean_surface_indian" grid_ref="gznl_T_3D" /> 132 <field field_ref="zosalind" name="zosalind" standard_name="zonal_mean_salinity_indian" grid_ref="gznl_T_3D" /> 133 <field field_ref="zomsfipc" name="zomsfipc" standard_name="meridional_streamfunction_indian_pacific" grid_ref="gznl_W_3D" /> 134 <field field_ref="zotemipc" name="zotemipc" standard_name="zonal_mean_temperature_indian_pacific" grid_ref="gznl_T_3D" /> 135 <field field_ref="zosrfipc" name="zosrfipc" standard_name="zonal_mean_surface_indian_pacific" grid_ref="gznl_T_3D" /> 136 <field field_ref="zosalipc" name="zosalipc" standard_name="zonal_mean_salinity_indian_pacific" grid_ref="gznl_T_3D" /> 137 <field field_ref="sophtadv" name="sophtadv" standard_name="advective_heat_transport" grid_ref="gznl_T_2D" /> 138 <field field_ref="sopstadv" name="sopstadv" standard_name="advective_salt_transport" grid_ref="gznl_T_2D" /> 139 <field field_ref="sophtldf" name="sophtldf" standard_name="diffusive_heat_transport" grid_ref="gznl_T_2D" /> 140 <field field_ref="sopstldf" name="sopstldf" standard_name="diffusive_salt_transport" grid_ref="gznl_T_2D" /> 141 <field field_ref="sophtadv_atl" name="sophtadv_atl" standard_name="advective_heat_transport_atlantic" grid_ref="gznl_T_2D" /> 142 <field field_ref="sopstadv_atl" name="sopstadv_atl" standard_name="advective_salt_transport_atlantic" grid_ref="gznl_T_2D" /> 143 <field field_ref="sophtldf_atl" name="sophtldf_atl" standard_name="diffusive_heat_transport_atlantic" grid_ref="gznl_T_2D" /> 144 <field field_ref="sopstldf_atl" name="sopstldf_atl" standard_name="diffusive_salt_transport_atlantic" grid_ref="gznl_T_2D" /> 145 <field field_ref="sophtadv_pac" name="sophtadv_pac" standard_name="advective_heat_transport_pacific" grid_ref="gznl_T_2D" /> 146 <field field_ref="sopstadv_pac" name="sopstadv_pac" standard_name="advective_salt_transport_pacific" grid_ref="gznl_T_2D" /> 147 <field field_ref="sophtldf_pac" name="sophtldf_pac" standard_name="diffusive_heat_transport_pacific" grid_ref="gznl_T_2D" /> 148 <field field_ref="sopstldf_pac" name="sopstldf_pac" standard_name="diffusive_salt_transport_pacific" grid_ref="gznl_T_2D" /> 149 <field field_ref="sophtadv_ind" name="sophtadv_ind" standard_name="advective_heat_transport_indian" grid_ref="gznl_T_2D" /> 150 <field field_ref="sopstadv_ind" name="sopstadv_ind" standard_name="advective_salt_transport_indian" grid_ref="gznl_T_2D" /> 151 <field field_ref="sophtldf_ind" name="sophtldf_ind" standard_name="diffusive_heat_transport_indian" grid_ref="gznl_T_2D" /> 152 <field field_ref="sopstldf_ind" name="sopstldf_ind" standard_name="diffusive_salt_transport_indian" grid_ref="gznl_T_2D" /> 153 <field field_ref="sophtadv_ipc" name="sophtadv_ipc" standard_name="advective_heat_transport_indopacific" grid_ref="gznl_T_2D" /> 154 <field field_ref="sopstadv_ipc" name="sopstadv_ipc" standard_name="advective_salt_transport_indopacific" grid_ref="gznl_T_2D" /> 155 <field field_ref="sophtldf_ipc" name="sophtldf_ipc" standard_name="diffusive_heat_transport_indopacific" grid_ref="gznl_T_2D" /> 156 <field field_ref="sopstldf_ipc" name="sopstldf_ipc" standard_name="diffusive_salt_transport_indopacific" grid_ref="gznl_T_2D" /> 157 </file> 158 72 159 <file id="file6" name_suffix="_ptrc_T" description="Medusa sms variables" > 73 160 <field_group group_ref="groupMEDUSA" /> … … 89 176 90 177 91 < !-- file_group id="1m" output_freq="1mo" split_freq="1mo" output_level="10" enabled=".TRUE." --><!-- 1mo files -->92 < file_group id="1m" output_freq="1mo" output_level="10" enabled=".FALSE."> <!-- 1mo files -->178 <file_group id="1m" output_freq="1mo" split_freq="1mo" output_level="10" enabled=".FALSE."> <!-- 1mo files --> 179 <!-- <file_group id="1m" output_freq="1mo" output_level="10" enabled=".TRUE."> - 1mo files --> 93 180 94 181 <file id="file9" name_suffix="_grid_T" description="ocean T grid variables" > 95 <field_group group_ref="SPIN-UP_T" /> 96 </file> 97 182 <field field_ref="e3t" name="e3t" standard_name="cell_thickness" /> 183 <field field_ref="toce" name="votemper" operation="average" freq_op="1mo" > @toce_e3t / @e3t </field> 184 <field field_ref="toce" name="votemper2" operation="average" freq_op="1mo" > @toce2_e3t / @e3t </field> 185 <field field_ref="soce" name="vosaline" operation="average" freq_op="1mo" > @soce_e3t / @e3t </field> 186 <field field_ref="soce" name="vosaline2" operation="average" freq_op="1mo" > @soce2_e3t / @e3t </field> 187 <field field_ref="ttrd_totad" name="opottempadvect" operation="average" freq_op="1mo" > @ttrd_totad_e3t * 1026.0 * 3991.86795711963 </field> 188 <field field_ref="ttrd_iso" name="opottemppmdiff" operation="average" freq_op="1mo" > @ttrd_iso_e3t * 1026.0 * 3991.86795711963 </field> 189 <field field_ref="ttrd_zdfp" name="opottempdiff" operation="average" freq_op="1mo" > @ttrd_zdfp_e3t * 1026.0 * 3991.86795711963 </field> 190 <field field_ref="ttrd_evd" name="ttrd_evd" operation="average" freq_op="1mo" > @ttrd_evd_e3t / @e3t </field> 191 <field field_ref="ttrd_qns" name="ttrd_qns" operation="average" freq_op="1mo" > @ttrd_qns_e3t / @e3t_surf </field> 192 <field field_ref="ttrd_qsr" name="rsdoabsorb" operation="average" freq_op="1mo" > @ttrd_qsr_e3t * 1026.0 * 3991.86795711963 </field> 193 <field field_ref="ttrd_bbl" name="ttrd_bbl" operation="average" freq_op="1mo" > @ttrd_bbl_e3t / @e3t </field> 194 <field field_ref="ttrd_tot" name="opottemptend" operation="average" freq_op="1mo" > @ttrd_tot_e3t * 1026.0 * 3991.86795711963 </field> 195 <field field_ref="strd_totad" name="osaltadvect" operation="average" freq_op="1mo" > @strd_totad_e3t * 1026.0 * 0.001 </field> 196 <field field_ref="strd_iso" name="osaltpmdiff" operation="average" freq_op="1mo" > @strd_iso_e3t * 1026.0 * 0.001 </field> 197 <field field_ref="strd_zdfp" name="osaltdiff" operation="average" freq_op="1mo" > @strd_zdfp_e3t * 1026.0 * 0.001 </field> 198 <field field_ref="strd_evd" name="strd_evd" operation="average" freq_op="1mo" > @strd_evd_e3t / @e3t </field> 199 <field field_ref="strd_tot" name="osalttend" operation="average" freq_op="1mo" > @strd_tot_e3t * 1026.0 * 0.001 </field> 200 <field field_ref="strd_bbl" name="strd_bbl" operation="average" freq_op="1mo" > @strd_bbl_e3t / @e3t </field> 201 <field field_ref="ssh" name="sossheig" standard_name="sea_surface_height_above_geoid" /> 202 <field field_ref="ssh2" name="zossq" standard_name="square_of_sea_surface_height_above_geoid" /> 203 <field field_ref="empmr" name="sowaflup" standard_name="water_flux_out_of_sea_ice_and_sea_water" /> 204 <field field_ref="saltflx" name="sosafldo" standard_name="salt_flux_into_sea_water" /> 205 <field field_ref="qt" name="sohefldo" standard_name="surface_downward_heat_flux_in_sea_water" /> 206 <field field_ref="qsr" name="soshfldo" standard_name="downwelling_photosynthetic_radiative_flux_in_sea_water" /> 207 <field field_ref="mldkz5" name="somixhgt" standard_name="ocean_mixed_layer_thickness_defined_by_vertical_tracer_diffusivity" /> 208 <field field_ref="mldr10_1" name="somxl010" standard_name="ocean_mixed_layer_thickness_defined_by_sigma_theta" /> 209 <field field_ref="mldzint_1" name="sokaraml" standard_name="ocean_mixed_layer_thickness" /> 210 <field field_ref="mldzint_2" name="somlddbm" standard_name="ocean_mixed_layer_thickness" /> 211 <field field_ref="ice_cover" name="soicecov" standard_name="sea_ice_area_fraction" /> 212 <field field_ref="wspd" name="sowindsp" standard_name="wind_speed" /> 213 <field field_ref="qhcisf" name="sohflisf" standard_name="ice_shelf_heat_flux" /> 214 <field field_ref="qlatisf" name="soqlatisf" standard_name="ice_shelf_latent_heat_flux" /> 215 <field field_ref="fwfisf" name="sowflisf" standard_name="ice_shelf_melting" /> 216 <field field_ref="berg_total_melt" name="berg_total_melt" standard_name="" /> 217 <field field_ref="berg_total_heat_flux" name="berg_total_heat_flux" standard_name="" /> 218 <field field_ref="runoffs" name="sorunoff" standard_name="water_flux_into_sea_water_from_rivers" /> 219 <field field_ref="hflx_rain_cea" name="hfrainds" standard_name="temperature_flux_due_to_rainfall_expressed_as_heat_flux_into_sea_water"/> 220 <field field_ref="hflx_evap_cea" name="hfevapds" standard_name="temperature_flux_due_to_evaporation_expressed_as_heat_flux_out_of_sea_water"/> 221 <field field_ref="evap_ao_cea" name="evap_ao_cea" standard_name="water_evaporation_flux" /> 222 <field field_ref="rain" name="pr" standard_name="rainfall_flux" /> 223 <field field_ref="snowpre" name="snowpre" standard_name="snowfall_flux" /> 224 <field field_ref="snow_ao_cea" name="snow_ao_cea" standard_name="snowfall_flux_over_ocean" /> 225 <field field_ref="snow_ai_cea" name="snow_ai_cea" standard_name="snowfall_flux_over_sea_ice" /> 226 </file> 227 98 228 <file id="file10" name_suffix="_grid_U" description="ocean U grid variables" > 99 <field_group group_ref="SPIN-UP_U" /> 100 </file> 101 229 <field field_ref="e3u" name="e3u" standard_name="cell_thickness" /> 230 <field field_ref="uoce" name="vozocrtx" operation="average" freq_op="1mo" > @uoce_e3u / @e3u </field> 231 <field field_ref="uoce" name="vozocrtx2" operation="average" freq_op="1mo" > @uoce2_e3u / @e3u </field> 232 <field field_ref="ut" name="ut" operation="average" freq_op="1mo" > @ut_e3u / @e3u </field> 233 <field field_ref="us" name="us" operation="average" freq_op="1mo" > @us_e3u / @e3u </field> 234 <field field_ref="utau" name="sozotaux" standard_name="surface_downward_x_stress" /> 235 </file> 236 102 237 <file id="file11" name_suffix="_grid_V" description="ocean V grid variables" > 103 <field_group group_ref="SPIN-UP_V" /> 104 </file> 105 238 <field field_ref="e3v" name="e3v" standard_name="cell_thickness" /> 239 <field field_ref="voce" name="vomecrty" operation="average" freq_op="1mo" > @voce_e3v / @e3v </field> 240 <field field_ref="voce" name="vomecrty2" operation="average" freq_op="1mo" > @voce2_e3v / @e3v </field> 241 <field field_ref="vt" name="vt" operation="average" freq_op="1mo" > @vt_e3v / @e3v </field> 242 <field field_ref="vs" name="vs" operation="average" freq_op="1mo" > @vs_e3v / @e3v </field> 243 <field field_ref="vtau" name="sometauy" standard_name="surface_downward_y_stress" /> 244 </file> 245 106 246 <file id="file12" name_suffix="_grid_W" description="ocean W grid variables" > 107 <field_group group_ref="SPIN-UP_W" /> 108 </file> 109 247 <field field_ref="woce" name="vovecrtz" standard_name="upward_sea_water_velocity" /> 248 <field field_ref="wt" name="wt" standard_name="" /> 249 <field field_ref="ws" name="ws" standard_name="" /> 250 <field field_ref="avt" name="votkeavt" standard_name="ocean_vertical_heat_diffusivity" /> 251 <field field_ref="avm" name="votkeavm" standard_name="ocean_vertical_momentum_diffusivity" /> 252 <field field_ref="aht2d" name="soleahtw" standard_name="vertical_component_of_ocean_xy_tracer_diffusivity" /> 253 </file> 254 110 255 <file id="file13" name_suffix="_grid_W" description="ocean W grid variables" > 111 <field_group group_ref="SPIN-UP_diaptr" /> 112 </file> 113 256 <field field_ref="zomsfglo" name="zomsfglo" standard_name="meridional_streamfunction_global" grid_ref="gznl_W_3D" /> 257 <field field_ref="zotemglo" name="zotemglo" standard_name="zonal_mean_temperature_global" grid_ref="gznl_T_3D" /> 258 <field field_ref="zosrfglo" name="zosrfglo" standard_name="zonal_mean_surface_global" grid_ref="gznl_T_3D" /> 259 <field field_ref="zosalglo" name="zosalglo" standard_name="zonal_mean_salinity_global" grid_ref="gznl_T_3D" /> 260 <field field_ref="zomsfatl" name="zomsfatl" standard_name="meridional_streamfunction_atlantic" grid_ref="gznl_W_3D" /> 261 <field field_ref="zotematl" name="zotematl" standard_name="zonal_mean_temperature_atlantic" grid_ref="gznl_T_3D" /> 262 <field field_ref="zosrfatl" name="zosrfatl" standard_name="zonal_mean_surface_atlantic" grid_ref="gznl_T_3D" /> 263 <field field_ref="zosalatl" name="zosalatl" standard_name="zonal_mean_salinity_atlantic" grid_ref="gznl_T_3D" /> 264 <field field_ref="zomsfpac" name="zomsfpac" standard_name="meridional_streamfunction_pacific" grid_ref="gznl_W_3D" /> 265 <field field_ref="zotempac" name="zotempac" standard_name="zonal_mean_temperature_pacific" grid_ref="gznl_T_3D" /> 266 <field field_ref="zosrfpac" name="zosrfpac" standard_name="zonal_mean_surface_pacific" grid_ref="gznl_T_3D" /> 267 <field field_ref="zosalpac" name="zosalpac" standard_name="zonal_mean_salinity_pacific" grid_ref="gznl_T_3D" /> 268 <field field_ref="zomsfind" name="zomsfind" standard_name="meridional_streamfunction_indian" grid_ref="gznl_W_3D" /> 269 <field field_ref="zotemind" name="zotemind" standard_name="zonal_mean_temperature_indian" grid_ref="gznl_T_3D" /> 270 <field field_ref="zosrfind" name="zosrfind" standard_name="zonal_mean_surface_indian" grid_ref="gznl_T_3D" /> 271 <field field_ref="zosalind" name="zosalind" standard_name="zonal_mean_salinity_indian" grid_ref="gznl_T_3D" /> 272 <field field_ref="zomsfipc" name="zomsfipc" standard_name="meridional_streamfunction_indian_pacific" grid_ref="gznl_W_3D" /> 273 <field field_ref="zotemipc" name="zotemipc" standard_name="zonal_mean_temperature_indian_pacific" grid_ref="gznl_T_3D" /> 274 <field field_ref="zosrfipc" name="zosrfipc" standard_name="zonal_mean_surface_indian_pacific" grid_ref="gznl_T_3D" /> 275 <field field_ref="zosalipc" name="zosalipc" standard_name="zonal_mean_salinity_indian_pacific" grid_ref="gznl_T_3D" /> 276 <field field_ref="sophtadv" name="sophtadv" standard_name="advective_heat_transport" grid_ref="gznl_T_2D" /> 277 <field field_ref="sopstadv" name="sopstadv" standard_name="advective_salt_transport" grid_ref="gznl_T_2D" /> 278 <field field_ref="sophtldf" name="sophtldf" standard_name="diffusive_heat_transport" grid_ref="gznl_T_2D" /> 279 <field field_ref="sopstldf" name="sopstldf" standard_name="diffusive_salt_transport" grid_ref="gznl_T_2D" /> 280 <field field_ref="sophtadv_atl" name="sophtadv_atl" standard_name="advective_heat_transport_atlantic" grid_ref="gznl_T_2D" /> 281 <field field_ref="sopstadv_atl" name="sopstadv_atl" standard_name="advective_salt_transport_atlantic" grid_ref="gznl_T_2D" /> 282 <field field_ref="sophtldf_atl" name="sophtldf_atl" standard_name="diffusive_heat_transport_atlantic" grid_ref="gznl_T_2D" /> 283 <field field_ref="sopstldf_atl" name="sopstldf_atl" standard_name="diffusive_salt_transport_atlantic" grid_ref="gznl_T_2D" /> 284 <field field_ref="sophtadv_pac" name="sophtadv_pac" standard_name="advective_heat_transport_pacific" grid_ref="gznl_T_2D" /> 285 <field field_ref="sopstadv_pac" name="sopstadv_pac" standard_name="advective_salt_transport_pacific" grid_ref="gznl_T_2D" /> 286 <field field_ref="sophtldf_pac" name="sophtldf_pac" standard_name="diffusive_heat_transport_pacific" grid_ref="gznl_T_2D" /> 287 <field field_ref="sopstldf_pac" name="sopstldf_pac" standard_name="diffusive_salt_transport_pacific" grid_ref="gznl_T_2D" /> 288 <field field_ref="sophtadv_ind" name="sophtadv_ind" standard_name="advective_heat_transport_indian" grid_ref="gznl_T_2D" /> 289 <field field_ref="sopstadv_ind" name="sopstadv_ind" standard_name="advective_salt_transport_indian" grid_ref="gznl_T_2D" /> 290 <field field_ref="sophtldf_ind" name="sophtldf_ind" standard_name="diffusive_heat_transport_indian" grid_ref="gznl_T_2D" /> 291 <field field_ref="sopstldf_ind" name="sopstldf_ind" standard_name="diffusive_salt_transport_indian" grid_ref="gznl_T_2D" /> 292 <field field_ref="sophtadv_ipc" name="sophtadv_ipc" standard_name="advective_heat_transport_indopacific" grid_ref="gznl_T_2D" /> 293 <field field_ref="sopstadv_ipc" name="sopstadv_ipc" standard_name="advective_salt_transport_indopacific" grid_ref="gznl_T_2D" /> 294 <field field_ref="sophtldf_ipc" name="sophtldf_ipc" standard_name="diffusive_heat_transport_indopacific" grid_ref="gznl_T_2D" /> 295 <field field_ref="sopstldf_ipc" name="sopstldf_ipc" standard_name="diffusive_salt_transport_indopacific" grid_ref="gznl_T_2D" /> 296 </file> 297 114 298 <file id="file14" name_suffix="_ptrc_T" description="Medusa sms variables" > 115 299 <field_group group_ref="groupMEDUSA" /> … … 147 331 ============================================================================================================ 148 332 --> 149 333 150 334 <axis_definition> 151 335 <axis id="deptht" standard_name="depth" long_name="Vertical T levels" unit="m" positive="down" /> … … 155 339 <axis id="nfloat" long_name="Float number" unit="-" /> 156 340 <axis id="icbcla" long_name="Iceberg class" unit="-" /> 341 <axis id="ncatice" long_name="Ice category" unit="1" /> 342 <axis id="iax_20C" long_name="20 degC isotherm" unit="degC" /> 343 <axis id="iax_28C" long_name="28 degC isotherm" unit="degC" /> 344 <axis id="deptht_surface" axis_ref="deptht" > 345 <zoom_axis begin=" 1 " n=" 1 " /> 346 </axis> 347 157 348 </axis_definition> 158 349 159 350 <domain_definition src="./domain_def.xml"/> 160 351 161 <grid_definition> 162 <grid id="grid_T_2D" domain_ref="grid_T"/> 163 <grid id="grid_T_3D" domain_ref="grid_T" axis_ref="deptht"/> 164 <grid id="grid_U_2D" domain_ref="grid_U"/> 165 <grid id="grid_U_3D" domain_ref="grid_U" axis_ref="depthu"/> 166 <grid id="grid_V_2D" domain_ref="grid_V"/> 167 <grid id="grid_V_3D" domain_ref="grid_V" axis_ref="depthv"/> 168 <grid id="grid_W_2D" domain_ref="grid_W"/> 169 <grid id="grid_W_3D" domain_ref="grid_W" axis_ref="depthw"/> 170 <grid id="grid_T_2D_zoom" domain_ref="zoomT"/> 171 <grid id="grid_T_3D_zoom" domain_ref="zoomT" axis_ref="deptht"/> 172 <grid id="grid_U_3D_zoom" domain_ref="zoomU" axis_ref="depthu"/> 173 <grid id="grid_V_3D_zoom" domain_ref="zoomV" axis_ref="depthv"/> 174 <grid id="gznl_T_2D" domain_ref="gznl"/> 175 <grid id="gznl_T_3D" domain_ref="gznl" axis_ref="deptht"/> 176 <grid id="gznl_W_3D" domain_ref="gznl" axis_ref="depthw"/> 177 </grid_definition> 352 <grid_definition> 353 <!-- --> 354 <grid id="grid_T_2D" > 355 <domain id="grid_T" /> 356 </grid> 357 <!-- --> 358 <grid id="grid_T_3D_ncatice" > 359 <domain id="grid_T" /> 360 <axis id="ncatice" /> 361 </grid> 362 <!-- --> 363 <grid id="grid_T_3D" > 364 <domain id="grid_T" /> 365 <axis id="deptht" /> 366 </grid> 367 <!-- --> 368 <grid id="grid_U_2D" > 369 <domain id="grid_U" /> 370 </grid> 371 <!-- --> 372 <grid id="grid_U_3D" > 373 <domain id="grid_U" /> 374 <axis id="depthu" /> 375 </grid> 376 <!-- --> 377 <grid id="grid_V_2D" > 378 <domain id="grid_V" /> 379 </grid> 380 <!-- --> 381 <grid id="grid_V_3D" > 382 <domain id="grid_V" /> 383 <axis id="depthv" /> 384 </grid> 385 <!-- --> 386 <grid id="grid_W_2D" > 387 <domain id="grid_W" /> 388 </grid> 389 <!-- --> 390 <grid id="grid_W_3D" > 391 <domain id="grid_W" /> 392 <axis id="depthw" /> 393 </grid> 394 <!-- --> 395 <grid id="grid_1point" > 396 <domain id="1point"/> 397 </grid> 398 <!-- --> 399 <grid id="grid_T_nfloat" > 400 <domain id="grid_T" /> 401 <axis id="nfloat" /> 402 </grid> 403 <!-- --> 404 <grid id="grid_EqT" > 405 <domain id="EqT" /> 406 </grid> 407 <!-- --> 408 <grid id="gznl_T_2D"> 409 <domain id="ptr" /> 410 </grid> 411 <!-- --> 412 <grid id="gznl_T_3D"> 413 <domain id="ptr" /> 414 <axis id="deptht" /> 415 </grid> 416 <!-- --> 417 <grid id="gznl_W_2D"> 418 <domain id="ptr" /> 419 </grid> 420 <!-- --> 421 <grid id="gznl_W_3D"> 422 <domain id="ptr" /> 423 <axis id="depthw" /> 424 </grid> 425 <!-- --> 426 <grid id="grid_T_surface_extract"> 427 <domain domain_ref="grid_T" /> 428 <axis axis_ref="deptht_surface" /> 429 </grid> 430 431 </grid_definition> 178 432 </context> 179 433 … … 183 437 <variable_definition> 184 438 185 <!-- 186 We must have buffer_size > jpi*jpj*jpk*8 (with jpi and jpj the subdomain size) 187 --> 188 <variable id="buffer_size" type="integer">30000000</variable> 189 <variable id="buffer_server_factor_size" type="integer">2</variable> 190 <variable id="info_level" type="integer">0</variable> 191 <variable id="using_server" type="boolean">true</variable> 192 <variable id="using_oasis" type="boolean">true</variable> 193 <variable id="oasis_codes_id" type="string" >toyoce</variable> 439 <variable id="info_level" type="int">0</variable> 440 <variable id="using_server" type="bool">false</variable> 441 <variable id="using_oasis" type="bool">true</variable> 442 <variable id="oasis_codes_id" type="string" >toyoce</variable> 194 443 195 444 </variable_definition> -
branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/MEDUSA/sms_medusa.F90
r6810 r7766 445 445 !!---------------------------------------------------------------------- 446 446 USE lib_mpp , ONLY: ctl_warn 447 INTEGER :: ierr( 7) ! Local variables447 INTEGER :: ierr(8) ! Local variables 448 448 !!---------------------------------------------------------------------- 449 449 ierr(:) = 0 -
branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcbio_medusa.F90
r7703 r7766 533 533 foxy_cons(:,:) = 0.0 !! oxygen consumption 534 534 foxy_anox(:,:) = 0.0 !! unrealised oxygen consumption 535 !! 535 536 # endif 536 537 ftot_n(:,:) = 0.0 !! N inventory … … 571 572 f_riv_c(:,:) = 0.0 !! riverine C input 572 573 f_riv_alk(:,:) = 0.0 !! riverine alk input 574 !! 575 !! Jpalm -- 06-03-2017 -- Forgotten var to init 576 f_omarg(:,:) = 0.0 !! 577 f_omcal(:,:) = 0.0 578 xFree(:,:) = 0.0 !! state variables for iron-ligand system 579 fcomm_resp(:,:) = 0.0 580 fprn_ml(:,:) = 0.0 !! mixed layer PP diagnostics 581 fprd_ml(:,:) = 0.0 !! mixed layer PP diagnostics 582 573 583 !! 574 584 !! allocate and initiate 2D diag … … 724 734 ENDIF 725 735 !! 726 !! skip FDT_XXXX, RG_XXXXF, FDS_XXXX, RGS XXXXF diagnostics here736 !! skip FDT_XXXX, RG_XXXXF, FDS_XXXX, RGS_XXXXF diagnostics here 727 737 !! 728 738 IF( med_diag%REMINN%dgsave ) THEN -
branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcdms_medusa.F90
r6719 r7766 5 5 !!====================================================================== 6 6 !! History : 7 !! - ! 2014-08 (J. Palmi éri - A. Yool) added for UKESM1 project7 !! - ! 2014-08 (J. Palmieri - A. Yool) added for UKESM1 project 8 8 !!---------------------------------------------------------------------- 9 9 #if defined key_medusa && defined key_roam … … 44 44 !! 45 45 !! Title : Calculates DMS ocean surface concentration 46 !! Author : Julien Palmi éri and Andrew Yool46 !! Author : Julien Palmieri and Andrew Yool 47 47 !! Date : 08/08/14 48 48 !! … … 118 118 ! 119 119 ! AXY (13/03/15): Simo & Dachs (2002) 120 fq1 = (-1 * log(mld)) + 5.7120 fq1 = (-1.0 * log(mld)) + 5.7 121 121 fq2 = (55.8 * cmr) + 0.6 122 122 if (cmr < 0.02) then … … 136 136 ! 137 137 ! AXY (13/03/15): Halloran et al. (2010) 138 fq1 = (-1 * log(mld)) + 5.7138 fq1 = (-1.0 * log(mld)) + 5.7 139 139 fq2 = (55.8 * cmr) + 0.6 140 140 fq3 = (90.0 / mld) … … 179 179 180 180 181 -
branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcnam_medusa.F90
r7709 r7766 84 84 & xburial, & 85 85 & jriver_n,jriver_si,jriver_c,jriver_alk,jriver_dep, & 86 & friver_dep, &87 86 & xsdiss, & 88 & vsed,xhr, &89 87 & sedlam,sedlostpoc,jpkb,jdms,jdms_input,jdms_model 90 88 #if defined key_roam 91 89 NAMELIST/natroam/ xthetaphy,xthetazoo,xthetanit, & 92 & xthetarem,xo2min, & 93 & f3_pH,f3_h2co3,f3_hco3,f3_co3,f3_omcal,f3_omarg, & 94 & f2_ccd_cal,f2_ccd_arg 90 & xthetarem,xo2min 95 91 #endif 96 92 NAMELIST/natopt/xkg0,xkr0,xkgp,xkrp,xlg,xlr,rpig -
branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcopt_medusa.F90
r5841 r7766 156 156 157 157 ! Euphotic layer depth 158 !! Jpalm -- 06-03-2017 -- add init xze, to avoid halo problems within the 159 !! writing process 160 xze(:,:) = 0.0 158 161 DO jj = 1, jpj 159 162 DO ji = 1, jpi -
branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcsed_medusa.F90
r6844 r7766 41 41 42 42 PUBLIC trc_sed_medusa ! called in ??? 43 PUBLIC trc_sed_medusa_sbc 43 PUBLIC trc_sed_medusa_sbc 44 PUBLIC trc_sed_medusa_dust 44 45 45 46 !! * Module variables … … 126 127 127 128 !! AXY (20/11/14): alter this to report on first MEDUSA call 128 !! IF( kt == nit000 ) THEN129 129 IF( kt == nittrc000 ) THEN 130 130 IF(lwp) WRITE(numout,*) … … 133 133 IF(lwp) WRITE(numout,*) ' kt =',kt 134 134 ENDIF 135 136 !! AXY (04/11/13): replace this with a call in trc_ini_medusa137 !! AXY (25/02/10)138 !! call routine for populating CCD array if this is the first time-step139 !! IF( kt == nittrc000 ) CALL medusa_ccd( kt )140 141 !! AXY (04/11/13): replace this with a call in trc_ini_medusa142 !! AXY (26/01/12)143 !! call routine for populating river arrays if this is the first time-step144 !! IF( kt == nittrc000 ) CALL medusa_river( kt )145 146 !! AXY (10/02/09)147 !!IF( (jnt == 1) .and. (bdustfer) ) CALL trc_sed_medusa_sbc( kt )148 149 !! JPALM -- 31-03-2016 -- rewrite trc_sed_medusa_sbc.150 !! IF (kt == nittrc000 ) CALL trc_sed_medusa_sbc151 152 !! JPALM -- 20-07-2016 -- adapt dust forcing fields reading and conversion153 !! To read dust dep in kg-dust/m2/s instead of g-Fe/m2/month154 !! So all forcings and coupling dust dep are in the same SI units155 !! and then convert in mmol-Fe/m2/day156 157 IF( bdustfer ) THEN158 !! IF( kt == nittrc000 .OR. ( kt /= nittrc000 .AND. ntimes_dust > 1 ) ) THEN159 CALL fld_read( kt, 1, sf_dust )160 dust(:,:) = sf_dust(1)%fnow(:,:,1)161 !! ENDIF162 ELSEIF (lk_oasis) THEN163 dust = Dust_in_cpl164 ELSE165 dust(:,:) = 0.0166 ENDIF167 !!168 zirondep(:,:) = 0.e0 !! Initialisation of deposition variables169 zirondep(:,:) = dust(:,:) * Fe_dust_mratio / xfe_mass * 1.e6 * 86400. !! mmol-Fe/m2/d170 171 !! JPALM -- 20-07-2016 -- Zirondep and zsidep are not used.172 !! So comment out the following lines. but keep them173 !! as we may want to used them later on174 !!================================================175 !!176 !! zirondep(:,:,:) = 0.e0 !! Initialisation of deposition variables177 !! zsidep (:,:) = 0.e0178 !!179 !! Iron and Si deposition at the surface180 !! -------------------------------------181 !!182 !! DO jj = 1, jpj183 !! DO ji = 1, jpi184 !! zirondep(ji,jj,1) = (dustsolub * dust(ji,jj) / (55.85 * rmtss) + 3.e-10 / ryyss) &185 !! & * rfact2 / fse3t(ji,jj,1)186 !! zsidep (ji,jj) = 8.8 * 0.075 * dust(ji,jj) * rfact2 / &187 !! & (fse3t(ji,jj,1) * 28.1 * rmtss)188 !! END DO189 !! END DO190 135 191 136 ! sedimentation of detrital nitrogen : upstream scheme … … 412 357 END SUBROUTINE trc_sed_medusa_sbc 413 358 359 !! ====================================================================== 360 !! ====================================================================== 361 !! ====================================================================== 362 363 !! AXY & JPALM (28/02/17) 364 365 SUBROUTINE trc_sed_medusa_dust( kt ) 366 !!--------------------------------------------------------------------- 367 !! *** ROUTINE trc_sed_medusa_dust *** 368 !! 369 !! ** Purpose : compute current dust *before* trc_bio_medusa call 370 !! 371 !! ** Method : does what it says on the tin 372 !!--------------------------------------------------------------------- 373 INTEGER, INTENT( in ) :: kt ! ocean time-step index 374 375 !! AXY (20/11/14): alter this to report on first MEDUSA call 376 IF( kt == nittrc000 ) THEN 377 IF(lwp) WRITE(numout,*) 378 IF(lwp) WRITE(numout,*) ' trc_sed_medusa_dust: MEDUSA dust timestep' 379 IF(lwp) WRITE(numout,*) ' ~~~~~~~' 380 IF(lwp) WRITE(numout,*) ' kt =',kt 381 ENDIF 382 383 !! AXY (04/11/13): replace this with a call in trc_ini_medusa 384 !! AXY (25/02/10) 385 !! call routine for populating CCD array if this is the first time-step 386 !! IF( kt == nittrc000 ) CALL medusa_ccd( kt ) 387 388 !! AXY (04/11/13): replace this with a call in trc_ini_medusa 389 !! AXY (26/01/12) 390 !! call routine for populating river arrays if this is the first time-step 391 !! IF( kt == nittrc000 ) CALL medusa_river( kt ) 392 393 !! AXY (10/02/09) 394 !! IF( (jnt == 1) .and. (bdustfer) ) CALL trc_sed_medusa_sbc( kt ) 395 396 !! JPALM -- 31-03-2016 -- rewrite trc_sed_medusa_sbc. 397 !! IF (kt == nittrc000 ) CALL trc_sed_medusa_sbc 398 399 !! JPALM -- 20-07-2016 -- adapt dust forcing fields reading and conversion 400 !! To read dust dep in kg-dust/m2/s instead of g-Fe/m2/month 401 !! So all forcings and coupling dust dep are in the same SI units 402 !! and then convert in mmol-Fe/m2/day 403 404 IF( bdustfer ) THEN 405 CALL fld_read( kt, 1, sf_dust ) 406 dust(:,:) = sf_dust(1)%fnow(:,:,1) 407 ELSEIF (lk_oasis) THEN 408 dust = Dust_in_cpl 409 ELSE 410 dust(:,:) = 0.0 411 ENDIF 412 !! 413 zirondep(:,:) = 0.e0 !! Initialisation of deposition variables 414 zirondep(:,:) = dust(:,:) * Fe_dust_mratio / xfe_mass * 1.e6 * 86400. !! mmol-Fe/m2/d 415 416 !! JPALM -- 20-07-2016 -- Zirondep and zsidep are not used. 417 !! So comment out the following lines. but keep them 418 !! as we may want to used them later on 419 !!================================================ 420 !! 421 !! zirondep(:,:,:) = 0.e0 !! Initialisation of deposition variables 422 !! zsidep (:,:) = 0.e0 423 !! 424 !! Iron and Si deposition at the surface 425 !! ------------------------------------- 426 !! 427 !! DO jj = 1, jpj 428 !! DO ji = 1, jpi 429 !! zirondep(ji,jj,1) = (dustsolub * dust(ji,jj) / (55.85 * rmtss) + 3.e-10 / ryyss) & 430 !! & * rfact2 / fse3t(ji,jj,1) 431 !! zsidep (ji,jj) = 8.8 * 0.075 * dust(ji,jj) * rfact2 / & 432 !! & (fse3t(ji,jj,1) * 28.1 * rmtss) 433 !! END DO 434 !! END DO 435 436 END SUBROUTINE trc_sed_medusa_dust 437 414 438 #else 415 439 !!====================================================================== … … 425 449 !!====================================================================== 426 450 END MODULE trcsed_medusa 451 -
branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcsms_medusa.F90
r5841 r7766 58 58 ENDIF 59 59 60 CALL trc_avg_medusa( kt ) 60 CALL trc_avg_medusa( kt ) ! rolling average module 61 61 # if defined key_debug_medusa 62 IF(lwp) WRITE(numout,*) ' MEDUSA done trc_avg_medusa' 63 CALL flush(numout) 62 IF(lwp) WRITE(numout,*) ' MEDUSA done trc_avg_medusa' 63 CALL flush(numout) 64 # endif 65 66 CALL trc_opt_medusa( kt ) ! optical model 67 # if defined key_debug_medusa 68 IF(lwp) WRITE(numout,*) ' MEDUSA done trc_opt_medusa' 69 CALL flush(numout) 64 70 # endif 65 71 66 CALL trc_opt_medusa( kt ) ! optical model 72 !! AXY & JPALM (28/02/17): call dust before trc_bio_medusa (because of coupling) 73 CALL trc_sed_medusa_dust( kt ) ! dust submodel 67 74 # if defined key_debug_medusa 68 IF(lwp) WRITE(numout,*) ' MEDUSA done trc_opt_medusa'69 75 IF(lwp) WRITE(numout,*) ' MEDUSA done trc_sed_medusa_dust' 76 CALL flush(numout) 70 77 # endif 71 78 … … 75 82 CALL flush(numout) 76 83 # else 77 CALL trc_bio_medusa( kt ) 78 # if defined key_debug_medusa79 80 81 # endif82 83 CALL trc_sed_medusa( kt ) 84 # if defined key_debug_medusa85 86 87 # endif84 CALL trc_bio_medusa( kt ) ! biological model 85 # if defined key_debug_medusa 86 IF(lwp) WRITE(numout,*) ' MEDUSA done trc_bio_medusa' 87 CALL flush(numout) 88 # endif 89 90 CALL trc_sed_medusa( kt ) ! sedimentation model 91 # if defined key_debug_medusa 92 IF(lwp) WRITE(numout,*) ' MEDUSA done trc_sed_medusa' 93 CALL flush(numout) 94 # endif 88 95 # endif 89 96 … … 103 110 !!====================================================================== 104 111 END MODULE trcsms_medusa 112 -
branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/trcrst.F90
r7498 r7766 686 686 !!---------------------------------------------------------------------- 687 687 INTEGER :: jk, jn 688 REAL(wp) :: ztraf, zmin, zmax, zmean, zdrift 688 REAL(wp) :: ztraf, zmin, zmax, zmean, zdrift, areasf 689 689 REAL(wp), DIMENSION(jpi,jpj) :: zvol 690 690 !!---------------------------------------------------------------------- … … 696 696 ENDIF 697 697 ! 698 zvol(:,:) = e1e2t(:,:) * fse3t_a(:,:,1) * tmask(:,:,1) 698 zvol(:,:) = e1e2t(:,:) * fse3t_a(:,:,1) * tmask(:,:,1) 699 areasf = glob_sum(zvol(:,:)) 699 700 DO jn = 1, jptra 700 701 ztraf = glob_sum( tra(:,:,1,jn) * zvol(:,:) ) … … 705 706 CALL mpp_max( zmax ) ! max over the global domain 706 707 END IF 707 zmean = ztraf / area tot708 zmean = ztraf / areasf 708 709 IF(lwp) WRITE(numout,9001) jn, TRIM( ctrcnm(jn) ), zmean, zmin, zmax 709 710 END DO … … 734 735 zvol(:,:) = e1e2t(:,:) * fse3t_a(:,:,1) * tmask(:,:,1) 735 736 ztraf = glob_sum( dgtr(:,:) * zvol(:,:) ) 736 areasf = glob_sum(e1e2t(:,:) * tmask(:,:,1) ) 737 !! areasf = glob_sum(e1e2t(:,:) * tmask(:,:,1) ) 738 areasf = glob_sum(zvol(:,:)) 737 739 zmin = MINVAL( dgtr(:,:), mask= ((tmask(:,:,1).NE.0.)) ) 738 740 zmax = MAXVAL( dgtr(:,:), mask= ((tmask(:,:,1).NE.0.)) ) … … 741 743 CALL mpp_max( zmax ) ! max over the global domain 742 744 END IF 743 zmean = ztraf / area tot745 zmean = ztraf / areasf 744 746 IF(lwp) WRITE(numout,9002) TRIM( names ), zmean, zmin, zmax 745 747 ! -
branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/trcstp.F90
r7709 r7766 146 146 ztrai = ztrai + glob_sum( trn(:,:,:,jn) * cvol(:,:,:) ) 147 147 END DO 148 IF( lwp) WRITE(numstr,9300) kt, ztrai / areatot148 IF( numstr /= -1 ) WRITE(numstr,9300) kt, ztrai / areatot 149 149 9300 FORMAT(i10,e18.10) 150 150 ENDIF
Note: See TracChangeset
for help on using the changeset viewer.
