Changeset 13193
- Timestamp:
- 2020-07-01T15:42:06+02:00 (4 years ago)
- Location:
- NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3
- Files:
-
- 1 deleted
- 117 edited
- 2 copied
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3
- Property svn:externals
-
old new 8 8 9 9 # SETTE 10 ^/utils/CI/sette@ HEADsette10 ^/utils/CI/sette@12931 sette
-
- Property svn:externals
-
NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/cfgs/AGRIF_DEMO/EXPREF/1_namelist_cfg
r12724 r13193 81 81 ln_blk = .true. ! Bulk formulation (T => fill namsbc_blk ) 82 82 ! Sea-ice : 83 nn_ice = 2 ! =2 or 3 automatically for SI3 or CICE ("key_si3" or "key_cice") 84 ! except in AGRIF zoom where it has to be specified 83 nn_ice = 2 ! =0 no ice boundary condition 84 ! ! =1 use observed ice-cover ( => fill namsbc_iif ) 85 ! ! =2 or 3 for SI3 and CICE, respectively 85 86 ! Misc. options of sbc : 86 87 ln_traqsr = .true. ! Light penetration in the ocean (T => fill namtra_qsr) -
NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/cfgs/AGRIF_DEMO/EXPREF/2_namelist_cfg
r12724 r13193 78 78 ln_blk = .true. ! Bulk formulation (T => fill namsbc_blk ) 79 79 ! Sea-ice : 80 nn_ice = 2 ! =0 Use SI3 model 80 nn_ice = 2 ! =0 no ice boundary condition 81 ! ! =1 use observed ice-cover ( => fill namsbc_iif ) 82 ! ! =2 or 3 for SI3 and CICE, respectively 81 83 ! Misc. options of sbc : 82 84 ln_traqsr = .true. ! Light penetration in the ocean (T => fill namtra_qsr) -
NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/cfgs/AGRIF_DEMO/EXPREF/3_namelist_cfg
r12724 r13193 78 78 ln_blk = .true. ! Bulk formulation (T => fill namsbc_blk ) 79 79 ! Sea-ice : 80 nn_ice = 2 ! =0 Use SI3 model 80 nn_ice = 2 ! =0 no ice boundary condition 81 ! ! =1 use observed ice-cover ( => fill namsbc_iif ) 82 ! ! =2 or 3 for SI3 and CICE, respectively 81 83 ! Misc. options of sbc : 82 84 ln_traqsr = .true. ! Light penetration in the ocean (T => fill namtra_qsr) -
NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/cfgs/AGRIF_DEMO/EXPREF/namelist_cfg
r12724 r13193 81 81 ln_blk = .true. ! Bulk formulation (T => fill namsbc_blk ) 82 82 ! Sea-ice : 83 nn_ice = 2 ! =2 or 3 automatically for SI3 or CICE ("key_si3" or "key_cice") 84 ! except in AGRIF zoom where it has to be specified 83 nn_ice = 2 ! =0 no ice boundary condition 84 ! ! =1 use observed ice-cover ( => fill namsbc_iif ) 85 ! ! =2 or 3 for SI3 and CICE, respectively 85 86 ! Misc. options of sbc : 86 87 ln_traqsr = .true. ! Light penetration in the ocean (T => fill namtra_qsr) -
NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/cfgs/C1D_PAPA/EXPREF/file_def_nemo-oce.xml
r9799 r13193 53 53 <file id="file4" name_suffix="_grid_W" description="ocean W grid variables" > 54 54 <field field_ref="e3w" /> 55 <field field_ref="woce" name="wo" />56 55 <field field_ref="avt" name="difvho" /> 57 56 </file> -
NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/cfgs/C1D_PAPA/EXPREF/namelist_cfg
r12724 r13193 49 49 &namdom ! time and space domain 50 50 !----------------------------------------------------------------------- 51 ln_linssh = .true. ! =T linear free surface ==>> model level are fixed in time 52 ! 51 53 rn_Dt = 360. ! time step for the dynamics and tracer 52 54 / … … 358 360 &namdyn_spg ! surface pressure gradient (default: NO selection) 359 361 !----------------------------------------------------------------------- 360 ln_dynspg_ts = .true. ! split-explicit free surface361 ln_bt_fw = .false. ! Forward integration of barotropic Eqs.362 ln_bt_av = .true. ! Time filtering of barotropic variables363 362 / 364 363 !----------------------------------------------------------------------- -
NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/cfgs/C1D_PAPA/MY_SRC/usrdef_zgr.F90
r12377 r13193 30 30 PUBLIC usr_def_zgr ! called by domzgr.F90 31 31 32 !! * Substitutions 33 # include "do_loop_substitute.h90" 32 34 !!---------------------------------------------------------------------- 33 35 !! NEMO/OCE 4.0 , NEMO Consortium (2018) … … 157 159 pe3vw(:,:,jk) = pe3w_1d (jk) 158 160 END DO 159 DO jj = 1, jpj ! bottom scale factors and depth at T- and W-points 160 DO ji = 1, jpi 161 ik = k_bot(ji,jj) 162 pdepw(ji,jj,ik+1) = MIN( zht(ji,jj) , pdepw_1d(ik+1) ) 163 pe3t (ji,jj,ik ) = pdepw(ji,jj,ik+1) - pdepw(ji,jj,ik) 164 pe3t (ji,jj,ik+1) = pe3t (ji,jj,ik ) 165 ! 166 pdept(ji,jj,ik ) = pdepw(ji,jj,ik ) + pe3t (ji,jj,ik ) * 0.5_wp 167 pdept(ji,jj,ik+1) = pdepw(ji,jj,ik+1) + pe3t (ji,jj,ik+1) * 0.5_wp 168 pe3w (ji,jj,ik+1) = pdept(ji,jj,ik+1) - pdept(ji,jj,ik) ! = pe3t (ji,jj,ik ) 169 END DO 170 END DO 161 ! bottom scale factors and depth at T- and W-points 162 DO_2D_11_11 163 ik = k_bot(ji,jj) 164 pdepw(ji,jj,ik+1) = MIN( zht(ji,jj) , pdepw_1d(ik+1) ) 165 pe3t (ji,jj,ik ) = pdepw(ji,jj,ik+1) - pdepw(ji,jj,ik) 166 pe3t (ji,jj,ik+1) = pe3t (ji,jj,ik ) 167 ! 168 pdept(ji,jj,ik ) = pdepw(ji,jj,ik ) + pe3t (ji,jj,ik ) * 0.5_wp 169 pdept(ji,jj,ik+1) = pdepw(ji,jj,ik+1) + pe3t (ji,jj,ik+1) * 0.5_wp 170 pe3w (ji,jj,ik+1) = pdept(ji,jj,ik+1) - pdept(ji,jj,ik) ! = pe3t (ji,jj,ik ) 171 END_2D 171 172 ! ! bottom scale factors and depth at U-, V-, UW and VW-points 172 173 ! ! usually Computed as the minimum of neighbooring scale factors -
NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/cfgs/ORCA2_ICE_ABL/EXPREF/namelist_cfg
r12724 r13193 84 84 ln_abl = .true. ! ABL formulation (T => fill namsbc_abl ) 85 85 ! Sea-ice : 86 nn_ice = 2 ! =2 or 3 automatically for SI3 or CICE ("key_si3" or "key_cice") 87 ! except in AGRIF zoom where it has to be specified 86 nn_ice = 2 ! =0 no ice boundary condition 87 ! ! =1 use observed ice-cover ( => fill namsbc_iif ) 88 ! ! =2 or 3 for SI3 and CICE, respectively 88 89 ! Misc. options of sbc : 89 90 ln_traqsr = .true. ! Light penetration in the ocean (T => fill namtra_qsr) -
NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/cfgs/ORCA2_ICE_PISCES/EXPREF/namelist_cfg
r12724 r13193 80 80 ln_blk = .true. ! Bulk formulation (T => fill namsbc_blk ) 81 81 ! Sea-ice : 82 nn_ice = 2 ! =2 or 3 automatically for SI3 or CICE ("key_si3" or "key_cice") 83 ! except in AGRIF zoom where it has to be specified 82 nn_ice = 2 ! =0 no ice boundary condition 83 ! ! =1 use observed ice-cover ( => fill namsbc_iif ) 84 ! ! =2 or 3 for SI3 and CICE, respectively 84 85 ! Misc. options of sbc : 85 86 ln_traqsr = .true. ! Light penetration in the ocean (T => fill namtra_qsr) -
NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/cfgs/ORCA2_ICE_PISCES/EXPREF/namelist_top_cfg
r12377 r13193 20 20 ! 21 21 ln_trcdta = .true. ! Initialisation from data input file (T) or not (F) 22 ln_trcbc = . true.! Enables Boundary conditions22 ln_trcbc = .false. ! Enables Boundary conditions 23 23 ! ! ! ! ! ! 24 24 ! ! name ! title of the field ! units ! init ! sbc ! cbc ! obc ! -
NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/cfgs/ORCA2_OFF_PISCES/EXPREF/namelist_top_cfg
r12377 r13193 20 20 ! 21 21 ln_trcdta = .true. ! Initialisation from data input file (T) or not (F) 22 ln_trcbc = . true.! Enables Boundary conditions22 ln_trcbc = .false. ! Enables Boundary conditions 23 23 ! ! ! ! ! ! 24 24 ! ! name ! title of the field ! units ! init ! sbc ! cbc ! obc ! -
NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/cfgs/ORCA2_SAS_ICE/EXPREF/namelist_cfg
r12377 r13193 59 59 nn_fsbc = 1 ! frequency of SBC module call 60 60 ln_blk = .true. ! Bulk formulation (T => fill namsbc_blk ) 61 nn_ice = 2 ! =2 sea-ice model ("key_SI3" or "key_cice") 61 nn_ice = 2 ! =0 no ice boundary condition 62 ! ! =1 use observed ice-cover ( => fill namsbc_iif ) 63 ! ! =2 or 3 for SI3 and CICE, respectively 62 64 / 63 65 !----------------------------------------------------------------------- -
NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/cfgs/SHARED/field_def_nemo-oce.xml
r12377 r13193 1 <?xml version="1.0"?> 1 <?xml version="1.0"?> 2 2 <!-- $id$ --> 3 3 … … 16 16 Configuration of multiple-linear-regression analysis (diamlr) 17 17 ===================================================================================================== 18 18 19 19 This field group configures diamlr for tidal harmonic analysis of field 20 20 ssh: in addition to a regressor for fitting the mean value (diamlr_r101), … … 73 73 74 74 </field_group> 75 76 <!-- 75 76 <!-- 77 77 ============================================================================================================ 78 78 = definition of all existing variables = … … 101 101 </field_group> 102 102 103 <!-- 103 <!-- 104 104 ============================================================================================================ 105 105 Physical ocean model variables … … 108 108 109 109 <!-- T grid --> 110 110 111 111 <field_group id="grid_T" grid_ref="grid_T_2D" > 112 112 <field id="e3t" long_name="T-cell thickness" standard_name="cell_thickness" unit="m" grid_ref="grid_T_3D" /> 113 113 <field id="e3ts" long_name="T-cell thickness" field_ref="e3t" standard_name="cell_thickness" unit="m" grid_ref="grid_T_SFC"/> 114 114 <field id="e3t_0" long_name="Initial T-cell thickness" standard_name="ref_cell_thickness" unit="m" grid_ref="grid_T_3D" /> 115 <field id="e3tb" long_name="bottom T-cell thickness" standard_name="bottom_cell_thickness" unit="m" grid_ref="grid_T_2D"/> 115 <field id="e3tb" long_name="bottom T-cell thickness" standard_name="bottom_cell_thickness" unit="m" grid_ref="grid_T_2D"/> 116 116 <field id="e3t_300" field_ref="e3t" grid_ref="grid_T_zoom_300" detect_missing_value="true" /> 117 117 <field id="e3t_vsum300" field_ref="e3t_300" grid_ref="grid_T_vsum" detect_missing_value="true" /> 118 118 <field id="masscello" long_name="Sea Water Mass per unit area" standard_name="sea_water_mass_per_unit_area" unit="kg/m2" grid_ref="grid_T_3D"/> 119 <field id="volcello" long_name="Ocean Volume" standard_name="ocean_volume" unit="m3" grid_ref="grid_T_3D"/> 119 <field id="volcello" long_name="Ocean Volume" standard_name="ocean_volume" unit="m3" grid_ref="grid_T_3D"/> 120 120 <field id="toce" long_name="temperature" standard_name="sea_water_potential_temperature" unit="degC" grid_ref="grid_T_3D"/> 121 121 <field id="toce_e3t" long_name="temperature (thickness weighted)" unit="degC" grid_ref="grid_T_3D" > toce * e3t </field > … … 146 146 <field id="sst_cs" long_name="Delta SST of cool skin" unit="degC" /> 147 147 <field id="temp_3m" long_name="temperature at 3m" unit="degC" /> 148 148 149 149 <field id="sss" long_name="sea surface salinity" standard_name="sea_surface_salinity" unit="1e-3" /> 150 150 <field id="sss2" long_name="square of sea surface salinity" unit="1e-6" > sss * sss </field > … … 152 152 <field id="sssmin" long_name="min of sea surface salinity" field_ref="sss" operation="minimum" /> 153 153 <field id="sbs" long_name="sea bottom salinity" unit="0.001" /> 154 <field id="somint" long_name="vertical integral of salinity times density" standard_name="integral_wrt_depth_of_product_of_density_and_salinity" unit="(kg m2) x (1e-3)" /> 155 156 <field id="taubot" long_name="bottom stress module" unit="N/m2" /> 154 <field id="somint" long_name="vertical integral of salinity times density" standard_name="integral_wrt_depth_of_product_of_density_and_salinity" unit="(kg m2) x (1e-3)" /> 155 156 <field id="taubot" long_name="bottom stress module" unit="N/m2" /> 157 157 158 158 <!-- Case EOS = TEOS-10 : output potential temperature --> … … 295 295 <field id="us_y" long_name="j component of Stokes drift" unit="m/s" /> 296 296 </field_group> 297 297 298 298 <!-- SBC --> 299 299 <field_group id="SBC" > <!-- time step automaticaly defined based on nn_fsbc --> … … 311 311 <field id="precip" long_name="Total precipitation" standard_name="precipitation_flux" unit="kg/m2/s" /> 312 312 <field id="wclosea" long_name="closed sea empmr correction" standard_name="closea_empmr" unit="kg/m2/s" /> 313 313 314 314 <field id="qt" long_name="Net Downward Heat Flux" standard_name="surface_downward_heat_flux_in_sea_water" unit="W/m2" /> 315 315 <field id="qns" long_name="non solar Downward Heat Flux" unit="W/m2" /> … … 321 321 <field id="taum" long_name="wind stress module" standard_name="magnitude_of_surface_downward_stress" unit="N/m2" /> 322 322 <field id="wspd" long_name="wind speed module" standard_name="wind_speed" unit="m/s" /> 323 323 324 324 <!-- * variable relative to atmospheric pressure forcing : available with ln_apr_dyn --> 325 325 <field id="ssh_ib" long_name="Inverse barometer sea surface height" standard_name="sea_surface_height_correction_due_to_air_pressure_at_low_frequency" unit="m" /> … … 369 369 <field id="taum_oce" long_name="wind stress module over open ocean" standard_name="magnitude_of_surface_downward_stress" unit="N/m2" /> 370 370 371 <!-- variables computed by the bulk parameterization algorithms (ln_blk) --> 372 <field id="Cd_oce" long_name="Drag coefficient over open ocean" standard_name="drag_coefficient_water" unit="" /> 373 <field id="Ce_oce" long_name="Evaporaion coefficient over open ocean" standard_name="evap_coefficient_water" unit="" /> 374 <field id="Ch_oce" long_name="Sensible heat coefficient over open ocean" standard_name="sensible_heat_coefficient_water" unit="" /> 375 <field id="theta_zt" long_name="Potential air temperature at z=zt" standard_name="potential_air_temperature_at_zt" unit="degC" /> 376 <field id="q_zt" long_name="Specific air humidity at z=zt" standard_name="specific_air_humidity_at_zt" unit="kg/kg" /> 377 <field id="theta_zu" long_name="Potential air temperature at z=zu" standard_name="potential_air_temperature_at_zu" unit="degC" /> 378 <field id="q_zu" long_name="Specific air humidity at z=zu" standard_name="specific_air_humidity_at_zu" unit="kg/kg" /> 379 <field id="ssq" long_name="Saturation specific humidity of air at z=0" standard_name="surface_air_saturation_spec_humidity" unit="kg/kg" /> 380 <field id="wspd_blk" long_name="Bulk wind speed at z=zu" standard_name="bulk_wind_speed_at_zu" unit="m/s" /> 381 <!-- ln_blk + key_si3 --> 382 <field id="Cd_ice" long_name="Drag coefficient over ice" standard_name="drag_coefficient_ice" unit="" /> 383 <field id="Ce_ice" long_name="Evaporaion coefficient over ice" standard_name="evap_coefficient_ice" unit="" /> 384 <field id="Ch_ice" long_name="Sensible heat coefficient over ice" standard_name="sensible_heat_coefficient_ice" unit="" /> 385 371 386 <!-- available key_oasis3 --> 372 387 <field id="snow_ao_cea" long_name="Snow over ice-free ocean (cell average)" standard_name="snowfall_flux" unit="kg/m2/s" /> … … 405 420 <!-- ice field (nn_ice=1) --> 406 421 <field id="ice_cover" long_name="Ice fraction" standard_name="sea_ice_area_fraction" unit="1" /> 407 422 408 423 <!-- dilution --> 409 424 <field id="emp_x_sst" long_name="Concentration/Dilution term on SST" unit="kg*degC/m2/s" /> 410 <field id="emp_x_sss" long_name="Concentration/Dilution term on SSS" unit="kg*1e-3/m2/s" /> 425 <field id="emp_x_sss" long_name="Concentration/Dilution term on SSS" unit="kg*1e-3/m2/s" /> 411 426 <field id="rnf_x_sst" long_name="Runoff term on SST" unit="kg*degC/m2/s" /> 412 427 <field id="rnf_x_sss" long_name="Runoff term on SSS" unit="kg*1e-3/m2/s" /> 413 428 414 429 <!-- sbcssm variables --> 415 430 <field id="sst_m" unit="degC" /> … … 422 437 423 438 </field_group> 424 439 425 440 426 441 </field_group> <!-- SBC --> 427 442 428 443 <!-- ABL --> 429 444 <field_group id="ABL" > <!-- time step automaticaly defined based on nn_fsbc --> … … 456 471 <field id="uz1_dta" long_name="DTA i-horizontal velocity" standard_name="dta_x_velocity" unit="m/s" /> 457 472 <field id="vz1_dta" long_name="DTA j-horizontal velocity" standard_name="dta_y_velocity" unit="m/s" /> 458 <field id="uvz1_dta" long_name="DTA wind speed module" standard_name="dta_wind_speed" unit="m/s" > sqrt( uz1_dta^2 + vz1_dta^2 ) </field> 473 <field id="uvz1_dta" long_name="DTA wind speed module" standard_name="dta_wind_speed" unit="m/s" > sqrt( uz1_dta^2 + vz1_dta^2 ) </field> 459 474 <field id="tz1_dta" long_name="DTA potential temperature" standard_name="dta_theta" unit="K" /> 460 475 <field id="qz1_dta" long_name="DTA specific humidity" standard_name="dta_qspe" unit="kg/kg" /> … … 462 477 <field id="uz1_geo" long_name="GEO i-horizontal velocity" standard_name="geo_x_velocity" unit="m/s" /> 463 478 <field id="vz1_geo" long_name="GEO j-horizontal velocity" standard_name="geo_y_velocity" unit="m/s" /> 464 <field id="uvz1_geo" long_name="GEO wind speed module" standard_name="geo_wind_speed" unit="m/s" > sqrt( uz1_geo^2 + vz1_geo^2 ) </field> 479 <field id="uvz1_geo" long_name="GEO wind speed module" standard_name="geo_wind_speed" unit="m/s" > sqrt( uz1_geo^2 + vz1_geo^2 ) </field> 465 480 </field_group> 466 481 467 482 </field_group> <!-- ABL --> 468 483 469 484 470 485 <!-- U grid --> 471 486 472 487 <field_group id="grid_U" grid_ref="grid_U_2D"> 473 488 <field id="e2u" long_name="U-cell width in meridional direction" standard_name="cell_width" unit="m" /> … … 478 493 <field id="uoce_e3u" long_name="ocean current along i-axis (thickness weighted)" unit="m/s" grid_ref="grid_U_3D" > uoce * e3u </field> 479 494 <field id="uoce_e3u_vsum" long_name="ocean current along i-axis * e3u summed on the vertical" field_ref="uoce_e3u" unit="m3/s" grid_ref="grid_U_vsum"/> 480 <field id="uocetr_vsum" long_name="ocean transport along i-axis summed on the vertical" field_ref="e2u" unit="m3/s"> this * uoce_e3u_vsum </field> 495 <field id="uocetr_vsum" long_name="ocean transport along i-axis summed on the vertical" field_ref="e2u" unit="m3/s"> this * uoce_e3u_vsum </field> 481 496 482 497 <field id="uocetr_vsum_op" long_name="ocean current along i-axis * e3u * e2u summed on the vertical" read_access="true" freq_op="1mo" field_ref="e2u" unit="m3/s"> @uocetr_vsum </field> 483 <field id="uocetr_vsum_cumul" long_name="ocean current along i-axis * e3u * e2u cumulated from southwest point" freq_offset="_reset_" operation="instant" freq_op="1mo" unit="m3/s" /> 498 <field id="uocetr_vsum_cumul" long_name="ocean current along i-axis * e3u * e2u cumulated from southwest point" freq_offset="_reset_" operation="instant" freq_op="1mo" unit="m3/s" /> 484 499 <field id="msftbarot" long_name="ocean_barotropic_mass_streamfunction" unit="kg s-1" > uocetr_vsum_cumul * $rau0 </field> 485 500 … … 534 549 <field id="udiff_salttr" long_name="ocean diffusion salt transport along i-axis" standard_name="ocean_salt_x_transport_due_to_diffusion" unit="1e-3*kg/s" /> 535 550 </field_group> 536 551 537 552 <!-- V grid --> 538 553 539 554 <field_group id="grid_V" grid_ref="grid_V_2D"> 540 555 <field id="e1v" long_name="V-cell width in longitudinal direction" standard_name="cell_width" unit="m" /> … … 593 608 <field id="vdiff_salttr" long_name="ocean diffusion salt transport along j-axis" standard_name="ocean_salt_y_transport_due_to_diffusion" unit="1e-3*kg/s" /> 594 609 </field_group> 595 610 596 611 <!-- W grid --> 597 612 598 613 <field_group id="grid_W" grid_ref="grid_W_3D"> 599 614 <field id="e3w" long_name="W-cell thickness" standard_name="cell_thickness" unit="m" /> 600 615 <field id="woce" long_name="ocean vertical velocity" standard_name="upward_sea_water_velocity" unit="m/s" /> 601 <field id="woce_e3w" long_name="ocean vertical velocity * e3w" unit="m2/s" > woce * e3w </field> 616 <field id="woce_e3w" long_name="ocean vertical velocity * e3w" unit="m2/s" > woce * e3w </field> 602 617 <field id="wocetr_eff" long_name="effective ocean vertical transport" unit="m3/s" /> 603 618 … … 609 624 610 625 <field id="avt" long_name="vertical eddy diffusivity" standard_name="ocean_vertical_heat_diffusivity" unit="m2/s" /> 611 <field id="avt_e3w" long_name="vertical heat diffusivity * e3w" unit="m3/s" > avt * e3w </field> 626 <field id="avt_e3w" long_name="vertical heat diffusivity * e3w" unit="m3/s" > avt * e3w </field> 612 627 <field id="logavt" long_name="logarithm of vertical eddy diffusivity" standard_name="ocean_vertical_heat_diffusivity" unit="m2/s" /> 613 628 <field id="avm" long_name="vertical eddy viscosity" standard_name="ocean_vertical_momentum_diffusivity" unit="m2/s" /> 614 <field id="avm_e3w" long_name="vertical eddy viscosity * e3w" unit="m3/s" > avm * e3w </field> 629 <field id="avm_e3w" long_name="vertical eddy viscosity * e3w" unit="m3/s" > avm * e3w </field> 615 630 616 631 <!-- avs: /= avt with ln_zdfddm=T --> 617 632 <field id="avs" long_name="salt vertical eddy diffusivity" standard_name="ocean_vertical_salt_diffusivity" unit="m2/s" /> 618 <field id="avs_e3w" long_name="vertical salt diffusivity * e3w" unit="m3/s" > avs * e3w </field> 633 <field id="avs_e3w" long_name="vertical salt diffusivity * e3w" unit="m3/s" > avs * e3w </field> 619 634 <field id="logavs" long_name="logarithm of salt vertical eddy diffusivity" standard_name="ocean_vertical_heat_diffusivity" unit="m2/s" /> 620 635 621 636 <!-- avt_evd and avm_evd: available with ln_zdfevd --> 622 637 <field id="avt_evd" long_name="convective enhancement of vertical diffusivity" standard_name="ocean_vertical_tracer_diffusivity_due_to_convection" unit="m2/s" /> 623 <field id="avt_evd_e3w" long_name="convective enhancement to vertical diffusivity * e3w " unit="m3/s" > avt_evd * e3w </field> 638 <field id="avt_evd_e3w" long_name="convective enhancement to vertical diffusivity * e3w " unit="m3/s" > avt_evd * e3w </field> 624 639 <field id="avm_evd" long_name="convective enhancement of vertical viscosity" standard_name="ocean_vertical_momentum_diffusivity_due_to_convection" unit="m2/s" /> 625 640 … … 634 649 <field id="wstokes" long_name="Stokes Drift vertical velocity" standard_name="upward_StokesDrift_velocity" unit="m/s" /> 635 650 636 <!-- variables available with diaar5 --> 651 <!-- variables available with diaar5 --> 637 652 <field id="w_masstr" long_name="vertical mass transport" standard_name="upward_ocean_mass_transport" unit="kg/s" /> 638 653 <field id="w_masstr2" long_name="square of vertical mass transport" standard_name="square_of_upward_ocean_mass_transport" unit="kg2/s2" /> 639 654 640 655 </field_group> 641 656 642 657 <!-- F grid --> 643 658 <!-- AGRIF sponge --> … … 694 709 </field_group> 695 710 696 711 697 712 <!-- variables available with ln_floats --> 698 713 … … 709 724 <!-- variables available with iceberg trajectories --> 710 725 711 <field_group id="icbvar" domain_ref="grid_T" > 726 <field_group id="icbvar" domain_ref="grid_T" > 712 727 <field id="berg_melt" long_name="icb melt rate of icebergs" unit="kg/m2/s" /> 713 728 <field id="berg_melt_hcflx" long_name="icb heat flux to ocean due to melting heat content" unit="J/m2/s" /> … … 727 742 </field_group> 728 743 729 <!-- Poleward transport : ptr --> 730 <field_group id="diaptr" > 744 <!-- Poleward transport : ptr --> 745 <field_group id="diaptr" > 731 746 <field id="zomsf" long_name="Overturning Stream-Function : All basins" unit="Sv" grid_ref="grid_znl_W_3D" /> 732 747 <field id="zotem" long_name="Zonal Mean Temperature : All basins" unit="degree_C" grid_ref="grid_znl_T_3D" /> … … 736 751 <field id="sopstove" long_name="Overturning Salt Transport: All basins" unit="Giga g/s" grid_ref="grid_znl_T_2D" /> 737 752 <field id="sophtbtr" long_name="Barotropic Heat Transport: All basins" unit="PW" grid_ref="grid_znl_T_2D" /> 738 <field id="sopstbtr" long_name="Barotropic Salt Transport: All basins" unit="Giga g/s" grid_ref="grid_znl_T_2D" /> 753 <field id="sopstbtr" long_name="Barotropic Salt Transport: All basins" unit="Giga g/s" grid_ref="grid_znl_T_2D" /> 739 754 <field id="sophtadv" long_name="Advective Heat Transport: All basins" unit="PW" grid_ref="grid_znl_T_2D" /> 740 755 <field id="sopstadv" long_name="Advective Salt Transport: All basins" unit="Giga g/s" grid_ref="grid_znl_T_2D" /> … … 742 757 <field id="sopstldf" long_name="Diffusive Salt Transport: All basins" unit="Giga g/s" grid_ref="grid_znl_T_2D" /> 743 758 <field id="sophtvtr" long_name="Heat Transport : All basins" unit="PW" grid_ref="grid_znl_T_2D" /> 744 <field id="sopstvtr" long_name="Salt Transport : All basins" unit="Giga g/s" grid_ref="grid_znl_T_2D" /> 759 <field id="sopstvtr" long_name="Salt Transport : All basins" unit="Giga g/s" grid_ref="grid_znl_T_2D" /> 745 760 <field id="sophteiv" long_name="Heat Transport from mesoscale eddy advection: All basins" unit="PW" grid_ref="grid_znl_T_2D" /> 746 761 <field id="sopsteiv" long_name="Salt Transport from mesoscale eddy advection : All basins" unit="Giga g/s" grid_ref="grid_znl_T_2D" /> … … 758 773 759 774 760 <!-- 775 <!-- 761 776 ============================================================================================================ 762 777 Physical ocean model trend diagnostics : temperature, KE, PE, momentum … … 899 914 <field id="ketrd_zdf" long_name="ke-trend: vertical diffusion" unit="W/s^3" /> 900 915 <field id="ketrd_tau" long_name="ke-trend: wind stress " unit="W/s^3" grid_ref="grid_T_2D" /> 901 <field id="ketrd_bfr" long_name="ke-trend: bottom friction (explicit)" unit="W/s^3" /> 902 <field id="ketrd_bfri" long_name="ke-trend: bottom friction (implicit)" unit="W/s^3" /> 903 <field id="ketrd_atf" long_name="ke-trend: asselin time filter trend" unit="W/s^3" /> 916 <field id="ketrd_bfr" long_name="ke-trend: bottom friction (explicit)" unit="W/s^3" /> 917 <field id="ketrd_bfri" long_name="ke-trend: bottom friction (implicit)" unit="W/s^3" /> 918 <field id="ketrd_atf" long_name="ke-trend: asselin time filter trend" unit="W/s^3" /> 904 919 <field id="ketrd_convP2K" long_name="ke-trend: conversion (potential to kinetic)" unit="W/s^3" /> 905 <field id="KE" long_name="kinetic energy: u(n)*u(n+1)/2" unit="W/s^2" /> 920 <field id="KE" long_name="kinetic energy: u(n)*u(n+1)/2" unit="W/s^2" /> 906 921 907 922 <!-- variables available when explicit lateral mixing is used (ln_dynldf_OFF=F) --> 908 <field id="dispkexyfo" long_name="KE-trend: lateral mixing induced dissipation" standard_name="ocean_kinetic_energy_dissipation_per_unit_area_due_to_xy_friction" unit="W/m^2" grid_ref="grid_T_2D" /> 909 <field id="dispkevfo" long_name="KE-trend: vertical mixing induced dissipation" standard_name="ocean_kinetic_energy_dissipation_per_unit_area_due_to_vertical_friction" unit="W/m^2" grid_ref="grid_T_2D" /> 923 <field id="dispkexyfo" long_name="KE-trend: lateral mixing induced dissipation" standard_name="ocean_kinetic_energy_dissipation_per_unit_area_due_to_xy_friction" unit="W/m^2" grid_ref="grid_T_2D" /> 924 <field id="dispkevfo" long_name="KE-trend: vertical mixing induced dissipation" standard_name="ocean_kinetic_energy_dissipation_per_unit_area_due_to_vertical_friction" unit="W/m^2" grid_ref="grid_T_2D" /> 910 925 <!-- variables available with ln_traadv_eiv=T and ln_diaeiv=T --> 911 <field id="eketrd_eiv" long_name="EKE-trend due to parameterized eddy advection" standard_name="tendency_of_ocean_eddy_kinetic_energy_content_due_to_parameterized_eddy_advection" unit="W/m^2" grid_ref="grid_T_2D" /> 926 <field id="eketrd_eiv" long_name="EKE-trend due to parameterized eddy advection" standard_name="tendency_of_ocean_eddy_kinetic_energy_content_due_to_parameterized_eddy_advection" unit="W/m^2" grid_ref="grid_T_2D" /> 912 927 913 928 <!-- variables available with ln_PE_trd --> … … 926 941 <field id="petrd_bbc" long_name="pe-trend: geothermal heating" unit="W/m^3" /> 927 942 <field id="petrd_atf" long_name="pe-trend: asselin time filter" unit="W/m^3" /> 928 <field id="PEanom" long_name="potential energy anomaly" unit="1" /> 929 <field id="alphaPE" long_name="partial deriv. of PEanom wrt T" unit="degC-1" /> 930 <field id="betaPE" long_name="partial deriv. of PEanom wrt S" unit="1e3" /> 943 <field id="PEanom" long_name="potential energy anomaly" unit="1" /> 944 <field id="alphaPE" long_name="partial deriv. of PEanom wrt T" unit="degC-1" /> 945 <field id="betaPE" long_name="partial deriv. of PEanom wrt S" unit="1e3" /> 931 946 </field_group> 932 947 … … 945 960 <field id="utrd_zdf" long_name="i-trend: vertical diffusion" unit="m/s^2" /> 946 961 <field id="utrd_tau" long_name="i-trend: wind stress " unit="m/s^2" grid_ref="grid_U_2D" /> 947 <field id="utrd_bfr" long_name="i-trend: bottom friction (explicit)" unit="m/s^2" /> 948 <field id="utrd_bfri" long_name="i-trend: bottom friction (implicit)" unit="m/s^2" /> 949 <field id="utrd_tot" long_name="i-trend: total momentum trend before atf" unit="m/s^2" /> 950 <field id="utrd_atf" long_name="i-trend: asselin time filter trend" unit="m/s^2" /> 962 <field id="utrd_bfr" long_name="i-trend: bottom friction (explicit)" unit="m/s^2" /> 963 <field id="utrd_bfri" long_name="i-trend: bottom friction (implicit)" unit="m/s^2" /> 964 <field id="utrd_tot" long_name="i-trend: total momentum trend before atf" unit="m/s^2" /> 965 <field id="utrd_atf" long_name="i-trend: asselin time filter trend" unit="m/s^2" /> 951 966 </field_group> 952 967 … … 965 980 <field id="vtrd_zdf" long_name="j-trend: vertical diffusion" unit="m/s^2" /> 966 981 <field id="vtrd_tau" long_name="j-trend: wind stress " unit="m/s^2" grid_ref="grid_V_2D" /> 967 <field id="vtrd_bfr" long_name="j-trend: bottom friction (explicit)" unit="m/s^2" /> 968 <field id="vtrd_bfri" long_name="j-trend: bottom friction (implicit)" unit="m/s^2" /> 969 <field id="vtrd_tot" long_name="j-trend: total momentum trend before atf" unit="m/s^2" /> 970 <field id="vtrd_atf" long_name="j-trend: asselin time filter trend" unit="m/s^2" /> 982 <field id="vtrd_bfr" long_name="j-trend: bottom friction (explicit)" unit="m/s^2" /> 983 <field id="vtrd_bfri" long_name="j-trend: bottom friction (implicit)" unit="m/s^2" /> 984 <field id="vtrd_tot" long_name="j-trend: total momentum trend before atf" unit="m/s^2" /> 985 <field id="vtrd_atf" long_name="j-trend: asselin time filter trend" unit="m/s^2" /> 971 986 </field_group> 972 987 973 988 974 <!-- 989 <!-- 975 990 ============================================================================================================ 976 991 Definitions for iodef_demo.xml … … 990 1005 <field field_ref="strd_zdfp_li" name="osaltdiff" /> 991 1006 </field_group> 992 1007 993 1008 <field_group id="mooring" > 994 1009 <field field_ref="toce" name="thetao" long_name="sea_water_potential_temperature" /> … … 999 1014 <field field_ref="avt" name="difvho" long_name="ocean_vertical_heat_diffusivity" /> 1000 1015 <field field_ref="avm" name="difvmo" long_name="ocean_vertical_momentum_diffusivity" /> 1001 1016 1002 1017 <field field_ref="sst" name="tos" long_name="sea_surface_temperature" /> 1003 1018 <field field_ref="sst2" name="tossq" long_name="square_of_sea_surface_temperature" /> … … 1047 1062 <field field_ref="BLT" name="blt" long_name="Barrier Layer Thickness" /> 1048 1063 </field_group> 1049 1064 1050 1065 <field_group id="groupU" > 1051 1066 <field field_ref="uoce" name="uo" long_name="sea_water_x_velocity" /> 1052 1067 <field field_ref="utau" name="tauuo" long_name="surface_downward_x_stress" /> 1053 1068 </field_group> 1054 1069 1055 1070 <field_group id="groupV" > 1056 1071 <field field_ref="voce" name="vo" long_name="sea_water_y_velocity" /> 1057 1072 <field field_ref="vtau" name="tauvo" long_name="surface_downward_y_stress" /> 1058 1073 </field_group> 1059 1074 1060 1075 <field_group id="groupW" > 1061 1076 <field field_ref="woce" name="wo" long_name="ocean vertical velocity" /> … … 1100 1115 </field_group> 1101 1116 1102 <!-- 1117 <!-- 1103 1118 ============================================================================================================ 1104 1119 --> 1105 <!-- output variables for my configuration (example) --> 1106 1120 <!-- output variables for my configuration (example) --> 1121 1107 1122 <field_group id="myvarOCE" > 1108 <!-- grid T --> 1123 <!-- grid T --> 1109 1124 <field field_ref="e3t" name="e3t" long_name="vertical scale factor" /> 1110 1125 <field field_ref="sst" name="tos" long_name="sea_surface_temperature" /> 1111 1126 <field field_ref="sss" name="sos" long_name="sea_surface_salinity" /> 1112 1127 <field field_ref="ssh" name="zos" long_name="sea_surface_height_above_geoid" /> 1113 1114 <!-- grid U --> 1128 1129 <!-- grid U --> 1115 1130 <field field_ref="e3u" name="e3u" long_name="vertical scale factor" /> 1116 1131 <field field_ref="ssu" name="uos" long_name="sea_surface_x_velocity" /> 1117 1118 <!-- grid V --> 1132 1133 <!-- grid V --> 1119 1134 <field field_ref="e3v" name="e3v" long_name="vertical scale factor" /> 1120 <field field_ref="ssv" name="vos" long_name="sea_surface_y_velocity" /> 1121 </field_group> 1135 <field field_ref="ssv" name="vos" long_name="sea_surface_y_velocity" /> 1136 </field_group> 1122 1137 1123 1138 </field_definition> -
NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/cfgs/SHARED/namelist_pisces_ref
r12377 r13193 352 352 ! 353 353 cn_dir = './' ! root directory for the location of the dynamical files 354 ln_ironsed = . true. ! boolean for Fe input from sediments355 ln_ironice = . true. ! boolean for Fe input from sea ice356 ln_hydrofe = . true. ! boolean for from hydrothermal vents354 ln_ironsed = .false. ! boolean for Fe input from sediments 355 ln_ironice = .false. ! boolean for Fe input from sea ice 356 ln_hydrofe = .false. ! boolean for from hydrothermal vents 357 357 sedfeinput = 2.e-9 ! Coastal release of Iron 358 358 distcoast = 5.e3 ! Distance off the coast for Iron from sediments -
NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/cfgs/SHARED/namelist_ref
r12724 r13193 1134 1134 ! ! = 3 as =2 with distinct dissipative an mixing length scale 1135 1135 ln_mxl0 = .true. ! surface mixing length scale = F(wind stress) (T) or not (F) 1136 nn_mxlice = 0 ! type of scaling under sea-ice 1137 ! = 0 no scaling under sea-ice 1138 ! = 1 scaling with constant sea-ice thickness 1139 ! = 2 scaling with mean sea-ice thickness ( only with SI3 sea-ice model ) 1140 ! = 3 scaling with maximum sea-ice thickness 1141 rn_mxlice = 10. ! max constant ice thickness value when scaling under sea-ice ( nn_mxlice=1) 1136 1142 rn_mxl0 = 0.04 ! surface buoyancy lenght scale minimum value 1137 1143 ln_drg = .false. ! top/bottom friction added as boundary condition of TKE … … 1394 1400 &namctl ! Control prints (default: OFF) 1395 1401 !----------------------------------------------------------------------- 1396 sn_cfctl%l_glochk = .FALSE. ! Range sanity checks are local (F) or global (T). Set T for debugging only 1397 sn_cfctl%l_allon = .FALSE. ! IF T activate all options. If F deactivate all unless l_config is T 1398 sn_cfctl%l_config = .TRUE. ! IF .true. then control which reports are written with the following 1399 sn_cfctl%l_runstat = .TRUE. ! switches and which areas produce reports with the proc integer settings. 1400 sn_cfctl%l_trcstat = .FALSE. ! The default settings for the proc integers should ensure 1401 sn_cfctl%l_oceout = .FALSE. ! that all areas report. 1402 sn_cfctl%l_layout = .FALSE. ! 1403 sn_cfctl%l_prtctl = .FALSE. ! 1404 sn_cfctl%l_prttrc = .FALSE. ! 1405 sn_cfctl%l_oasout = .FALSE. ! 1406 sn_cfctl%procmin = 0 ! Minimum area number for reporting [default:0] 1407 sn_cfctl%procmax = 1000000 ! Maximum area number for reporting [default:1000000] 1408 sn_cfctl%procincr = 1 ! Increment for optional subsetting of areas [default:1] 1409 sn_cfctl%ptimincr = 1 ! Timestep increment for writing time step progress info 1410 nn_print = 0 ! level of print (0 no extra print) 1411 nn_ictls = 0 ! start i indice of control sum (use to compare mono versus 1412 nn_ictle = 0 ! end i indice of control sum multi processor runs 1413 nn_jctls = 0 ! start j indice of control over a subdomain) 1414 nn_jctle = 0 ! end j indice of control 1415 nn_isplt = 1 ! number of processors in i-direction 1416 nn_jsplt = 1 ! number of processors in j-direction 1417 ln_timing = .false. ! timing by routine write out in timing.output file 1418 ln_diacfl = .false. ! CFL diagnostics write out in cfl_diagnostics.ascii 1402 sn_cfctl%l_runstat = .TRUE. ! switches and which areas produce reports with the proc integer settings. 1403 sn_cfctl%l_trcstat = .FALSE. ! The default settings for the proc integers should ensure 1404 sn_cfctl%l_oceout = .FALSE. ! that all areas report. 1405 sn_cfctl%l_layout = .FALSE. ! 1406 sn_cfctl%l_prtctl = .FALSE. ! 1407 sn_cfctl%l_prttrc = .FALSE. ! 1408 sn_cfctl%l_oasout = .FALSE. ! 1409 sn_cfctl%procmin = 0 ! Minimum area number for reporting [default:0] 1410 sn_cfctl%procmax = 1000000 ! Maximum area number for reporting [default:1000000] 1411 sn_cfctl%procincr = 1 ! Increment for optional subsetting of areas [default:1] 1412 sn_cfctl%ptimincr = 1 ! Timestep increment for writing time step progress info 1413 nn_print = 0 ! level of print (0 no extra print) 1414 nn_ictls = 0 ! start i indice of control sum (use to compare mono versus 1415 nn_ictle = 0 ! end i indice of control sum multi processor runs 1416 nn_jctls = 0 ! start j indice of control over a subdomain) 1417 nn_jctle = 0 ! end j indice of control 1418 nn_isplt = 1 ! number of processors in i-direction 1419 nn_jsplt = 1 ! number of processors in j-direction 1420 ln_timing = .false. ! timing by routine write out in timing.output file 1421 ln_diacfl = .false. ! CFL diagnostics write out in cfl_diagnostics.ascii 1419 1422 / 1420 1423 !----------------------------------------------------------------------- -
NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/cfgs/WED025/EXPREF/file_def_nemo-ice.xml
r11844 r13193 78 78 </file> 79 79 80 <file id="file22" name_suffix="_SBC_scalar" description="scalar variables" enabled=".true." >81 <!-- global contents -->82 <field field_ref="ibgvol_tot" grid_ref="grid_1point" name="ibgvol_tot" />83 <field field_ref="sbgvol_tot" grid_ref="grid_1point" name="sbgvol_tot" />84 <field field_ref="ibgarea_tot" grid_ref="grid_1point" name="ibgarea_tot" />85 <field field_ref="ibgsalt_tot" grid_ref="grid_1point" name="ibgsalt_tot" />86 <field field_ref="ibgheat_tot" grid_ref="grid_1point" name="ibgheat_tot" />87 <field field_ref="sbgheat_tot" grid_ref="grid_1point" name="sbgheat_tot" />88 89 <!-- global drifts (conservation checks) -->90 <field field_ref="ibgvolume" grid_ref="grid_1point" name="ibgvolume" />91 <field field_ref="ibgsaltco" grid_ref="grid_1point" name="ibgsaltco" />92 <field field_ref="ibgheatco" grid_ref="grid_1point" name="ibgheatco" />93 <field field_ref="ibgheatfx" grid_ref="grid_1point" name="ibgheatfx" />94 95 <!-- global forcings -->96 <field field_ref="ibgfrcvoltop" grid_ref="grid_1point" name="ibgfrcvoltop" />97 <field field_ref="ibgfrcvolbot" grid_ref="grid_1point" name="ibgfrcvolbot" />98 <field field_ref="ibgfrctemtop" grid_ref="grid_1point" name="ibgfrctemtop" />99 <field field_ref="ibgfrctembot" grid_ref="grid_1point" name="ibgfrctembot" />100 <field field_ref="ibgfrcsal" grid_ref="grid_1point" name="ibgfrcsal" />101 <field field_ref="ibgfrchfxtop" grid_ref="grid_1point" name="ibgfrchfxtop" />102 <field field_ref="ibgfrchfxbot" grid_ref="grid_1point" name="ibgfrchfxbot" />103 </file>104 105 80 </file_group> 106 81 -
NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/cfgs/WED025/EXPREF/namelist_cfg
r12724 r13193 5 5 !! namelists 2 - Surface boundary (namsbc, namsbc_flx, namsbc_blk, namsbc_cpl, 6 6 !! namsbc_sas, namtra_qsr, namsbc_rnf, 7 !! nam sbc_isf, namsbc_iscpl, namsbc_apr,7 !! namisf, namsbc_apr, 8 8 !! namsbc_ssr, namsbc_wave, namberg) 9 9 !! 3 - lateral boundary (namlbc, namagrif, nambdy, nambdy_tide) … … 38 38 nn_it000 = 1 ! first time step 39 39 nn_itend = 26280 ! last time step (std 5475) 40 nn_date0 = 19760301 ! date at nit_0000 (format yyyymmdd) used if ln_rstart=F or (ln_rstart=T and nn_rstctl=0 or 1)40 nn_date0 = 20000101 ! date at nit_0000 (format yyyymmdd) used if ln_rstart=F or (ln_rstart=T and nn_rstctl=0 or 1) 41 41 ln_rstart = .false. ! start from rest (F) or from a restart file (T) 42 42 nn_rstctl = 2 ! restart control ==> activated only if ln_rstart=T … … 61 61 ln_tsd_init = .true. ! ocean initialisation 62 62 ln_tsd_dmp = .false. ! T-S restoring (see namtra_dmp) 63 63 64 64 cn_dir = './' ! root directory for the T-S data location 65 !___________!_____________________ ____!___________________!___________!_____________!________!___________!__________________!__________!_______________!66 ! ! file name 67 ! ! 68 sn_tem = ' dta_temp_WED025' , -12 , 'votemper', .true., .true. , 'yearly' , '' , '' , ''69 sn_sal = ' dta_sal_WED025' , -12 , 'vosaline', .true., .true. , 'yearly' , '' , '' , ''65 !___________!_____________________!___________________!___________!_____________!________!___________!__________________!__________!_______________! 66 ! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! 67 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! ! pairing ! filename ! 68 sn_tem = 'WED025_init_JRA_200001.nc', -12 , 'votemper', .false. , .true. , 'yearly' , '' , '' , '' 69 sn_sal = 'WED025_init_JRA_200001.nc', -12 , 'vosaline', .false. , .true. , 'yearly' , '' , '' , '' 70 70 / 71 71 !----------------------------------------------------------------------- … … 116 116 ln_blk = .true. ! Bulk formulation (T => fill namsbc_blk ) 117 117 ! Sea-ice : 118 nn_ice = 2 ! =0 no ice boundary condition 118 nn_ice = 2 ! =0 no ice boundary condition 119 119 ! ! =1 use observed ice-cover ( => fill namsbc_iif ) 120 ! ! =2 or 3 automatically for SI3 or CICE ("key_si3" or "key_cice") 121 ! ! except in AGRIF zoom where it has to be specified 120 ! ! =2 or 3 for SI3 and CICE, respectively 122 121 ln_ice_embd = .false. ! =T embedded sea-ice (pressure + mass and salt exchanges) 123 122 ! ! =F levitating ice (no pressure, mass and salt exchanges) 124 123 ! Misc. options of sbc : 125 124 ln_traqsr = .true. ! Light penetration in the ocean (T => fill namtra_qsr) 126 ln_dm2dc = . true.! daily mean to diurnal cycle on short wave125 ln_dm2dc = .false. ! daily mean to diurnal cycle on short wave 127 126 ln_ssr = .false. ! Sea Surface Restoring on T and/or S (T => fill namsbc_ssr) 128 127 nn_fwb = 0 ! FreshWater Budget: =0 unchecked … … 141 140 ln_NCAR = .true. ! "NCAR" algorithm (Large and Yeager 2008) 142 141 ln_COARE_3p0 = .false. ! "COARE 3.0" algorithm (Fairall et al. 2003) 143 ln_COARE_3p 5 = .false. ! "COARE 3.5" algorithm (Edson et al. 2013)144 ln_ECMWF = .false. ! "ECMWF" algorithm (IFS cycle 31)142 ln_COARE_3p6 = .false. ! "COARE 3.6" algorithm (Edson et al. 2013) 143 ln_ECMWF = .false. ! "ECMWF" algorithm (IFS cycle 45r1) 145 144 146 145 cn_dir = './' ! root directory for the bulk data location … … 148 147 ! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! 149 148 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! ! pairing ! filename ! 150 sn_wndi = 'u10_ core' , 6 , 'U_10_MOD', .true. , .false. , 'yearly' , 'weights_bicubic_core.nc' , 'Uwnd' , ''151 sn_wndj = 'v10_ core' , 6 , 'V_10_MOD', .true. , .false. , 'yearly' , 'weights_bicubic_core.nc' , 'Vwnd' , ''152 sn_qsr = ' qsw_core' , 24 , 'SWDN_MOD', .false. , .false. , 'yearly' , 'weights_bilin_core.nc' , '' , ''153 sn_qlw = ' qlw_core' , 24 , 'LWDN_MOD', .false. , .false. , 'yearly' , 'weights_bilin_core.nc' , '' , ''154 sn_tair = 't10_ core' , 6 , 'T_10_MOD', .true. , .false. , 'yearly' , 'weights_bilin_core.nc' , '' , ''155 sn_humi = 'q10_ core' , 6 , 'Q_10_MOD', .true. , .false. , 'yearly' , 'weights_bilin_core.nc' , '' , ''156 sn_prec = 'precip_ core' , -1 , 'TPRECIP', .true. , .false. , 'yearly' , 'weights_bilin_core.nc' , '' , ''157 sn_snow = 'snow_ core' , -1 , 'SNOW' , .true. , .false. , 'yearly' , 'weights_bilin_core.nc' , '' , ''158 sn_slp = 'slp_ core' , 6 , 'SLP' , .true. , .false. , 'yearly' , 'weights_bilin_core.nc' , '' , ''149 sn_wndi = 'u10_JRA' , 3 , 'uas_10m' , .true. , .false. , 'yearly' , 'weights_bicubic_JRA.nc' , 'Uwnd' , '' 150 sn_wndj = 'v10_JRA' , 3 , 'vas_10m' , .true. , .false. , 'yearly' , 'weights_bicubic_JRA.nc' , 'Vwnd' , '' 151 sn_qsr = 'rsds_JRA' , 3 , 'rsds' , .true. , .false. , 'yearly' , 'weights_bilin_JRA.nc' , '' , '' 152 sn_qlw = 'rlds_JRA' , 3 , 'rlds' , .true. , .false. , 'yearly' , 'weights_bilin_JRA.nc' , '' , '' 153 sn_tair = 't10_JRA' , 3 , 'tas_10m' , .true. , .false. , 'yearly' , 'weights_bilin_JRA.nc' , '' , '' 154 sn_humi = 'q10_JRA' , 3 , 'huss_10m', .true. , .false. , 'yearly' , 'weights_bilin_JRA.nc' , '' , '' 155 sn_prec = 'precip_JRA' , 3 , 'prto' , .true. , .false. , 'yearly' , 'weights_bilin_JRA.nc' , '' , '' 156 sn_snow = 'snow_JRA' , 3 , 'prsn' , .true. , .false. , 'yearly' , 'weights_bilin_JRA.nc' , '' , '' 157 sn_slp = 'slp_JRA' , 3 , 'psl' , .true. , .false. , 'yearly' , 'weights_bilin_JRA.nc' , '' , '' 159 158 / 160 159 !----------------------------------------------------------------------- … … 201 200 ! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! 202 201 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! ! pairing ! filename ! 203 sn_rnf = ' runoff_WED025' , -1 , 'runoff' , .true. , .false., 'yearly' , '' , '' , ''202 sn_rnf = 'WED025_icb' , -1 , 'runoff' , .true. , .false., 'yearly' , '' , '' , '' 204 203 / 205 204 !----------------------------------------------------------------------- … … 221 220 cn_isfcav_mlt = '3eq' ! ice shelf melting formulation (spe/2eq/3eq/oasis) 222 221 ! ! spe = fwfisf is read from a forcing field 223 ! ! 2eq = ISOMIP like: 2 equations formulation (Hunter et al., 2006 )224 ! ! 3eq = ISOMIP+ like: 3 equations formulation (Asay-Davis et al., 201 5)222 ! ! 2eq = ISOMIP like: 2 equations formulation (Hunter et al., 2006 for a short description) 223 ! ! 3eq = ISOMIP+ like: 3 equations formulation (Asay-Davis et al., 2016 for a short description) 225 224 ! ! oasis = fwfisf is given by oasis and pattern by file sn_isfcav_fwf 226 225 ! ! cn_isfcav_mlt = 2eq or 3eq cases: 227 226 cn_gammablk = 'vel' ! scheme to compute gammat/s (spe,ad15,hj99) 228 ! ! ad15 = velocity dependend Gamma (u* * gammat/s) (Jenkins et al. 2010) 229 ! ! hj99 = velocity and stability dependent Gamma (Holland et al. 1999) 230 rn_gammat0 = 1.4e-2 ! gammat coefficient used in blk formula 231 rn_gammas0 = 4.e-4 ! gammas coefficient used in blk formula 227 ! ! spe = constant transfert velocity (rn_gammat0, rn_gammas0) 228 ! ! vel = velocity dependent transfert velocity (u* * gammat/s) (Asay-Davis et al. 2016 for a short description) 229 ! ! vel_stab = velocity and stability dependent transfert coeficient (Holland et al. 1999 for a complete description) 230 rn_gammat0 = 1.4e-2 ! gammat coefficient used in spe, vel and vel_stab gamma computation method 231 rn_gammas0 = 4.0e-4 ! gammas coefficient used in spe, vel and vel_stab gamma computation method 232 232 ! 233 233 rn_htbl = 30. ! thickness of the top boundary layer (Losh et al. 2008) … … 255 255 sn_isfpar_zmin = 'isfmlt_par', -12. , 'sozisfmin' , .false. , .true. , 'yearly' , '' , '' , '' 256 256 !* 'spe' and 'oasis' case 257 sn_isfpar_fwf = 'isfmlt_par' , -12. , 257 sn_isfpar_fwf = 'isfmlt_par' , -12. ,'sofwfisf' , .false. , .true. , 'yearly' , '' , '' , '' 258 258 !* 'bg03' case 259 sn_isfpar_Leff = 'isfmlt_par', 0. , 259 sn_isfpar_Leff = 'isfmlt_par', 0. ,'Leff' , .false. , .true. , 'yearly' , '' , '' , '' 260 260 ! 261 261 ! ---------------- ice sheet coupling ------------------------------- … … 300 300 ln_tide = .true. ! Activate tides 301 301 ln_tide_pot = .false. ! use tidal potential forcing 302 clname(1) = 'M2' ! name of constituent - all tidal components must be set in namelist_cfg303 clname(2) = 'S2'304 clname(3) = 'K1'305 clname(4) = 'O1'302 sn_tide_cnames(1) = 'M2' ! name of constituent - all tidal components must be set in namelist_cfg 303 sn_tide_cnames(2) = 'S2' 304 sn_tide_cnames(3) = 'K1' 305 sn_tide_cnames(4) = 'O1' 306 306 / 307 307 !----------------------------------------------------------------------- … … 340 340 ! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! 341 341 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! ! pairing ! filename ! 342 bn_ssh = ' bdyT_ssh_WED025' , -1 , 'sossheig' , .true. , .false., 'yearly' , '' , '' , ''343 bn_u2d = ' bdyU_u2d_WED025' , -1 , 'vobtcrtx' , .true. , .false., 'yearly' , '' , '' , ''344 bn_v2d = ' bdyV_u2d_WED025' , -1 , 'vobtcrty' , .true. , .false., 'yearly' , '' , '' , ''345 bn_u3d = ' bdyU_u3d_WED025' , -1 , 'vozocrtx' , .true. , .false., 'yearly' , '' , '' , ''346 bn_v3d = ' bdyV_u3d_WED025' , -1 , 'vomecrty' , .true. , .false., 'yearly' , '' , '' , ''347 bn_tem = ' bdyT_tra_WED025' , -1 , 'votemper' , .true. , .false., 'yearly' , '' , '' , ''348 bn_sal = ' bdyT_tra_WED025' , -1 , 'vosaline' , .true. , .false., 'yearly' , '' , '' , ''342 bn_ssh = 'WED025_bdyT_ssh' , -1 , 'sossheig' , .true. , .false., 'yearly' , '' , '' , '' 343 bn_u2d = 'WED025_bdyU_u2d' , -1 , 'vobtcrtx' , .true. , .false., 'yearly' , '' , '' , '' 344 bn_v2d = 'WED025_bdyV_u2d' , -1 , 'vobtcrty' , .true. , .false., 'yearly' , '' , '' , '' 345 bn_u3d = 'WED025_bdyU_u3d' , -1 , 'vozocrtx' , .true. , .false., 'yearly' , '' , '' , '' 346 bn_v3d = 'WED025_bdyV_u3d' , -1 , 'vomecrty' , .true. , .false., 'yearly' , '' , '' , '' 347 bn_tem = 'WED025_bdyT_tra' , -1 , 'votemper' , .true. , .false., 'yearly' , '' , '' , '' 348 bn_sal = 'WED025_bdyT_tra' , -1 , 'vosaline' , .true. , .false., 'yearly' , '' , '' , '' 349 349 !* for si3 350 bn_a_i = ' bdyT_ice_WED025' , -1 , 'ileadfra' , .true. , .false., 'yearly' , '' , '' , ''351 bn_h_i = ' bdyT_ice_WED025' , -1 , 'iicethic' , .true. , .false., 'yearly' , '' , '' , ''352 bn_h_s = ' bdyT_ice_WED025' , -1 , 'isnowthi' , .true. , .false., 'yearly' , '' , '' , ''350 bn_a_i = 'WED025_bdyT_ice' , -1 , 'ileadfra' , .true. , .false., 'yearly' , '' , '' , '' 351 bn_h_i = 'WED025_bdyT_ice' , -1 , 'iicethic' , .true. , .false., 'yearly' , '' , '' , '' 352 bn_h_s = 'WED025_bdyT_ice' , -1 , 'isnowthi' , .true. , .false., 'yearly' , '' , '' , '' 353 353 / 354 354 !----------------------------------------------------------------------- 355 355 &nambdy_tide ! tidal forcing at open boundaries (default: OFF) 356 356 !----------------------------------------------------------------------- 357 filtide = ' bdytide_WED025_' ! file name root of tidal forcing files357 filtide = 'WED025_bdytide_' ! file name root of tidal forcing files 358 358 / 359 359 … … 658 658 &namctl ! Control prints (default: OFF) 659 659 !----------------------------------------------------------------------- 660 ln_ctl = .FALSE. ! Toggle all report printing on/off (T/F); Ignored if sn_cfctl%l_config is T 661 sn_cfctl%l_config = .TRUE. ! IF .true. then control which reports are written with the following 662 sn_cfctl%l_runstat = .FALSE. ! switches and which areas produce reports with the proc integer settings. 663 sn_cfctl%l_trcstat = .FALSE. ! The default settings for the proc integers should ensure 664 sn_cfctl%l_oceout = .FALSE. ! that all areas report. 665 sn_cfctl%l_layout = .FALSE. ! 666 sn_cfctl%l_mppout = .FALSE. ! 667 sn_cfctl%l_mpptop = .FALSE. ! 668 sn_cfctl%procmin = 0 ! Minimum area number for reporting [default:0] 669 sn_cfctl%procmax = 1000000 ! Maximum area number for reporting [default:1000000] 670 sn_cfctl%procincr = 1 ! Increment for optional subsetting of areas [default:1] 671 sn_cfctl%ptimincr = 1 ! Timestep increment for writing time step progress info 672 nn_print = 0 ! level of print (0 no extra print) 673 nn_ictls = 0 ! start i indice of control sum (use to compare mono versus 674 nn_ictle = 0 ! end i indice of control sum multi processor runs 675 nn_jctls = 0 ! start j indice of control over a subdomain) 676 nn_jctle = 0 ! end j indice of control 677 nn_isplt = 1 ! number of processors in i-direction 678 nn_jsplt = 1 ! number of processors in j-direction 679 ln_timing = .true. ! timing by routine write out in timing.output file 680 ln_diacfl = .false. ! CFL diagnostics write out in cfl_diagnostics.ascii 660 sn_cfctl%l_runstat = .true. ! switches and which areas produce reports with the proc integer settings. 661 ln_timing = .true. ! timing by routine write out in timing.output file 681 662 / 682 663 !----------------------------------------------------------------------- -
NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/cfgs/WED025/EXPREF/namelist_ice_cfg
r11487 r13193 26 26 &namitd ! Ice discretization 27 27 !------------------------------------------------------------------------------ 28 ln_cat_hfn = .true. ! ice categories are defined by a function following rn_himean**(-0.05) 29 rn_himean = 2.0 ! expected domain-average ice thickness (m) 30 rn_himin = 0.01 ! minimum ice thickness (m) used in remapping 28 31 / 29 32 !------------------------------------------------------------------------------ 30 33 &namdyn ! Ice dynamics 31 34 !------------------------------------------------------------------------------ 35 ln_landfast_L16 = .true. ! landfast: parameterization from Lemieux 2016 32 36 / 33 37 !------------------------------------------------------------------------------ … … 42 46 &namdyn_adv ! Ice advection 43 47 !------------------------------------------------------------------------------ 48 ln_adv_Pra = .false. ! Advection scheme (Prather) 49 ln_adv_UMx = .true. ! Advection scheme (Ultimate-Macho) 50 nn_UMx = 5 ! order of the scheme for UMx (1-5 ; 20=centered 2nd order) 44 51 / 45 52 !------------------------------------------------------------------------------ … … 62 69 &namthd_do ! Ice growth in open water 63 70 !------------------------------------------------------------------------------ 71 rn_hinew = 0.02 ! thickness for new ice formation in open water (m), must be larger than rn_himin 72 ln_frazil = .true. ! Frazil ice parameterization (ice collection as a function of wind) 64 73 / 65 74 !------------------------------------------------------------------------------ … … 70 79 &namthd_pnd ! Melt ponds 71 80 !------------------------------------------------------------------------------ 81 ln_pnd = .true. ! activate melt ponds or not 82 ln_pnd_H12 = .true. ! activate evolutive melt ponds (from Holland et al 2012) 83 ln_pnd_alb = .true. ! melt ponds affect albedo or not 72 84 / 85 73 86 !------------------------------------------------------------------------------ 74 87 &namini ! Ice initialization 75 88 !------------------------------------------------------------------------------ 89 ln_iceini = .true. ! activate ice initialization (T) or not (F) 90 ln_iceini_file = .true. ! netcdf file provided for initialization (T) or not (F) 91 ! -- for ln_iceini_file = T 92 sn_hti = 'WED025_init_JRA_200001.nc', -12 ,'icethic_cea', .false. , .true., 'yearly' , '' , '', '' 93 sn_hts = 'WED025_init_JRA_200001.nc', -12 ,'icesnow_cea', .false. , .true., 'yearly' , '' , '', '' 94 sn_ati = 'WED025_init_JRA_200001.nc', -12 ,'ice_cover' , .false. , .true., 'yearly' , '' , '', '' 95 sn_smi = 'NOT USED' , -12 ,'smi' , .false. , .true., 'yearly' , '' , '', '' 96 sn_tmi = 'NOT USED' , -12 ,'tmi' , .false. , .true., 'yearly' , '' , '', '' 97 sn_tsu = 'NOT USED' , -12 ,'tsu' , .false. , .true., 'yearly' , '' , '', '' 98 sn_tms = 'NOT USED' , -12 ,'tms' , .false. , .true., 'yearly' , '' , '', '' 99 ! melt ponds (be careful, sn_apd is the pond concentration (not fraction), so it differs from rn_apd) 100 sn_apd = 'NOT USED' , -12 ,'apd' , .false. , .true., 'yearly' , '' , '', '' 101 sn_hpd = 'NOT USED' , -12 ,'hpd' , .false. , .true., 'yearly' , '' , '', '' 102 cn_dir='./' 76 103 / 77 104 !------------------------------------------------------------------------------ -
NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/src/ABL/ablmod.F90
r12724 r13193 592 592 ! !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 593 593 ! ! 8 *** Swap time indices for the next timestep 594 ! !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 595 nt_n = 1 + MOD( kt, 2)596 nt_a = 1 + MOD( kt+1, 2)597 ! 594 ! !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 595 nt_n = 1 + MOD( nt_n, 2) 596 nt_a = 1 + MOD( nt_a, 2) 597 ! 598 598 !--------------------------------------------------------------------------------------------------- 599 599 END SUBROUTINE abl_stp -
NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/src/ABL/par_abl.F90
r12724 r13193 29 29 LOGICAL , PUBLIC :: ln_smth_pblh !: smoothing of atmospheric PBL height 30 30 31 LOGICAL , PUBLIC :: ln_rstart_abl !: (de)activate abl restart 31 32 CHARACTER(len=256), PUBLIC :: cn_ablrst_in !: suffix of abl restart name (input) 32 33 CHARACTER(len=256), PUBLIC :: cn_ablrst_out !: suffix of abl restart name (output) -
NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/src/ABL/sbcabl.F90
r12724 r13193 68 68 LOGICAL :: lluldl 69 69 NAMELIST/namsbc_abl/ cn_dir, cn_dom, cn_ablrst_in, cn_ablrst_out, & 70 & cn_ablrst_indir, cn_ablrst_outdir, 70 & cn_ablrst_indir, cn_ablrst_outdir, ln_rstart_abl, & 71 71 & ln_hpgls_frc, ln_geos_winds, nn_dyn_restore, & 72 72 & rn_ldyn_min , rn_ldyn_max, rn_ltra_min, rn_ltra_max, & … … 263 263 264 264 ! Initialize the time index for now time (nt_n) and after time (nt_a) 265 nt_n = 1 + MOD( nit000 , 2) 266 nt_a = 1 + MOD( nit000+1, 2) 265 nt_n = 1; nt_a = 2 267 266 268 267 ! initialize ABL from data or restart 269 IF( ln_rstart ) THEN268 IF( ln_rstart_abl ) THEN 270 269 CALL abl_rst_read 271 270 ELSE … … 288 287 ENDIF 289 288 290 rhoa(:,:) = rho_air( tq_abl(:,:,2,nt_n,jp_ta), tq_abl(:,:,2,nt_n,jp_qa), sf(jp_slp)%fnow(:,:,1) ) !!GS: rhoa must be (re)computed here here to avoid division by zero in blk_ice_1 (TBI)291 292 289 END SUBROUTINE sbc_abl_init 293 290 … … 329 326 CALL fld_read( kt, nn_fsbc, sf ) ! input fields provided at the current time-step 330 327 331 !!------------------------------------------------------------------------------------------- 332 !! 2 - Compute Cd x ||U||, Ch x ||U||, Ce x ||U||, and SSQ using now fields 333 !!------------------------------------------------------------------------------------------- 334 335 CALL blk_oce_1( kt, u_abl(:,:,2,nt_n ), v_abl(:,:,2,nt_n ), & ! <<= in 336 & tq_abl(:,:,2,nt_n,jp_ta), tq_abl(:,:,2,nt_n,jp_qa), & ! <<= in 337 & sf(jp_slp )%fnow(:,:,1) , sst_m, ssu_m, ssv_m , & ! <<= in 338 & sf(jp_qsr )%fnow(:,:,1) , sf(jp_qlw )%fnow(:,:,1) , & ! <<= in 339 & tsk_m, zssq, zcd_du, zsen, zevp ) ! =>> out 340 341 #if defined key_si3 342 CALL blk_ice_1( u_abl(:,:,2,nt_n ), v_abl(:,:,2,nt_n ), & ! <<= in 343 & tq_abl(:,:,2,nt_n,jp_ta), tq_abl(:,:,2,nt_n,jp_qa), & ! <<= in 344 & sf(jp_slp)%fnow(:,:,1) , u_ice, v_ice, tm_su , & ! <<= in 345 & pseni=zseni, pevpi=zevpi, pssqi=zssqi, pcd_dui=zcd_dui ) ! <<= out 346 #endif 347 348 !!------------------------------------------------------------------------------------------- 349 !! 3 - Advance ABL variables from now (n) to after (n+1) 350 !!------------------------------------------------------------------------------------------- 351 352 CALL abl_stp( kt, tsk_m, ssu_m, ssv_m, zssq, & ! <<= in 353 & sf(jp_wndi)%fnow(:,:,:), sf(jp_wndj)%fnow(:,:,:), & ! <<= in 354 & sf(jp_tair)%fnow(:,:,:), sf(jp_humi)%fnow(:,:,:), & ! <<= in 355 & sf(jp_slp )%fnow(:,:,1), & ! <<= in 356 & sf(jp_hpgi)%fnow(:,:,:), sf(jp_hpgj)%fnow(:,:,:), & ! <<= in 357 & zcd_du, zsen, zevp, & ! <=> in/out 358 & wndm, utau, vtau, taum & ! =>> out 359 #if defined key_si3 360 & , tm_su, u_ice, v_ice, zssqi, zcd_dui & ! <<= in 361 & , zseni, zevpi, wndm_ice, ato_i & ! <<= in 362 & , utau_ice, vtau_ice & ! =>> out 363 #endif 364 & ) 365 !!------------------------------------------------------------------------------------------- 366 !! 4 - Finalize flux computation using ABL variables at (n+1), nt_n corresponds to (n+1) since 367 !! time swap is done in abl_stp 368 !!------------------------------------------------------------------------------------------- 369 370 CALL blk_oce_2( tq_abl(:,:,2,nt_n,jp_ta), & 371 & sf(jp_qsr )%fnow(:,:,1) , sf(jp_qlw )%fnow(:,:,1), & 372 & sf(jp_prec)%fnow(:,:,1) , sf(jp_snow)%fnow(:,:,1), & 373 & tsk_m, zsen, zevp ) 374 375 CALL abl_rst_opn( kt ) ! Open abl restart file (if necessary) 376 IF( lrst_abl ) CALL abl_rst_write( kt ) ! -- abl restart file 377 378 #if defined key_si3 379 ! Avoid a USE abl in icesbc module 380 sf(jp_tair)%fnow(:,:,1) = tq_abl(:,:,2,nt_n,jp_ta); sf(jp_humi)%fnow(:,:,1) = tq_abl(:,:,2,nt_n,jp_qa) 381 #endif 328 IF( MOD( kt - 1, nn_fsbc ) == 0 ) THEN 329 330 !!------------------------------------------------------------------------------------------- 331 !! 2 - Compute Cd x ||U||, Ch x ||U||, Ce x ||U||, and SSQ using now fields 332 !!------------------------------------------------------------------------------------------- 333 334 CALL blk_oce_1( kt, u_abl(:,:,2,nt_n ), v_abl(:,:,2,nt_n ), & ! <<= in 335 & tq_abl(:,:,2,nt_n,jp_ta), tq_abl(:,:,2,nt_n,jp_qa), & ! <<= in 336 & sf(jp_slp )%fnow(:,:,1) , sst_m, ssu_m, ssv_m , & ! <<= in 337 & sf(jp_qsr )%fnow(:,:,1) , sf(jp_qlw )%fnow(:,:,1) , & ! <<= in 338 & tsk_m, zssq, zcd_du, zsen, zevp ) ! =>> out 339 340 #if defined key_si3 341 CALL blk_ice_1( u_abl(:,:,2,nt_n ), v_abl(:,:,2,nt_n ), & ! <<= in 342 & tq_abl(:,:,2,nt_n,jp_ta), tq_abl(:,:,2,nt_n,jp_qa), & ! <<= in 343 & sf(jp_slp)%fnow(:,:,1) , u_ice, v_ice, tm_su , & ! <<= in 344 & pseni=zseni, pevpi=zevpi, pssqi=zssqi, pcd_dui=zcd_dui ) ! <<= out 345 #endif 346 347 !!------------------------------------------------------------------------------------------- 348 !! 3 - Advance ABL variables from now (n) to after (n+1) 349 !!------------------------------------------------------------------------------------------- 350 351 CALL abl_stp( kt, tsk_m, ssu_m, ssv_m, zssq, & ! <<= in 352 & sf(jp_wndi)%fnow(:,:,:), sf(jp_wndj)%fnow(:,:,:), & ! <<= in 353 & sf(jp_tair)%fnow(:,:,:), sf(jp_humi)%fnow(:,:,:), & ! <<= in 354 & sf(jp_slp )%fnow(:,:,1), & ! <<= in 355 & sf(jp_hpgi)%fnow(:,:,:), sf(jp_hpgj)%fnow(:,:,:), & ! <<= in 356 & zcd_du, zsen, zevp, & ! <=> in/out 357 & wndm, utau, vtau, taum & ! =>> out 358 #if defined key_si3 359 & , tm_su, u_ice, v_ice, zssqi, zcd_dui & ! <<= in 360 & , zseni, zevpi, wndm_ice, ato_i & ! <<= in 361 & , utau_ice, vtau_ice & ! =>> out 362 #endif 363 & ) 364 !!------------------------------------------------------------------------------------------- 365 !! 4 - Finalize flux computation using ABL variables at (n+1), nt_n corresponds to (n+1) since 366 !! time swap is done in abl_stp 367 !!------------------------------------------------------------------------------------------- 368 369 CALL blk_oce_2( tq_abl(:,:,2,nt_n,jp_ta), & 370 & sf(jp_qsr )%fnow(:,:,1) , sf(jp_qlw )%fnow(:,:,1), & 371 & sf(jp_prec)%fnow(:,:,1) , sf(jp_snow)%fnow(:,:,1), & 372 & tsk_m, zsen, zevp ) 373 374 CALL abl_rst_opn( kt ) ! Open abl restart file (if necessary) 375 IF( lrst_abl ) CALL abl_rst_write( kt ) ! -- abl restart file 376 377 #if defined key_si3 378 ! Avoid a USE abl in icesbc module 379 sf(jp_tair)%fnow(:,:,1) = tq_abl(:,:,2,nt_n,jp_ta); sf(jp_humi)%fnow(:,:,1) = tq_abl(:,:,2,nt_n,jp_qa) 380 #endif 381 END IF 382 382 383 383 END SUBROUTINE sbc_abl -
NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/src/ICE/iceistate.F90
r12724 r13193 179 179 ! 180 180 ! -- mandatory fields -- ! 181 zht_i_ini(:,:) = si(jp_hti)%fnow(:,:,1) 182 zht_s_ini(:,:) = si(jp_hts)%fnow(:,:,1) 183 zat_i_ini(:,:) = si(jp_ati)%fnow(:,:,1) 181 zht_i_ini(:,:) = si(jp_hti)%fnow(:,:,1) * tmask(:,:,1) 182 zht_s_ini(:,:) = si(jp_hts)%fnow(:,:,1) * tmask(:,:,1) 183 zat_i_ini(:,:) = si(jp_ati)%fnow(:,:,1) * tmask(:,:,1) 184 184 185 185 ! -- optional fields -- ! … … 219 219 & si(jp_hpd)%fnow(:,:,1) = ( rn_hpd_ini_n * zswitch + rn_hpd_ini_s * (1._wp - zswitch) ) * tmask(:,:,1) 220 220 ! 221 zsm_i_ini(:,:) = si(jp_smi)%fnow(:,:,1) 222 ztm_i_ini(:,:) = si(jp_tmi)%fnow(:,:,1) 223 zt_su_ini(:,:) = si(jp_tsu)%fnow(:,:,1) 224 ztm_s_ini(:,:) = si(jp_tms)%fnow(:,:,1) 225 zapnd_ini(:,:) = si(jp_apd)%fnow(:,:,1) 226 zhpnd_ini(:,:) = si(jp_hpd)%fnow(:,:,1) 221 zsm_i_ini(:,:) = si(jp_smi)%fnow(:,:,1) * tmask(:,:,1) 222 ztm_i_ini(:,:) = si(jp_tmi)%fnow(:,:,1) * tmask(:,:,1) 223 zt_su_ini(:,:) = si(jp_tsu)%fnow(:,:,1) * tmask(:,:,1) 224 ztm_s_ini(:,:) = si(jp_tms)%fnow(:,:,1) * tmask(:,:,1) 225 zapnd_ini(:,:) = si(jp_apd)%fnow(:,:,1) * tmask(:,:,1) 226 zhpnd_ini(:,:) = si(jp_hpd)%fnow(:,:,1) * tmask(:,:,1) 227 227 ! 228 228 ! change the switch for the following -
NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/src/OCE/BDY/bdydta.F90
r12724 r13193 92 92 INTEGER :: jbdy, jfld, jstart, jend, ib, jl ! dummy loop indices 93 93 INTEGER :: ii, ij, ik, igrd, ipl ! local integers 94 INTEGER, DIMENSION(jpbgrd) :: ilen195 94 TYPE(OBC_DATA) , POINTER :: dta_alias ! short cut 96 95 TYPE(FLD), DIMENSION(:), POINTER :: bf_alias … … 117 116 END DO 118 117 ENDIF 119 IF( dta_bdy(jbdy)%lneed_dyn2d .AND. ASSOCIATED(dta_bdy(jbdy)%u2d) ) THEN ! no SIZE with a unassociated pointer118 IF( ASSOCIATED(dta_bdy(jbdy)%u2d) ) THEN ! no SIZE with a unassociated pointer. v2d and u2d can differ on subdomain 120 119 igrd = 2 121 DO ib = 1, SIZE(dta_bdy(jbdy)%u2d) ! u2d is used only on the rim except if ln_full_vel = T, see bdy_dta_init120 DO ib = 1, SIZE(dta_bdy(jbdy)%u2d) ! u2d is used either over the whole bdy or only on the rim 122 121 ii = idx_bdy(jbdy)%nbi(ib,igrd) 123 122 ij = idx_bdy(jbdy)%nbj(ib,igrd) 124 123 dta_bdy(jbdy)%u2d(ib) = uu_b(ii,ij,Kmm) * umask(ii,ij,1) 125 124 END DO 125 ENDIF 126 IF( ASSOCIATED(dta_bdy(jbdy)%v2d) ) THEN ! no SIZE with a unassociated pointer. v2d and u2d can differ on subdomain 126 127 igrd = 3 127 DO ib = 1, SIZE(dta_bdy(jbdy)%v2d) ! v2d is used only on the rim except if ln_full_vel = T, see bdy_dta_init128 DO ib = 1, SIZE(dta_bdy(jbdy)%v2d) ! v2d is used either over the whole bdy or only on the rim 128 129 ii = idx_bdy(jbdy)%nbi(ib,igrd) 129 130 ij = idx_bdy(jbdy)%nbj(ib,igrd) … … 211 212 ! 212 213 ! if runoff condition: change river flow we read (in m3/s) into barotropic velocity (m/s) 213 IF( cn_tra(jbdy) == 'runoff' .AND. TRIM(bf_alias(jp_bdyu2d)%clrootname) /= 'NOT USED' ) THEN ! runoff and we read u/v2d214 IF( cn_tra(jbdy) == 'runoff' ) THEN ! runoff 214 215 ! 215 igrd = 2 ! zonal flow (m3/s) to barotropic zonal velocity (m/s) 216 DO ib = 1, idx_bdy(jbdy)%nblen(igrd) 217 ii = idx_bdy(jbdy)%nbi(ib,igrd) 218 ij = idx_bdy(jbdy)%nbj(ib,igrd) 219 dta_alias%u2d(ib) = dta_alias%u2d(ib) / ( e2u(ii,ij) * hu_0(ii,ij) ) 220 END DO 221 igrd = 3 ! meridional flow (m3/s) to barotropic meridional velocity (m/s) 222 DO ib = 1, idx_bdy(jbdy)%nblen(igrd) 223 ii = idx_bdy(jbdy)%nbi(ib,igrd) 224 ij = idx_bdy(jbdy)%nbj(ib,igrd) 225 dta_alias%v2d(ib) = dta_alias%v2d(ib) / ( e1v(ii,ij) * hv_0(ii,ij) ) 226 END DO 216 IF( ASSOCIATED(dta_bdy(jbdy)%u2d) ) THEN ! no SIZE with a unassociated pointer. v2d and u2d can differ on subdomain 217 igrd = 2 ! zonal flow (m3/s) to barotropic zonal velocity (m/s) 218 DO ib = 1, SIZE(dta_alias%u2d) ! u2d is used either over the whole bdy or only on the rim 219 ii = idx_bdy(jbdy)%nbi(ib,igrd) 220 ij = idx_bdy(jbdy)%nbj(ib,igrd) 221 dta_alias%u2d(ib) = dta_alias%u2d(ib) / ( e2u(ii,ij) * hu_0(ii,ij) ) 222 END DO 223 ENDIF 224 IF( ASSOCIATED(dta_bdy(jbdy)%v2d) ) THEN ! no SIZE with a unassociated pointer. v2d and u2d can differ on subdomain 225 igrd = 3 ! meridional flow (m3/s) to barotropic meridional velocity (m/s) 226 DO ib = 1, SIZE(dta_alias%v2d) ! v2d is used either over the whole bdy or only on the rim 227 ii = idx_bdy(jbdy)%nbi(ib,igrd) 228 ij = idx_bdy(jbdy)%nbj(ib,igrd) 229 dta_alias%v2d(ib) = dta_alias%v2d(ib) / ( e1v(ii,ij) * hv_0(ii,ij) ) 230 END DO 231 ENDIF 227 232 ENDIF 228 233 229 234 ! tidal harmonic forcing ONLY: initialise arrays 230 235 IF( nn_dyn2d_dta(jbdy) == 2 ) THEN ! we did not read ssh, u/v2d 231 IF( dta_alias%lneed_ssh .AND.ASSOCIATED(dta_alias%ssh) ) dta_alias%ssh(:) = 0._wp232 IF( dta_alias%lneed_dyn2d .AND.ASSOCIATED(dta_alias%u2d) ) dta_alias%u2d(:) = 0._wp233 IF( dta_alias%lneed_dyn2d .AND.ASSOCIATED(dta_alias%v2d) ) dta_alias%v2d(:) = 0._wp236 IF( ASSOCIATED(dta_alias%ssh) ) dta_alias%ssh(:) = 0._wp 237 IF( ASSOCIATED(dta_alias%u2d) ) dta_alias%u2d(:) = 0._wp 238 IF( ASSOCIATED(dta_alias%v2d) ) dta_alias%v2d(:) = 0._wp 234 239 ENDIF 235 240 … … 333 338 DO jbdy = 1, nb_bdy ! Tidal component added in ts loop 334 339 IF ( nn_dyn2d_dta(jbdy) .GE. 2 ) THEN 335 IF( cn_dyn2d(jbdy) == 'frs' ) THEN ; ilen1(:)=idx_bdy(jbdy)%nblen(:) 336 ELSE ; ilen1(:)=idx_bdy(jbdy)%nblenrim(:) 337 ENDIF 338 IF ( dta_bdy(jbdy)%lneed_ssh ) dta_bdy_s(jbdy)%ssh(1:ilen1(1)) = dta_bdy(jbdy)%ssh(1:ilen1(1)) 339 IF ( dta_bdy(jbdy)%lneed_dyn2d ) dta_bdy_s(jbdy)%u2d(1:ilen1(2)) = dta_bdy(jbdy)%u2d(1:ilen1(2)) 340 IF ( dta_bdy(jbdy)%lneed_dyn2d ) dta_bdy_s(jbdy)%v2d(1:ilen1(3)) = dta_bdy(jbdy)%v2d(1:ilen1(3)) 340 IF( ASSOCIATED(dta_bdy(jbdy)%ssh) ) dta_bdy_s(jbdy)%ssh(:) = dta_bdy(jbdy)%ssh(:) 341 IF( ASSOCIATED(dta_bdy(jbdy)%u2d) ) dta_bdy_s(jbdy)%u2d(:) = dta_bdy(jbdy)%u2d(:) 342 IF( ASSOCIATED(dta_bdy(jbdy)%v2d) ) dta_bdy_s(jbdy)%v2d(:) = dta_bdy(jbdy)%v2d(:) 341 343 ENDIF 342 344 END DO 343 345 ELSE ! Add tides if not split-explicit free surface else this is done in ts loop 344 346 ! 345 ! BDY: use pt_offset=1.0 as applied at the end of the step and bdy_dta_tides is referenced at the middle of the step346 347 CALL bdy_dta_tides( kt=kt, pt_offset = 1._wp ) 347 348 ENDIF … … 351 352 ! 352 353 END SUBROUTINE bdy_dta 353 354 354 355 355 356 SUBROUTINE bdy_dta_init … … 383 384 LOGICAL :: llneed ! 384 385 LOGICAL :: llread ! 386 LOGICAL :: llfullbdy ! 385 387 TYPE(FLD_N), DIMENSION(1), TARGET :: bn_tem, bn_sal, bn_u3d, bn_v3d ! must be an array to be used with fld_fill 386 388 TYPE(FLD_N), DIMENSION(1), TARGET :: bn_ssh, bn_u2d, bn_v2d ! informations about the fields to be read … … 497 499 igrd = 2 ! U point 498 500 ipk = 1 ! surface data 499 llneed = dta_bdy(jbdy)%lneed_dyn2d ! dta_bdy(jbdy)% sshwill be needed501 llneed = dta_bdy(jbdy)%lneed_dyn2d ! dta_bdy(jbdy)%u2d will be needed 500 502 llread = .NOT. ln_full_vel .AND. MOD(nn_dyn2d_dta(jbdy),2) == 1 ! don't get u2d from u3d and read NetCDF file 501 503 bf_alias => bf(jp_bdyu2d,jbdy:jbdy) ! alias for u2d structure of bdy number jbdy 502 504 bn_alias => bn_u2d ! alias for u2d structure of nambdy_dta 503 IF( ln_full_vel ) THEN ; iszdim = idx_bdy(jbdy)%nblen(igrd) ! will be computed from u3d -> need on the full bdy 504 ELSE ; iszdim = idx_bdy(jbdy)%nblenrim(igrd) ! used only on the rim 505 llfullbdy = ln_full_vel .OR. cn_dyn2d(jbdy) == 'frs' ! need u2d over the whole bdy or only over the rim? 506 IF( llfullbdy ) THEN ; iszdim = idx_bdy(jbdy)%nblen(igrd) 507 ELSE ; iszdim = idx_bdy(jbdy)%nblenrim(igrd) 505 508 ENDIF 506 509 ENDIF … … 509 512 igrd = 3 ! V point 510 513 ipk = 1 ! surface data 511 llneed = dta_bdy(jbdy)%lneed_dyn2d ! dta_bdy(jbdy)% sshwill be needed514 llneed = dta_bdy(jbdy)%lneed_dyn2d ! dta_bdy(jbdy)%v2d will be needed 512 515 llread = .NOT. ln_full_vel .AND. MOD(nn_dyn2d_dta(jbdy),2) == 1 ! don't get v2d from v3d and read NetCDF file 513 516 bf_alias => bf(jp_bdyv2d,jbdy:jbdy) ! alias for v2d structure of bdy number jbdy 514 517 bn_alias => bn_v2d ! alias for v2d structure of nambdy_dta 515 IF( ln_full_vel ) THEN ; iszdim = idx_bdy(jbdy)%nblen(igrd) ! will be computed from v3d -> need on the full bdy 516 ELSE ; iszdim = idx_bdy(jbdy)%nblenrim(igrd) ! used only on the rim 518 llfullbdy = ln_full_vel .OR. cn_dyn2d(jbdy) == 'frs' ! need v2d over the whole bdy or only over the rim? 519 IF( llfullbdy ) THEN ; iszdim = idx_bdy(jbdy)%nblen(igrd) 520 ELSE ; iszdim = idx_bdy(jbdy)%nblenrim(igrd) 517 521 ENDIF 518 522 ENDIF -
NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/src/OCE/BDY/bdyini.F90
r12377 r13193 19 19 USE oce ! ocean dynamics and tracers variables 20 20 USE dom_oce ! ocean space and time domain 21 USE sbc_oce , ONLY: nn_ice 21 22 USE bdy_oce ! unstructured open boundary conditions 22 23 USE bdydta ! open boundary cond. setting (bdy_dta_init routine) 23 24 USE bdytides ! open boundary cond. setting (bdytide_init routine) 24 25 USE tide_mod, ONLY: ln_tide ! tidal forcing 25 USE phycst 26 USE phycst , ONLY: rday 26 27 ! 27 28 USE in_out_manager ! I/O units … … 315 316 316 317 dta_bdy(ib_bdy)%lneed_ice = cn_ice(ib_bdy) /= 'none' 318 319 IF( dta_bdy(ib_bdy)%lneed_ice .AND. nn_ice /= 2 ) THEN 320 WRITE(ctmp1,*) 'bdy number ', ib_bdy,', needs ice model but nn_ice = ', nn_ice 321 CALL ctl_stop( ctmp1 ) 322 ENDIF 317 323 318 324 IF( lwp .AND. dta_bdy(ib_bdy)%lneed_ice ) THEN -
NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/src/OCE/BDY/bdytides.F90
r12724 r13193 65 65 !! namelist variables 66 66 !!------------------- 67 CHARACTER(len=80) :: filtide ! :Filename root for tidal input files68 LOGICAL :: ln_bdytide_2ddta ! :If true, read 2d harmonic data67 CHARACTER(len=80) :: filtide ! Filename root for tidal input files 68 LOGICAL :: ln_bdytide_2ddta ! If true, read 2d harmonic data 69 69 !! 70 INTEGER :: ib_bdy, itide, ib ! :dummy loop indices71 INTEGER :: ii, ij ! :dummy loop indices70 INTEGER :: ib_bdy, itide, ib ! dummy loop indices 71 INTEGER :: ii, ij ! dummy loop indices 72 72 INTEGER :: inum, igrd 73 INTEGER , DIMENSION(3) :: ilen0 !: length of boundary data (from OBC arrays)73 INTEGER :: isz ! bdy data size 74 74 INTEGER :: ios ! Local integer output status for namelist read 75 75 INTEGER :: nbdy_rdstart, nbdy_loc 76 CHARACTER(LEN=50) :: cerrmsg ! :error string77 CHARACTER(len=80) :: clfile ! :full file name for tidal input file78 REAL(wp),ALLOCATABLE, DIMENSION(:,:,:) :: dta_read ! :work space to read in tidal harmonics data79 REAL(wp),ALLOCATABLE, DIMENSION(:,:) :: ztr, zti ! :" " " " " " " "76 CHARACTER(LEN=50) :: cerrmsg ! error string 77 CHARACTER(len=80) :: clfile ! full file name for tidal input file 78 REAL(wp),ALLOCATABLE, DIMENSION(:,:,:) :: dta_read ! work space to read in tidal harmonics data 79 REAL(wp),ALLOCATABLE, DIMENSION(:,:) :: ztr, zti ! " " " " " " " " 80 80 !! 81 TYPE(TIDES_DATA), POINTER :: td !: local short cut 81 TYPE(TIDES_DATA), POINTER :: td ! local short cut 82 TYPE( OBC_DATA), POINTER :: dta ! local short cut 82 83 !! 83 84 NAMELIST/nambdy_tide/filtide, ln_bdytide_2ddta … … 93 94 IF( nn_dyn2d_dta(ib_bdy) >= 2 ) THEN 94 95 ! 95 td => tides(ib_bdy) 96 96 td => tides(ib_bdy) 97 dta => dta_bdy(ib_bdy) 98 97 99 ! Namelist nambdy_tide : tidal harmonic forcing at open boundaries 98 100 filtide(:) = '' … … 130 132 IF(lwp) WRITE(numout,*) ' ' 131 133 132 ! Allocate space for tidal harmonics data - get size from OBC data arrays 134 ! Allocate space for tidal harmonics data - get size from BDY data arrays 135 ! Allocate also slow varying data in the case of time splitting: 136 ! Do it anyway because at this stage knowledge of free surface scheme is unknown 133 137 ! ----------------------------------------------------------------------- 134 135 ! JC: If FRS scheme is used, we assume that tidal is needed over the whole 136 ! relaxation area 137 IF( cn_dyn2d(ib_bdy) == 'frs' ) THEN ; ilen0(:) = idx_bdy(ib_bdy)%nblen (:) 138 ELSE ; ilen0(:) = idx_bdy(ib_bdy)%nblenrim(:) 139 ENDIF 140 141 ALLOCATE( td%ssh0( ilen0(1), nb_harmo, 2 ) ) 142 ALLOCATE( td%ssh ( ilen0(1), nb_harmo, 2 ) ) 143 144 ALLOCATE( td%u0( ilen0(2), nb_harmo, 2 ) ) 145 ALLOCATE( td%u ( ilen0(2), nb_harmo, 2 ) ) 146 147 ALLOCATE( td%v0( ilen0(3), nb_harmo, 2 ) ) 148 ALLOCATE( td%v ( ilen0(3), nb_harmo, 2 ) ) 149 150 td%ssh0(:,:,:) = 0._wp 151 td%ssh (:,:,:) = 0._wp 152 td%u0 (:,:,:) = 0._wp 153 td%u (:,:,:) = 0._wp 154 td%v0 (:,:,:) = 0._wp 155 td%v (:,:,:) = 0._wp 156 138 IF( ASSOCIATED(dta%ssh) ) THEN ! we use bdy ssh on this mpi subdomain 139 isz = SIZE(dta%ssh) 140 ALLOCATE( td%ssh0( isz, nb_harmo, 2 ), td%ssh( isz, nb_harmo, 2 ), dta_bdy_s(ib_bdy)%ssh( isz ) ) 141 dta_bdy_s(ib_bdy)%ssh(:) = 0._wp ! needed? 142 ENDIF 143 IF( ASSOCIATED(dta%u2d) ) THEN ! we use bdy u2d on this mpi subdomain 144 isz = SIZE(dta%u2d) 145 ALLOCATE( td%u0 ( isz, nb_harmo, 2 ), td%u ( isz, nb_harmo, 2 ), dta_bdy_s(ib_bdy)%u2d( isz ) ) 146 dta_bdy_s(ib_bdy)%u2d(:) = 0._wp ! needed? 147 ENDIF 148 IF( ASSOCIATED(dta%v2d) ) THEN ! we use bdy v2d on this mpi subdomain 149 isz = SIZE(dta%v2d) 150 ALLOCATE( td%v0 ( isz, nb_harmo, 2 ), td%v ( isz, nb_harmo, 2 ), dta_bdy_s(ib_bdy)%v2d( isz ) ) 151 dta_bdy_s(ib_bdy)%v2d(:) = 0._wp ! needed? 152 ENDIF 153 154 ! fill td%ssh0, td%u0, td%v0 155 ! ----------------------------------------------------------------------- 157 156 IF( ln_bdytide_2ddta ) THEN 157 ! 158 158 ! It is assumed that each data file contains all complex harmonic amplitudes 159 159 ! given on the global domain (ie global, jpiglo x jpjglo) … … 162 162 ! 163 163 ! SSH fields 164 clfile = TRIM(filtide)//'_grid_T.nc' 165 CALL iom_open( clfile , inum ) 166 igrd = 1 ! Everything is at T-points here 167 DO itide = 1, nb_harmo 168 CALL iom_get( inum, jpdom_autoglo, TRIM(tide_harmonics(itide)%cname_tide)//'_z1', ztr(:,:) ) 169 CALL iom_get( inum, jpdom_autoglo, TRIM(tide_harmonics(itide)%cname_tide)//'_z2', zti(:,:) ) 170 DO ib = 1, ilen0(igrd) 171 ii = idx_bdy(ib_bdy)%nbi(ib,igrd) 172 ij = idx_bdy(ib_bdy)%nbj(ib,igrd) 173 IF( ii == 1 .OR. ii == jpi .OR. ij == 1 .OR. ij == jpj ) CYCLE ! to remove? 174 td%ssh0(ib,itide,1) = ztr(ii,ij) 175 td%ssh0(ib,itide,2) = zti(ii,ij) 176 END DO 177 END DO 178 CALL iom_close( inum ) 164 IF( ASSOCIATED(dta%ssh) ) THEN ! we use bdy ssh on this mpi subdomain 165 clfile = TRIM(filtide)//'_grid_T.nc' 166 CALL iom_open( clfile , inum ) 167 igrd = 1 ! Everything is at T-points here 168 DO itide = 1, nb_harmo 169 CALL iom_get( inum, jpdom_autoglo, TRIM(tide_harmonics(itide)%cname_tide)//'_z1', ztr(:,:) ) 170 CALL iom_get( inum, jpdom_autoglo, TRIM(tide_harmonics(itide)%cname_tide)//'_z2', zti(:,:) ) 171 DO ib = 1, SIZE(dta%ssh) 172 ii = idx_bdy(ib_bdy)%nbi(ib,igrd) 173 ij = idx_bdy(ib_bdy)%nbj(ib,igrd) 174 td%ssh0(ib,itide,1) = ztr(ii,ij) 175 td%ssh0(ib,itide,2) = zti(ii,ij) 176 END DO 177 END DO 178 CALL iom_close( inum ) 179 ENDIF 179 180 ! 180 181 ! U fields 181 clfile = TRIM(filtide)//'_grid_U.nc' 182 CALL iom_open( clfile , inum ) 183 igrd = 2 ! Everything is at U-points here 184 DO itide = 1, nb_harmo 185 CALL iom_get ( inum, jpdom_autoglo, TRIM(tide_harmonics(itide)%cname_tide)//'_u1', ztr(:,:) ) 186 CALL iom_get ( inum, jpdom_autoglo, TRIM(tide_harmonics(itide)%cname_tide)//'_u2', zti(:,:) ) 187 DO ib = 1, ilen0(igrd) 188 ii = idx_bdy(ib_bdy)%nbi(ib,igrd) 189 ij = idx_bdy(ib_bdy)%nbj(ib,igrd) 190 IF( ii == 1 .OR. ii == jpi .OR. ij == 1 .OR. ij == jpj ) CYCLE ! to remove? 191 td%u0(ib,itide,1) = ztr(ii,ij) 192 td%u0(ib,itide,2) = zti(ii,ij) 193 END DO 194 END DO 195 CALL iom_close( inum ) 182 IF( ASSOCIATED(dta%u2d) ) THEN ! we use bdy u2d on this mpi subdomain 183 clfile = TRIM(filtide)//'_grid_U.nc' 184 CALL iom_open( clfile , inum ) 185 igrd = 2 ! Everything is at U-points here 186 DO itide = 1, nb_harmo 187 CALL iom_get ( inum, jpdom_autoglo, TRIM(tide_harmonics(itide)%cname_tide)//'_u1', ztr(:,:) ) 188 CALL iom_get ( inum, jpdom_autoglo, TRIM(tide_harmonics(itide)%cname_tide)//'_u2', zti(:,:) ) 189 DO ib = 1, SIZE(dta%u2d) 190 ii = idx_bdy(ib_bdy)%nbi(ib,igrd) 191 ij = idx_bdy(ib_bdy)%nbj(ib,igrd) 192 td%u0(ib,itide,1) = ztr(ii,ij) 193 td%u0(ib,itide,2) = zti(ii,ij) 194 END DO 195 END DO 196 CALL iom_close( inum ) 197 ENDIF 196 198 ! 197 199 ! V fields 198 clfile = TRIM(filtide)//'_grid_V.nc' 199 CALL iom_open( clfile , inum ) 200 igrd = 3 ! Everything is at V-points here 201 DO itide = 1, nb_harmo 202 CALL iom_get ( inum, jpdom_autoglo, TRIM(tide_harmonics(itide)%cname_tide)//'_v1', ztr(:,:) ) 203 CALL iom_get ( inum, jpdom_autoglo, TRIM(tide_harmonics(itide)%cname_tide)//'_v2', zti(:,:) ) 204 DO ib = 1, ilen0(igrd) 205 ii = idx_bdy(ib_bdy)%nbi(ib,igrd) 206 ij = idx_bdy(ib_bdy)%nbj(ib,igrd) 207 IF( ii == 1 .OR. ii == jpi .OR. ij == 1 .OR. ij == jpj ) CYCLE ! to remove? 208 td%v0(ib,itide,1) = ztr(ii,ij) 209 td%v0(ib,itide,2) = zti(ii,ij) 210 END DO 211 END DO 212 CALL iom_close( inum ) 200 IF( ASSOCIATED(dta%v2d) ) THEN ! we use bdy v2d on this mpi subdomain 201 clfile = TRIM(filtide)//'_grid_V.nc' 202 CALL iom_open( clfile , inum ) 203 igrd = 3 ! Everything is at V-points here 204 DO itide = 1, nb_harmo 205 CALL iom_get ( inum, jpdom_autoglo, TRIM(tide_harmonics(itide)%cname_tide)//'_v1', ztr(:,:) ) 206 CALL iom_get ( inum, jpdom_autoglo, TRIM(tide_harmonics(itide)%cname_tide)//'_v2', zti(:,:) ) 207 DO ib = 1, SIZE(dta%v2d) 208 ii = idx_bdy(ib_bdy)%nbi(ib,igrd) 209 ij = idx_bdy(ib_bdy)%nbj(ib,igrd) 210 td%v0(ib,itide,1) = ztr(ii,ij) 211 td%v0(ib,itide,2) = zti(ii,ij) 212 END DO 213 END DO 214 CALL iom_close( inum ) 215 ENDIF 213 216 ! 214 217 DEALLOCATE( ztr, zti ) … … 218 221 ! Read tidal data only on bdy segments 219 222 ! 220 ALLOCATE( dta_read( MAXVAL( ilen0(1:3)), 1, 1 ) )223 ALLOCATE( dta_read( MAXVAL( idx_bdy(ib_bdy)%nblen(:) ), 1, 1 ) ) 221 224 ! 222 225 ! Open files and read in tidal forcing data … … 225 228 DO itide = 1, nb_harmo 226 229 ! ! SSH fields 227 clfile = TRIM(filtide)//TRIM(tide_harmonics(itide)%cname_tide)//'_grid_T.nc' 228 CALL iom_open( clfile, inum ) 229 CALL fld_map( inum, 'z1' , dta_read(1:ilen0(1),1:1,1:1) , 1, idx_bdy(ib_bdy)%nbmap(:,1) ) 230 td%ssh0(:,itide,1) = dta_read(1:ilen0(1),1,1) 231 CALL fld_map( inum, 'z2' , dta_read(1:ilen0(1),1:1,1:1) , 1, idx_bdy(ib_bdy)%nbmap(:,1) ) 232 td%ssh0(:,itide,2) = dta_read(1:ilen0(1),1,1) 233 CALL iom_close( inum ) 230 IF( ASSOCIATED(dta%ssh) ) THEN ! we use bdy ssh on this mpi subdomain 231 isz = SIZE(dta%ssh) 232 clfile = TRIM(filtide)//TRIM(tide_harmonics(itide)%cname_tide)//'_grid_T.nc' 233 CALL iom_open( clfile, inum ) 234 CALL fld_map( inum, 'z1', dta_read(1:isz,1:1,1:1) , 1, idx_bdy(ib_bdy)%nbmap(:,1) ) 235 td%ssh0(:,itide,1) = dta_read(1:isz,1,1) 236 CALL fld_map( inum, 'z2', dta_read(1:isz,1:1,1:1) , 1, idx_bdy(ib_bdy)%nbmap(:,1) ) 237 td%ssh0(:,itide,2) = dta_read(1:isz,1,1) 238 CALL iom_close( inum ) 239 ENDIF 234 240 ! ! U fields 235 clfile = TRIM(filtide)//TRIM(tide_harmonics(itide)%cname_tide)//'_grid_U.nc' 236 CALL iom_open( clfile, inum ) 237 CALL fld_map( inum, 'u1' , dta_read(1:ilen0(2),1:1,1:1) , 1, idx_bdy(ib_bdy)%nbmap(:,2) ) 238 td%u0(:,itide,1) = dta_read(1:ilen0(2),1,1) 239 CALL fld_map( inum, 'u2' , dta_read(1:ilen0(2),1:1,1:1) , 1, idx_bdy(ib_bdy)%nbmap(:,2) ) 240 td%u0(:,itide,2) = dta_read(1:ilen0(2),1,1) 241 CALL iom_close( inum ) 241 IF( ASSOCIATED(dta%u2d) ) THEN ! we use bdy u2d on this mpi subdomain 242 isz = SIZE(dta%u2d) 243 clfile = TRIM(filtide)//TRIM(tide_harmonics(itide)%cname_tide)//'_grid_U.nc' 244 CALL iom_open( clfile, inum ) 245 CALL fld_map( inum, 'u1', dta_read(1:isz,1:1,1:1) , 1, idx_bdy(ib_bdy)%nbmap(:,2) ) 246 td%u0(:,itide,1) = dta_read(1:isz,1,1) 247 CALL fld_map( inum, 'u2', dta_read(1:isz,1:1,1:1) , 1, idx_bdy(ib_bdy)%nbmap(:,2) ) 248 td%u0(:,itide,2) = dta_read(1:isz,1,1) 249 CALL iom_close( inum ) 250 ENDIF 242 251 ! ! V fields 243 clfile = TRIM(filtide)//TRIM(tide_harmonics(itide)%cname_tide)//'_grid_V.nc' 244 CALL iom_open( clfile, inum ) 245 CALL fld_map( inum, 'v1' , dta_read(1:ilen0(3),1:1,1:1) , 1, idx_bdy(ib_bdy)%nbmap(:,3) ) 246 td%v0(:,itide,1) = dta_read(1:ilen0(3),1,1) 247 CALL fld_map( inum, 'v2' , dta_read(1:ilen0(3),1:1,1:1) , 1, idx_bdy(ib_bdy)%nbmap(:,3) ) 248 td%v0(:,itide,2) = dta_read(1:ilen0(3),1,1) 249 CALL iom_close( inum ) 252 IF( ASSOCIATED(dta%v2d) ) THEN ! we use bdy v2d on this mpi subdomain 253 isz = SIZE(dta%v2d) 254 clfile = TRIM(filtide)//TRIM(tide_harmonics(itide)%cname_tide)//'_grid_V.nc' 255 CALL iom_open( clfile, inum ) 256 CALL fld_map( inum, 'v1', dta_read(1:isz,1:1,1:1) , 1, idx_bdy(ib_bdy)%nbmap(:,3) ) 257 td%v0(:,itide,1) = dta_read(1:isz,1,1) 258 CALL fld_map( inum, 'v2', dta_read(1:isz,1:1,1:1) , 1, idx_bdy(ib_bdy)%nbmap(:,3) ) 259 td%v0(:,itide,2) = dta_read(1:isz,1,1) 260 CALL iom_close( inum ) 261 ENDIF 250 262 ! 251 263 END DO ! end loop on tidal components … … 254 266 ! 255 267 ENDIF ! ln_bdytide_2ddta=.true. 256 !257 ! Allocate slow varying data in the case of time splitting:258 ! Do it anyway because at this stage knowledge of free surface scheme is unknown259 ALLOCATE( dta_bdy_s(ib_bdy)%ssh ( ilen0(1) ) )260 ALLOCATE( dta_bdy_s(ib_bdy)%u2d ( ilen0(2) ) )261 ALLOCATE( dta_bdy_s(ib_bdy)%v2d ( ilen0(3) ) )262 dta_bdy_s(ib_bdy)%ssh(:) = 0._wp263 dta_bdy_s(ib_bdy)%u2d(:) = 0._wp264 dta_bdy_s(ib_bdy)%v2d(:) = 0._wp265 268 ! 266 269 ENDIF ! nn_dyn2d_dta(ib_bdy) >= 2 … … 283 286 ! 284 287 LOGICAL :: lk_first_btstp ! =.TRUE. if time splitting and first barotropic step 285 INTEGER :: itide, ib_bdy, ib, igrd ! loop indices 286 INTEGER, DIMENSION(jpbgrd) :: ilen0 287 INTEGER, DIMENSION(1:jpbgrd) :: nblen, nblenrim ! short cuts 288 INTEGER :: itide, ib_bdy, ib ! loop indices 288 289 REAL(wp) :: z_arg, z_sarg, zramp, zoff, z_cost, z_sist, zt_offset 289 290 !!---------------------------------------------------------------------- … … 310 311 IF( nn_dyn2d_dta(ib_bdy) >= 2 ) THEN 311 312 ! 312 nblen(1:jpbgrd) = idx_bdy(ib_bdy)%nblen(1:jpbgrd)313 nblenrim(1:jpbgrd) = idx_bdy(ib_bdy)%nblenrim(1:jpbgrd)314 !315 IF( cn_dyn2d(ib_bdy) == 'frs' ) THEN ; ilen0(:) = nblen (:)316 ELSE ; ilen0(:) = nblenrim(:)317 ENDIF318 !319 313 ! We refresh nodal factors every day below 320 314 ! This should be done somewhere else … … 337 331 ! If time splitting, initialize arrays from slow varying open boundary data: 338 332 IF ( PRESENT(kit) ) THEN 339 IF ( dta_bdy(ib_bdy)%lneed_ssh ) dta_bdy(ib_bdy)%ssh(1:ilen0(1)) = dta_bdy_s(ib_bdy)%ssh(1:ilen0(1))340 IF ( dta_bdy(ib_bdy)%lneed_dyn2d ) dta_bdy(ib_bdy)%u2d(1:ilen0(2)) = dta_bdy_s(ib_bdy)%u2d(1:ilen0(2))341 IF ( dta_bdy(ib_bdy)%lneed_dyn2d ) dta_bdy(ib_bdy)%v2d(1:ilen0(3)) = dta_bdy_s(ib_bdy)%v2d(1:ilen0(3))333 IF ( ASSOCIATED(dta_bdy(ib_bdy)%ssh) ) dta_bdy(ib_bdy)%ssh(:) = dta_bdy_s(ib_bdy)%ssh(:) 334 IF ( ASSOCIATED(dta_bdy(ib_bdy)%u2d) ) dta_bdy(ib_bdy)%u2d(:) = dta_bdy_s(ib_bdy)%u2d(:) 335 IF ( ASSOCIATED(dta_bdy(ib_bdy)%v2d) ) dta_bdy(ib_bdy)%v2d(:) = dta_bdy_s(ib_bdy)%v2d(:) 342 336 ENDIF 343 337 ! … … 349 343 z_sist = zramp * SIN( z_sarg ) 350 344 ! 351 IF ( dta_bdy(ib_bdy)%lneed_ssh ) THEN 352 igrd=1 ! SSH on tracer grid 353 DO ib = 1, ilen0(igrd) 345 IF ( ASSOCIATED(dta_bdy(ib_bdy)%ssh) ) THEN ! SSH on tracer grid 346 DO ib = 1, SIZE(dta_bdy(ib_bdy)%ssh) 354 347 dta_bdy(ib_bdy)%ssh(ib) = dta_bdy(ib_bdy)%ssh(ib) + & 355 348 & ( tides(ib_bdy)%ssh(ib,itide,1)*z_cost + & … … 358 351 ENDIF 359 352 ! 360 IF ( dta_bdy(ib_bdy)%lneed_dyn2d ) THEN 361 igrd=2 ! U grid 362 DO ib = 1, ilen0(igrd) 353 IF ( ASSOCIATED(dta_bdy(ib_bdy)%u2d) ) THEN ! U grid 354 DO ib = 1, SIZE(dta_bdy(ib_bdy)%u2d) 363 355 dta_bdy(ib_bdy)%u2d(ib) = dta_bdy(ib_bdy)%u2d(ib) + & 364 356 & ( tides(ib_bdy)%u(ib,itide,1)*z_cost + & 365 357 & tides(ib_bdy)%u(ib,itide,2)*z_sist ) 366 358 END DO 367 igrd=3 ! V grid 368 DO ib = 1, ilen0(igrd) 359 ENDIF 360 ! 361 IF ( ASSOCIATED(dta_bdy(ib_bdy)%v2d) ) THEN ! V grid 362 DO ib = 1, SIZE(dta_bdy(ib_bdy)%v2d) 369 363 dta_bdy(ib_bdy)%v2d(ib) = dta_bdy(ib_bdy)%v2d(ib) + & 370 364 & ( tides(ib_bdy)%v(ib,itide,1)*z_cost + & … … 372 366 END DO 373 367 ENDIF 368 ! 374 369 END DO 375 END 370 ENDIF 376 371 END DO 377 372 ! … … 386 381 TYPE(TIDES_DATA), INTENT(inout) :: td ! tidal harmonics data 387 382 ! 388 INTEGER :: itide, igrd, ib ! dummy loop indices 389 INTEGER, DIMENSION(1) :: ilen0 ! length of boundary data (from OBC arrays) 383 INTEGER :: itide, isz, ib ! dummy loop indices 390 384 REAL(wp),ALLOCATABLE, DIMENSION(:) :: mod_tide, phi_tide 391 385 !!---------------------------------------------------------------------- 392 386 ! 393 igrd=1 394 ! SSH on tracer grid. 395 ilen0(1) = SIZE(td%ssh0(:,1,1)) 396 ! 397 ALLOCATE( mod_tide(ilen0(igrd)), phi_tide(ilen0(igrd)) ) 398 ! 399 DO itide = 1, nb_harmo 400 DO ib = 1, ilen0(igrd) 401 mod_tide(ib)=SQRT(td%ssh0(ib,itide,1)**2.+td%ssh0(ib,itide,2)**2.) 402 phi_tide(ib)=ATAN2(-td%ssh0(ib,itide,2),td%ssh0(ib,itide,1)) 387 IF( ASSOCIATED(td%ssh0) ) THEN ! SSH on tracer grid. 388 ! 389 isz = SIZE( td%ssh0, dim = 1 ) 390 ALLOCATE( mod_tide(isz), phi_tide(isz) ) 391 ! 392 DO itide = 1, nb_harmo 393 DO ib = 1, isz 394 mod_tide(ib)=SQRT( td%ssh0(ib,itide,1)*td%ssh0(ib,itide,1) + td%ssh0(ib,itide,2)*td%ssh0(ib,itide,2) ) 395 phi_tide(ib)=ATAN2(-td%ssh0(ib,itide,2),td%ssh0(ib,itide,1)) 396 END DO 397 DO ib = 1, isz 398 mod_tide(ib)=mod_tide(ib)*tide_harmonics(itide)%f 399 phi_tide(ib)=phi_tide(ib)+tide_harmonics(itide)%v0+tide_harmonics(itide)%u 400 END DO 401 DO ib = 1, isz 402 td%ssh(ib,itide,1)= mod_tide(ib)*COS(phi_tide(ib)) 403 td%ssh(ib,itide,2)=-mod_tide(ib)*SIN(phi_tide(ib)) 404 END DO 403 405 END DO 404 DO ib = 1 , ilen0(igrd) 405 mod_tide(ib)=mod_tide(ib)*tide_harmonics(itide)%f 406 phi_tide(ib)=phi_tide(ib)+tide_harmonics(itide)%v0+tide_harmonics(itide)%u 407 ENDDO 408 DO ib = 1 , ilen0(igrd) 409 td%ssh(ib,itide,1)= mod_tide(ib)*COS(phi_tide(ib)) 410 td%ssh(ib,itide,2)=-mod_tide(ib)*SIN(phi_tide(ib)) 411 ENDDO 412 END DO 413 ! 414 DEALLOCATE( mod_tide, phi_tide ) 406 ! 407 DEALLOCATE( mod_tide, phi_tide ) 408 ! 409 ENDIF 415 410 ! 416 411 END SUBROUTINE tide_init_elevation … … 424 419 TYPE(TIDES_DATA), INTENT(inout) :: td ! tidal harmonics data 425 420 ! 426 INTEGER :: itide, igrd, ib ! dummy loop indices 427 INTEGER, DIMENSION(3) :: ilen0 ! length of boundary data (from OBC arrays) 421 INTEGER :: itide, isz, ib ! dummy loop indices 428 422 REAL(wp),ALLOCATABLE, DIMENSION(:) :: mod_tide, phi_tide 429 423 !!---------------------------------------------------------------------- 430 424 ! 431 ilen0(2) = SIZE(td%u0(:,1,1)) 432 ilen0(3) = SIZE(td%v0(:,1,1)) 433 ! 434 igrd=2 ! U grid. 435 ! 436 ALLOCATE( mod_tide(ilen0(igrd)) , phi_tide(ilen0(igrd)) ) 437 ! 438 DO itide = 1, nb_harmo 439 DO ib = 1, ilen0(igrd) 440 mod_tide(ib)=SQRT(td%u0(ib,itide,1)**2.+td%u0(ib,itide,2)**2.) 441 phi_tide(ib)=ATAN2(-td%u0(ib,itide,2),td%u0(ib,itide,1)) 425 IF( ASSOCIATED(td%u0) ) THEN ! U grid. we use bdy u2d on this mpi subdomain 426 ! 427 isz = SIZE( td%u0, dim = 1 ) 428 ALLOCATE( mod_tide(isz), phi_tide(isz) ) 429 ! 430 DO itide = 1, nb_harmo 431 DO ib = 1, isz 432 mod_tide(ib)=SQRT( td%u0(ib,itide,1)*td%u0(ib,itide,1) + td%u0(ib,itide,2)*td%u0(ib,itide,2) ) 433 phi_tide(ib)=ATAN2(-td%u0(ib,itide,2),td%u0(ib,itide,1)) 434 END DO 435 DO ib = 1, isz 436 mod_tide(ib)=mod_tide(ib)*tide_harmonics(itide)%f 437 phi_tide(ib)=phi_tide(ib)+tide_harmonics(itide)%v0 + tide_harmonics(itide)%u 438 END DO 439 DO ib = 1, isz 440 td%u(ib,itide,1)= mod_tide(ib)*COS(phi_tide(ib)) 441 td%u(ib,itide,2)=-mod_tide(ib)*SIN(phi_tide(ib)) 442 END DO 442 443 END DO 443 DO ib = 1, ilen0(igrd) 444 mod_tide(ib)=mod_tide(ib)*tide_harmonics(itide)%f 445 phi_tide(ib)=phi_tide(ib)+tide_harmonics(itide)%v0 + tide_harmonics(itide)%u 446 ENDDO 447 DO ib = 1, ilen0(igrd) 448 td%u(ib,itide,1)= mod_tide(ib)*COS(phi_tide(ib)) 449 td%u(ib,itide,2)=-mod_tide(ib)*SIN(phi_tide(ib)) 450 ENDDO 451 END DO 452 ! 453 DEALLOCATE( mod_tide , phi_tide ) 454 ! 455 igrd=3 ! V grid. 456 ! 457 ALLOCATE( mod_tide(ilen0(igrd)) , phi_tide(ilen0(igrd)) ) 458 459 DO itide = 1, nb_harmo 460 DO ib = 1, ilen0(igrd) 461 mod_tide(ib)=SQRT(td%v0(ib,itide,1)**2.+td%v0(ib,itide,2)**2.) 462 phi_tide(ib)=ATAN2(-td%v0(ib,itide,2),td%v0(ib,itide,1)) 444 ! 445 DEALLOCATE( mod_tide, phi_tide ) 446 ! 447 ENDIF 448 ! 449 IF( ASSOCIATED(td%v0) ) THEN ! V grid. we use bdy u2d on this mpi subdomain 450 ! 451 isz = SIZE( td%v0, dim = 1 ) 452 ALLOCATE( mod_tide(isz), phi_tide(isz) ) 453 ! 454 DO itide = 1, nb_harmo 455 DO ib = 1, isz 456 mod_tide(ib)=SQRT( td%v0(ib,itide,1)*td%v0(ib,itide,1) + td%v0(ib,itide,2)*td%v0(ib,itide,2) ) 457 phi_tide(ib)=ATAN2(-td%v0(ib,itide,2),td%v0(ib,itide,1)) 458 END DO 459 DO ib = 1, isz 460 mod_tide(ib)=mod_tide(ib)*tide_harmonics(itide)%f 461 phi_tide(ib)=phi_tide(ib)+tide_harmonics(itide)%v0 + tide_harmonics(itide)%u 462 END DO 463 DO ib = 1, isz 464 td%v(ib,itide,1)= mod_tide(ib)*COS(phi_tide(ib)) 465 td%v(ib,itide,2)=-mod_tide(ib)*SIN(phi_tide(ib)) 466 END DO 463 467 END DO 464 DO ib = 1, ilen0(igrd) 465 mod_tide(ib)=mod_tide(ib)*tide_harmonics(itide)%f 466 phi_tide(ib)=phi_tide(ib)+tide_harmonics(itide)%v0 + tide_harmonics(itide)%u 467 ENDDO 468 DO ib = 1, ilen0(igrd) 469 td%v(ib,itide,1)= mod_tide(ib)*COS(phi_tide(ib)) 470 td%v(ib,itide,2)=-mod_tide(ib)*SIN(phi_tide(ib)) 471 ENDDO 472 END DO 473 ! 474 DEALLOCATE( mod_tide, phi_tide ) 475 ! 476 END SUBROUTINE tide_init_velocities 468 ! 469 DEALLOCATE( mod_tide, phi_tide ) 470 ! 471 ENDIF 472 ! 473 END SUBROUTINE tide_init_velocities 477 474 478 475 !!====================================================================== -
NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/src/OCE/C1D/step_c1d.F90
r12680 r13193 27 27 PRIVATE 28 28 29 PUBLIC stp_c1d ! called by opa.F9029 PUBLIC stp_c1d ! called by nemogcm.F90 30 30 31 31 !!---------------------------------------------------------------------- … … 56 56 ! 57 57 INTEGER :: jk ! dummy loop indice 58 INTEGER :: indic ! error indicator if < 059 58 !! --------------------------------------------------------------------- 60 61 indic = 0 ! reset to no error condition62 59 IF( kstp == nit000 ) CALL iom_init( "nemo") ! iom_put initialization (must be done after nemo_init for AGRIF+XIOS+OASIS) 63 60 IF( kstp /= nit000 ) CALL day( kstp ) ! Calendar (day was already called at nit000 in day_init) … … 88 85 !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 89 86 CALL dia_wri( kstp, Nnn ) ! ocean model: outputs 90 IF( lk_diahth )CALL dia_hth( kstp, Nnn ) ! Thermocline depth (20°C)87 CALL dia_hth( kstp, Nnn ) ! Thermocline depth (20°C) 91 88 92 89 … … 111 108 CALL eos( ts(:,:,:,:,Nnn), rhd, rhop, gdept_0(:,:,:) ) ! now potential density for zdfmxl 112 109 IF( ln_zdfnpc ) CALL tra_npc( kstp, Nnn, Nrhs, ts, Naa ) ! applied non penetrative convective adjustment on (t,s) 113 CALL tra_atf( kstp, Nbb, Nnn, Nrhs, Naa, ts ) ! time filtering of "now" tracer fields 114 115 110 CALL tra_atf( kstp, Nbb, Nnn, Naa, ts ) ! time filtering of "now" tracer arrays 116 111 117 112 !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> … … 139 134 ! Control and restarts 140 135 !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 141 CALL stp_ctl( kstp, Nnn , indic)136 CALL stp_ctl( kstp, Nnn ) 142 137 IF( kstp == nit000 ) CALL iom_close( numror ) ! close input ocean restart file 143 138 IF( lrst_oce ) CALL rst_write( kstp, Nbb, Nnn ) ! write output ocean restart file 144 139 ! 145 140 #if defined key_iomput 146 IF( kstp == nitend .OR. indic <0 ) CALL xios_context_finalize() ! needed for XIOS141 IF( kstp == nitend .OR. nstop > 0 ) CALL xios_context_finalize() ! needed for XIOS 147 142 ! 148 143 #endif -
NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/src/OCE/DIA/diaar5.F90
r13156 r13193 76 76 REAL(wp) :: zaw, zbw, zrw 77 77 ! 78 REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: zarea_ssh , zbotpres! 2D workspace78 REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: zarea_ssh, zbotpres ! 2D workspace 79 79 REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: z2d, zpe ! 2D workspace 80 REAL(wp), ALLOCATABLE, DIMENSION(:,:,:) :: z3d, zrhd , zrhop, ztpot, zgdept ! 3D workspace (zgdept: needed to use the substitute)80 REAL(wp), ALLOCATABLE, DIMENSION(:,:,:) :: z3d, zrhd, ztpot, zgdept ! 3D workspace (zgdept: needed to use the substitute) 81 81 REAL(wp), ALLOCATABLE, DIMENSION(:,:,:,:) :: ztsn ! 4D workspace 82 82 … … 88 88 IF( l_ar5 ) THEN 89 89 ALLOCATE( zarea_ssh(jpi,jpj), zbotpres(jpi,jpj), z2d(jpi,jpj) ) 90 ALLOCATE( zrhd(jpi,jpj,jpk) , zrhop(jpi,jpj,jpk))90 ALLOCATE( zrhd(jpi,jpj,jpk) ) 91 91 ALLOCATE( ztsn(jpi,jpj,jpk,jpts) ) 92 92 zarea_ssh(:,:) = e1e2t(:,:) * ssh(:,:,Kmm) … … 163 163 164 164 ! ! steric sea surface height 165 CALL eos( ts(:,:,:,:,Kmm), zrhd, zrhop, zgdept ) ! now in situ and potential density166 zrhop(:,:,jpk) = 0._wp167 CALL iom_put( 'rhop', zrhop )168 !169 165 zbotpres(:,:) = 0._wp ! no atmospheric surface pressure, levitating sea-ice 170 166 DO jk = 1, jpkm1 171 zbotpres(:,:) = zbotpres(:,:) + e3t(:,:,jk,Kmm) * zrhd(:,:,jk)167 zbotpres(:,:) = zbotpres(:,:) + e3t(:,:,jk,Kmm) * rhd(:,:,jk) 172 168 END DO 173 169 IF( ln_linssh ) THEN … … 176 172 DO jj = 1,jpj 177 173 iks = mikt(ji,jj) 178 zbotpres(ji,jj) = zbotpres(ji,jj) + ssh(ji,jj,Kmm) * zrhd(ji,jj,iks) + riceload(ji,jj)174 zbotpres(ji,jj) = zbotpres(ji,jj) + ssh(ji,jj,Kmm) * rhd(ji,jj,iks) + riceload(ji,jj) 179 175 END DO 180 176 END DO 181 177 ELSE 182 zbotpres(:,:) = zbotpres(:,:) + ssh(:,:,Kmm) * zrhd(:,:,1)178 zbotpres(:,:) = zbotpres(:,:) + ssh(:,:,Kmm) * rhd(:,:,1) 183 179 END IF 184 180 END IF … … 303 299 IF( l_ar5 ) THEN 304 300 DEALLOCATE( zarea_ssh , zbotpres, z2d ) 305 DEALLOCATE( zrhd , zrhop )306 301 DEALLOCATE( ztsn ) 307 302 ENDIF … … 377 372 IF( iom_use( 'voltot' ) .OR. iom_use( 'sshtot' ) .OR. iom_use( 'sshdyn' ) .OR. & 378 373 & iom_use( 'masstot' ) .OR. iom_use( 'temptot' ) .OR. iom_use( 'saltot' ) .OR. & 379 & iom_use( 'botpres' ) .OR. iom_use( 'sshthster' ) .OR. iom_use( 'sshsteric' ) ) L_ar5 = .TRUE. 374 & iom_use( 'botpres' ) .OR. iom_use( 'sshthster' ) .OR. iom_use( 'sshsteric' ) .OR. & 375 & iom_use( 'rhop' ) ) L_ar5 = .TRUE. 380 376 381 377 IF( l_ar5 ) THEN -
NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/src/OCE/DIA/diawri.F90
r12911 r13193 191 191 CALL iom_put( "sbs", z2d ) ! bottom salinity 192 192 ENDIF 193 194 CALL iom_put( "rhop", rhop(:,:,:) ) ! 3D potential density (sigma0) 193 195 194 196 IF ( iom_use("taubot") ) THEN ! bottom stress … … 1033 1035 CALL iom_close( inum ) 1034 1036 ENDIF 1037 ! 1035 1038 #endif 1036 1037 1039 END SUBROUTINE dia_wri_state 1038 1040 -
NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/src/OCE/DOM/dom_oce.F90
r12731 r13193 17 17 !!---------------------------------------------------------------------- 18 18 !! Agrif_Root : dummy function used when lk_agrif=F 19 !! Agrif_Fixed : dummy function used when lk_agrif=F 19 20 !! Agrif_CFixed : dummy function used when lk_agrif=F 20 21 !! dom_oce_alloc : dynamical allocation of dom_oce arrays … … 243 244 END FUNCTION Agrif_Root 244 245 246 INTEGER FUNCTION Agrif_Fixed() 247 Agrif_Fixed = 0 248 END FUNCTION Agrif_Fixed 249 245 250 CHARACTER(len=3) FUNCTION Agrif_CFixed() 246 251 Agrif_CFixed = '0' -
NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/src/OCE/DOM/dommsk.F90
r12680 r13193 261 261 ENDIF 262 262 END DO 263 #if defined key_agrif264 IF( .NOT. AGRIF_Root() ) THEN265 IF ((nbondi == 1).OR.(nbondi == 2)) fmask(nlci-1 , : ,jk) = 0.e0 ! east266 IF ((nbondi == -1).OR.(nbondi == 2)) fmask(1 , : ,jk) = 0.e0 ! west267 IF ((nbondj == 1).OR.(nbondj == 2)) fmask(: ,nlcj-1 ,jk) = 0.e0 ! north268 IF ((nbondj == -1).OR.(nbondj == 2)) fmask(: ,1 ,jk) = 0.e0 ! south269 ENDIF270 #endif271 263 END DO 272 264 ! -
NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/src/OCE/DOM/domvvl.F90
r12724 r13193 915 915 e3t(:,:,:,Kbb) = e3t(:,:,:,Kmm) 916 916 917 DO ji = 1, jpi 918 DO jj = 1, jpj 919 IF ( ht_0(ji,jj) .LE. 0.0 .AND. NINT( ssmask(ji,jj) ) .EQ. 1) THEN 920 CALL ctl_stop( 'dom_vvl_rst: ht_0 must be positive at potentially wet points' ) 921 ENDIF 922 END DO 923 END DO 917 DO_2D_11_11 918 IF ( ht_0(ji,jj) .LE. 0.0 .AND. NINT( ssmask(ji,jj) ) .EQ. 1) THEN 919 CALL ctl_stop( 'dom_vvl_rst: ht_0 must be positive at potentially wet points' ) 920 ENDIF 921 END_2D 924 922 ! 925 923 ELSE -
NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/src/OCE/DOM/istate.F90
r12911 r13193 24 24 USE dom_oce ! ocean space and time domain 25 25 USE daymod ! calendar 26 USE divhor ! horizontal divergence (div_hor routine)27 26 USE dtatsd ! data temperature and salinity (dta_tsd routine) 28 27 USE dtauvd ! data: U & V current (dta_uvd routine) … … 126 125 uu (:,:,:,Kmm) = uu (:,:,:,Kbb) 127 126 vv (:,:,:,Kmm) = vv (:,:,:,Kbb) 128 hdiv(:,:,jpk) = 0._wp ! bottom divergence set one for 0 to zero at jpk level129 CALL div_hor( 0, Kbb, Kmm ) ! compute interior hdiv value130 !!gm hdiv(:,:,:) = 0._wp131 127 132 128 !!gm POTENTIAL BUG : -
NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/src/OCE/DYN/divhor.F90
r12590 r13193 85 85 END_3D 86 86 ! 87 #if defined key_agrif88 IF( .NOT. Agrif_Root() ) THEN89 IF( nbondi == -1 .OR. nbondi == 2 ) hdiv( 2 , : ,:) = 0._wp ! west90 IF( nbondi == 1 .OR. nbondi == 2 ) hdiv( nlci-1, : ,:) = 0._wp ! east91 IF( nbondj == -1 .OR. nbondj == 2 ) hdiv( : , 2 ,:) = 0._wp ! south92 IF( nbondj == 1 .OR. nbondj == 2 ) hdiv( : ,nlcj-1,:) = 0._wp ! north93 ENDIF94 #endif95 !96 87 IF( ln_rnf ) CALL sbc_rnf_div( hdiv, Kmm ) !== runoffs ==! (update hdiv field) 97 88 ! -
NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/src/OCE/DYN/dynldf_lap_blp.F90
r12606 r13193 75 75 DO_2D_01_01 76 76 ! ! ahm * e3 * curl (computed from 1 to jpim1/jpjm1) 77 !!gm open question here : e3f at before or now ? probably now... 78 !!gm note that ahmf has already been multiplied by fmask 79 zcur(ji-1,jj-1) = & 80 & ahmf(ji-1,jj-1,jk) * e3f(ji-1,jj-1,jk) * r1_e1e2f(ji-1,jj-1) & 81 & * ( e2v(ji ,jj-1) * pv(ji ,jj-1,jk) - e2v(ji-1,jj-1) * pv(ji-1,jj-1,jk) & 82 & - e1u(ji-1,jj ) * pu(ji-1,jj ,jk) + e1u(ji-1,jj-1) * pu(ji-1,jj-1,jk) ) 77 zcur(ji-1,jj-1) = ahmf(ji-1,jj-1,jk) * e3f(ji-1,jj-1,jk) * r1_e1e2f(ji-1,jj-1) & ! ahmf already * by fmask 78 & * ( e2v(ji ,jj-1) * pv(ji ,jj-1,jk) - e2v(ji-1,jj-1) * pv(ji-1,jj-1,jk) & 79 & - e1u(ji-1,jj ) * pu(ji-1,jj ,jk) + e1u(ji-1,jj-1) * pu(ji-1,jj-1,jk) ) 83 80 ! ! ahm * div (computed from 2 to jpi/jpj) 84 !!gm note that ahmt has already been multiplied by tmask 85 zdiv(ji,jj) = & 86 & ahmt(ji,jj,jk) * r1_e1e2t(ji,jj) / e3t(ji,jj,jk,Kbb) & 87 & * ( e2u(ji,jj)*e3u(ji,jj,jk,Kbb) * pu(ji,jj,jk) & 88 & - e2u(ji-1,jj)*e3u(ji-1,jj,jk,Kbb) * pu(ji-1,jj,jk) & 89 & + e1v(ji,jj)*e3v(ji,jj,jk,Kbb) * pv(ji,jj,jk) & 90 & - e1v(ji,jj-1)*e3v(ji,jj-1,jk,Kbb) * pv(ji,jj-1,jk) ) 81 zdiv(ji,jj) = ahmt(ji,jj,jk) * r1_e1e2t(ji,jj) / e3t(ji,jj,jk,Kbb) & ! ahmt already * by tmask 82 & * ( e2u(ji,jj)*e3u(ji,jj,jk,Kbb) * pu(ji,jj,jk) - e2u(ji-1,jj)*e3u(ji-1,jj,jk,Kbb) * pu(ji-1,jj,jk) & 83 & + e1v(ji,jj)*e3v(ji,jj,jk,Kbb) * pv(ji,jj,jk) - e1v(ji,jj-1)*e3v(ji,jj-1,jk,Kbb) * pv(ji,jj-1,jk) ) 91 84 END_2D 92 85 ! 93 86 DO_2D_00_00 94 pu_rhs(ji,jj,jk) = pu_rhs(ji,jj,jk) + zsign * ( & 95 & - ( zcur(ji ,jj) - zcur(ji,jj-1) ) * r1_e2u(ji,jj) & 96 & / e3u(ji,jj,jk,Kmm) & 97 & + ( zdiv(ji+1,jj) - zdiv(ji,jj ) ) * r1_e1u(ji,jj) ) 87 pu_rhs(ji,jj,jk) = pu_rhs(ji,jj,jk) + zsign * umask(ji,jj,jk) * ( & ! * by umask is mandatory for dyn_ldf_blp use 88 & - ( zcur(ji ,jj) - zcur(ji,jj-1) ) * r1_e2u(ji,jj) / e3u(ji,jj,jk,Kmm) & 89 & + ( zdiv(ji+1,jj) - zdiv(ji,jj ) ) * r1_e1u(ji,jj) ) 98 90 ! 99 pv_rhs(ji,jj,jk) = pv_rhs(ji,jj,jk) + zsign * ( & 100 & ( zcur(ji,jj ) - zcur(ji-1,jj) ) * r1_e1v(ji,jj) & 101 & / e3v(ji,jj,jk,Kmm) & 102 & + ( zdiv(ji,jj+1) - zdiv(ji ,jj) ) * r1_e2v(ji,jj) ) 91 pv_rhs(ji,jj,jk) = pv_rhs(ji,jj,jk) + zsign * vmask(ji,jj,jk) * ( & ! * by vmask is mandatory for dyn_ldf_blp use 92 & ( zcur(ji,jj ) - zcur(ji-1,jj) ) * r1_e1v(ji,jj) / e3v(ji,jj,jk,Kmm) & 93 & + ( zdiv(ji,jj+1) - zdiv(ji ,jj) ) * r1_e2v(ji,jj) ) 103 94 END_2D 104 95 ! ! =============== -
NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/src/OCE/DYN/dynvor.F90
r12590 r13193 820 820 DO_3D_10_10( 1, jpk ) 821 821 IF( tmask(ji,jj+1,jk) + tmask(ji+1,jj+1,jk) & 822 & + tmask(ji,jj ,jk) + tmask(ji+1,jj +1,jk) == 3._wp ) fmask(ji,jj,jk) = 1._wp822 & + tmask(ji,jj ,jk) + tmask(ji+1,jj ,jk) == 3._wp ) fmask(ji,jj,jk) = 1._wp 823 823 END_3D 824 824 ! -
NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/src/OCE/DYN/sshwzv.F90
r13167 r13193 215 215 ENDIF 216 216 ! 217 #if defined key_agrif 218 IF( .NOT. AGRIF_Root() ) THEN 219 IF ((nbondi == 1).OR.(nbondi == 2)) pww(nlci-1 , : ,:) = 0.e0 ! east 220 IF ((nbondi == -1).OR.(nbondi == 2)) pww(2 , : ,:) = 0.e0 ! west 221 IF ((nbondj == 1).OR.(nbondj == 2)) pww(: ,nlcj-1 ,:) = 0.e0 ! north 222 IF ((nbondj == -1).OR.(nbondj == 2)) pww(: ,2 ,:) = 0.e0 ! south 223 ENDIF 224 #endif 217 #if defined key_agrif 218 IF( .NOT. AGRIF_Root() ) THEN 219 ! Mask vertical velocity at first/last columns/row 220 ! inside computational domain (cosmetic) 221 ! --- West --- ! 222 DO ji = mi0(2), mi1(2) 223 DO jj = 1, jpj 224 pww(ji,jj,:) = 0._wp 225 ENDDO 226 ENDDO 227 ! 228 ! --- East --- ! 229 DO ji = mi0(jpiglo-1), mi1(jpiglo-1) 230 DO jj = 1, jpj 231 pww(ji,jj,:) = 0._wp 232 ENDDO 233 ENDDO 234 ! 235 ! --- South --- ! 236 DO jj = mj0(2), mj1(2) 237 DO ji = 1, jpi 238 pww(ji,jj,:) = 0._wp 239 ENDDO 240 ENDDO 241 ! 242 ! --- North --- ! 243 DO jj = mj0(jpjglo-1), mj1(jpjglo-1) 244 DO ji = 1, jpi 245 pww(ji,jj,:) = 0._wp 246 ENDDO 247 ENDDO 248 ENDIF 249 #endif 225 250 ! 226 251 IF( ln_timing ) CALL timing_stop('wzv') -
NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/src/OCE/ICB/icbrst.F90
r12724 r13193 188 188 ! 189 189 INTEGER :: jn ! dummy loop index 190 INTEGER :: idg ! number of digits 190 191 INTEGER :: ix_dim, iy_dim, ik_dim, in_dim 191 192 CHARACTER(len=256) :: cl_path 192 193 CHARACTER(len=256) :: cl_filename 193 CHARACTER(len=256) :: cl_kt 194 CHARACTER(len=8 ) :: cl_kt 195 CHARACTER(LEN=12 ) :: clfmt ! writing format 194 196 TYPE(iceberg), POINTER :: this 195 197 TYPE(point) , POINTER :: pt … … 211 213 ! file name 212 214 WRITE(cl_kt, '(i8.8)') kt 213 cl_filename = TRIM(cexper)//"_"// TRIM(ADJUSTL(cl_kt))//"_"//TRIM(cn_icbrst_out)215 cl_filename = TRIM(cexper)//"_"//cl_kt//"_"//TRIM(cn_icbrst_out) 214 216 IF( lk_mpp ) THEN 215 WRITE(cl_filename,'(A,"_",I4.4,".nc")') TRIM(cl_filename), narea-1 217 idg = MAX( INT(LOG10(REAL(MAX(1,jpnij-1),wp))) + 1, 4 ) ! how many digits to we need to write? min=4, max=9 218 WRITE(clfmt, "('(a,a,i', i1, '.', i1, ',a)')") idg, idg ! '(a,a,ix.x,a)' 219 WRITE(cl_filename, clfmt) TRIM(cl_filename), '_', narea-1, '.nc' 216 220 ELSE 217 WRITE(cl_filename,'( A,".nc")') TRIM(cl_filename)221 WRITE(cl_filename,'(a,a)') TRIM(cl_filename), '.nc' 218 222 ENDIF 219 223 -
NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/src/OCE/ICB/icbtrj.F90
r12724 r13193 62 62 ! 63 63 INTEGER :: iret, iyear, imonth, iday 64 INTEGER :: idg ! number of digits 64 65 REAL(wp) :: zfjulday, zsec 65 66 CHARACTER(len=80) :: cl_filename 66 CHARACTER(LEN=20) :: cldate_ini, cldate_end 67 CHARACTER(LEN=12) :: clfmt ! writing format 68 CHARACTER(LEN=8 ) :: cldate_ini, cldate_end 67 69 TYPE(iceberg), POINTER :: this 68 70 TYPE(point) , POINTER :: pt … … 80 82 81 83 ! define trajectory output name 82 IF ( lk_mpp ) THEN ; WRITE(cl_filename,'("trajectory_icebergs_",A,"-",A,"_",I4.4,".nc")') & 83 & TRIM(ADJUSTL(cldate_ini)), TRIM(ADJUSTL(cldate_end)), narea-1 84 ELSE ; WRITE(cl_filename,'("trajectory_icebergs_",A,"-",A ,".nc")') & 85 & TRIM(ADJUSTL(cldate_ini)), TRIM(ADJUSTL(cldate_end)) 84 cl_filename = 'trajectory_icebergs_'//cldate_ini//'-'//cldate_end 85 IF ( lk_mpp ) THEN 86 idg = MAX( INT(LOG10(REAL(MAX(1,jpnij-1),wp))) + 1, 4 ) ! how many digits to we need to write? min=4, max=9 87 WRITE(clfmt, "('(a,a,i', i1, '.', i1, ',a)')") idg, idg ! '(a,a,ix.x,a)' 88 WRITE(cl_filename, clfmt) TRIM(cl_filename), '_', narea-1, '.nc' 89 ELSE 90 WRITE(cl_filename,'(a,a)') TRIM(cl_filename), '.nc' 86 91 ENDIF 87 92 IF( lwp .AND. nn_verbose_level >= 0 ) WRITE(numout,'(2a)') 'icebergs, icb_trj_init: creating ',TRIM(cl_filename) -
NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/src/OCE/IOM/in_out_manager.F90
r12377 r13193 100 100 !!---------------------------------------------------------------------- 101 101 TYPE :: sn_ctl !: structure for control over output selection 102 LOGICAL :: l_glochk = .FALSE. !: range sanity checks are local (F) or global (T)103 ! Use global setting for debugging only;104 ! local breaches will still be reported105 ! and stop the code in most cases.106 LOGICAL :: l_allon = .FALSE. !: overall control; activate all following output options107 LOGICAL :: l_config = .FALSE. !: activate/deactivate finer control108 ! Note if l_config is True then sn_cfctl%l_allon is ignored.109 ! Otherwise setting sn_cfctl%l_allon T/F is equivalent to110 ! setting all the following logicals in this structure T/F111 ! and disabling subsetting of processors112 102 LOGICAL :: l_runstat = .FALSE. !: Produce/do not produce run.stat file (T/F) 113 103 LOGICAL :: l_trcstat = .FALSE. !: Produce/do not produce tracer.stat file (T/F) … … 169 159 INTEGER :: no_print = 0 !: optional argument of fld_fill (if present, suppress some control print) 170 160 INTEGER :: nstop = 0 !: error flag (=number of reason for a premature stop run) 161 !$AGRIF_DO_NOT_TREAT 162 INTEGER :: ngrdstop = -1 !: grid number having nstop > 1 163 !$AGRIF_END_DO_NOT_TREAT 171 164 INTEGER :: nwarn = 0 !: warning flag (=number of warning found during the run) 172 165 CHARACTER(lc) :: ctmp1, ctmp2, ctmp3 !: temporary characters 1 to 3 -
NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/src/OCE/IOM/iom_def.F90
r12724 r13193 33 33 INTEGER, PARAMETER, PUBLIC :: jpmax_vars = 1200 !: maximum number of variables in one file 34 34 INTEGER, PARAMETER, PUBLIC :: jpmax_dims = 4 !: maximum number of dimensions for one variable 35 INTEGER, PARAMETER, PUBLIC :: jpmax_digits = 5!: maximum number of digits for the cpu number in the file name35 INTEGER, PARAMETER, PUBLIC :: jpmax_digits = 9 !: maximum number of digits for the cpu number in the file name 36 36 37 37 -
NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/src/OCE/IOM/iom_nf90.F90
r12724 r13193 62 62 CHARACTER(LEN=256) :: clinfo ! info character 63 63 CHARACTER(LEN=256) :: cltmp ! temporary character 64 CHARACTER(LEN=12 ) :: clfmt ! writing format 64 65 CHARACTER(LEN=3 ) :: clcomp ! name of component calling iom_nf90_open 66 INTEGER :: idg ! number of digits 65 67 INTEGER :: iln ! lengths of character 66 68 INTEGER :: istop ! temporary storage of nstop … … 109 111 IF( ldwrt ) THEN !* the file should be open in write mode so we create it... 110 112 IF( jpnij > 1 ) THEN 111 WRITE(cltmp,'(a,a,i4.4,a)') cdname(1:iln-1), '_', narea-1, '.nc' 113 idg = MAX( INT(LOG10(REAL(MAX(1,jpnij-1),wp))) + 1, 4 ) ! how many digits to we need to write? min=4, max=9 114 WRITE(clfmt, "('(a,a,i', i1, '.', i1, ',a)')") idg, idg ! '(a,a,ix.x,a)' 115 WRITE(cltmp,clfmt) cdname(1:iln-1), '_', narea-1, '.nc' 112 116 cdname = TRIM(cltmp) 113 117 ENDIF -
NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/src/OCE/ISF/isfdiags.F90
r12616 r13193 89 89 REAL(wp), DIMENSION(jpi,jpj), INTENT(in) :: phtbl, pfrac ! thickness of the tbl and fraction of last cell affected by the tbl 90 90 REAL(wp), DIMENSION(jpi,jpj), INTENT(in) :: pvar2d ! 2d var to map in 3d 91 CHARACTER(LEN= 256), INTENT(in) :: cdvar91 CHARACTER(LEN=*), INTENT(in) :: cdvar 92 92 !!--------------------------------------------------------------------- 93 93 INTEGER :: ji, jj, jk ! loop indices -
NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/src/OCE/LBC/lib_mpp.F90
r12724 r13193 1112 1112 CHARACTER(len=*), INTENT(in ), OPTIONAL :: cd2, cd3, cd4, cd5 1113 1113 CHARACTER(len=*), INTENT(in ), OPTIONAL :: cd6, cd7, cd8, cd9, cd10 1114 ! 1115 CHARACTER(LEN=8) :: clfmt ! writing format 1116 INTEGER :: inum 1114 1117 !!---------------------------------------------------------------------- 1115 1118 ! 1116 1119 nstop = nstop + 1 1117 1120 ! 1118 ! force to open ocean.output file if not already opened 1119 IF( numout == 6 ) CALL ctl_opn( numout, 'ocean.output', 'APPEND', 'FORMATTED', 'SEQUENTIAL', -1, 6, .FALSE. ) 1121 IF( cd1 == 'STOP' .AND. narea /= 1 ) THEN ! Immediate stop: add an arror message in 'ocean.output' file 1122 CALL ctl_opn( inum, 'ocean.output', 'APPEND', 'FORMATTED', 'SEQUENTIAL', -1, 6, .FALSE. ) 1123 WRITE(inum,*) 1124 WRITE(inum,*) ' ==>>> Look for "E R R O R" messages in all existing *ocean.output* files' 1125 CLOSE(inum) 1126 ENDIF 1127 IF( numout == 6 ) THEN ! force to open ocean.output file if not already opened 1128 CALL ctl_opn( numout, 'ocean.output', 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, -1, .FALSE., narea ) 1129 ENDIF 1120 1130 ! 1121 1131 WRITE(numout,*) … … 1145 1155 WRITE(numout,*) 'huge E-R-R-O-R : immediate stop' 1146 1156 WRITE(numout,*) 1157 CALL FLUSH(numout) 1158 CALL SLEEP(60) ! make sure that all output and abort files are written by all cores. 60s should be enough... 1147 1159 CALL mppstop( ld_abort = .true. ) 1148 1160 ENDIF … … 1207 1219 ! 1208 1220 CHARACTER(len=80) :: clfile 1221 CHARACTER(LEN=10) :: clfmt ! writing format 1209 1222 INTEGER :: iost 1223 INTEGER :: idg ! number of digits 1210 1224 !!---------------------------------------------------------------------- 1211 1225 ! … … 1214 1228 clfile = TRIM(cdfile) 1215 1229 IF( PRESENT( karea ) ) THEN 1216 IF( karea > 1 ) WRITE(clfile, "(a,'_',i4.4)") TRIM(clfile), karea-1 1230 IF( karea > 1 ) THEN 1231 ! Warning: jpnij is maybe not already defined when calling ctl_opn -> use mppsize instead of jpnij 1232 idg = MAX( INT(LOG10(REAL(MAX(1,mppsize-1),wp))) + 1, 4 ) ! how many digits to we need to write? min=4, max=9 1233 WRITE(clfmt, "('(a,a,i', i1, '.', i1, ')')") idg, idg ! '(a,a,ix.x)' 1234 WRITE(clfile, clfmt) TRIM(clfile), '_', karea-1 1235 ENDIF 1217 1236 ENDIF 1218 1237 #if defined key_agrif -
NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/src/OCE/LBC/mpp_loc_generic.h90
r10716 r13193 32 32 REAL(wp) , INTENT( out) :: pmin ! Global minimum of ptab 33 33 INDEX_TYPE(:) ! index of minimum in global frame 34 # if defined key_mpp_mpi35 34 ! 36 35 INTEGER :: ierror, ii, idim … … 56 55 ! 57 56 kindex(1) = mig( ilocs(1) ) 58 # 57 #if defined DIM_2d || defined DIM_3d /* avoid warning when kindex has 1 element */ 59 58 kindex(2) = mjg( ilocs(2) ) 60 # 61 # 59 #endif 60 #if defined DIM_3d /* avoid warning when kindex has 2 elements */ 62 61 kindex(3) = ilocs(3) 63 # 62 #endif 64 63 ! 65 64 DEALLOCATE (ilocs) 66 65 ! 67 66 index0 = kindex(1)-1 ! 1d index starting at 0 68 # 67 #if defined DIM_2d || defined DIM_3d /* avoid warning when kindex has 1 element */ 69 68 index0 = index0 + jpiglo * (kindex(2)-1) 70 # 71 # 69 #endif 70 #if defined DIM_3d /* avoid warning when kindex has 2 elements */ 72 71 index0 = index0 + jpiglo * jpjglo * (kindex(3)-1) 73 # 72 #endif 74 73 END IF 75 74 zain(1,:) = zmin … … 77 76 ! 78 77 IF( ln_timing ) CALL tic_tac(.TRUE., ld_global = .TRUE.) 78 #if defined key_mpp_mpi 79 79 CALL MPI_ALLREDUCE( zain, zaout, 1, MPI_2DOUBLE_PRECISION, MPI_OPERATION ,MPI_COMM_OCE, ierror) 80 #else 81 zaout(:,:) = zain(:,:) 82 #endif 80 83 IF( ln_timing ) CALL tic_tac(.FALSE., ld_global = .TRUE.) 81 84 ! 82 85 pmin = zaout(1,1) 83 86 index0 = NINT( zaout(2,1) ) 84 # 87 #if defined DIM_3d /* avoid warning when kindex has 2 elements */ 85 88 kindex(3) = index0 / (jpiglo*jpjglo) 86 89 index0 = index0 - kindex(3) * (jpiglo*jpjglo) 87 # 88 # 90 #endif 91 #if defined DIM_2d || defined DIM_3d /* avoid warning when kindex has 1 element */ 89 92 kindex(2) = index0 / jpiglo 90 93 index0 = index0 - kindex(2) * jpiglo 91 # 94 #endif 92 95 kindex(1) = index0 93 96 kindex(:) = kindex(:) + 1 ! start indices at 1 94 #else95 kindex = 0 ; pmin = 0.96 WRITE(*,*) 'ROUTINE_LOC: You should not have seen this print! error?'97 #endif98 97 99 98 END SUBROUTINE ROUTINE_LOC -
NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/src/OCE/OBS/obs_grid.F90
r10068 r13193 684 684 & fhistx1, fhistx2, fhisty1, fhisty2 685 685 REAL(wp) :: histtol 686 686 CHARACTER(LEN=26) :: clfmt ! writing format 687 INTEGER :: idg ! number of digits 688 687 689 IF (ln_grid_search_lookup) THEN 688 690 … … 709 711 710 712 IF ( ln_grid_global ) THEN 711 WRITE(cfname, FMT="(A,'_',A)") & 712 & TRIM(cn_gridsearchfile), 'global.nc' 713 WRITE(cfname, FMT="(A,'_',A)") TRIM(cn_gridsearchfile), 'global.nc' 713 714 ELSE 714 WRITE(cfname, FMT="(A,'_',I4.4,'of',I4.4,'by',I4.4,'.nc')") & 715 & TRIM(cn_gridsearchfile), nproc, jpni, jpnj 715 idg = MAX( INT(LOG10(REAL(jpnij,wp))) + 1, 4 ) ! how many digits to we need to write? min=4, max=9 716 ! define the following format: "(a,a,ix.x,a,ix.x,a,ix.x,a)" 717 WRITE(clfmt, "('(a,a,i', i1, '.', i1',a,i', i1, '.', i1',a,i', i1, '.', i1',a)')") idg, idg, idg, idg, idg, idg 718 WRITE(cfname, clfmt ) TRIM(cn_gridsearchfile),'_', nproc,'of', jpni,'by', jpnj,'.nc' 716 719 ENDIF 717 720 -
NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/src/OCE/OBS/obs_write.F90
r12377 r13193 86 86 CHARACTER(LEN=40) :: clfname 87 87 CHARACTER(LEN=10) :: clfiletype 88 CHARACTER(LEN=12) :: clfmt ! writing format 89 INTEGER :: idg ! number of digits 88 90 INTEGER :: ilevel 89 91 INTEGER :: jvar … … 181 183 fbdata%caddname(1) = 'Hx' 182 184 183 WRITE(clfname, FMT="(A,'_fdbk_',I4.4,'.nc')") TRIM(clfiletype), nproc 185 idg = MAX( INT(LOG10(REAL(jpnij,wp))) + 1, 4 ) ! how many digits to we need to write? min=4, max=9 186 WRITE(clfmt, "('(a,a,i', i1, '.', i1, ',a)')") idg, idg ! '(a,a,ix.x,a)' 187 WRITE(clfname,clfmt) TRIM(clfiletype), '_fdbk_', nproc, '.nc' 184 188 185 189 IF(lwp) THEN … … 326 330 CHARACTER(LEN=10) :: clfiletype 327 331 CHARACTER(LEN=12), PARAMETER :: cpname = 'obs_wri_surf' 332 CHARACTER(LEN=12) :: clfmt ! writing format 333 INTEGER :: idg ! number of digits 328 334 INTEGER :: jo 329 335 INTEGER :: ja … … 453 459 fbdata%caddname(1) = 'Hx' 454 460 455 WRITE(clfname, FMT="(A,'_fdbk_',I4.4,'.nc')") TRIM(clfiletype), nproc 461 idg = MAX( INT(LOG10(REAL(jpnij,wp))) + 1, 4 ) ! how many digits to we need to write? min=4, max=9 462 WRITE(clfmt, "('(a,a,i', i1, '.', i1, ',a)')") idg, idg ! '(a,a,ix.x,a)' 463 WRITE(clfname,clfmt) TRIM(clfiletype), '_fdbk_', nproc, '.nc' 456 464 457 465 IF(lwp) THEN -
NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/src/OCE/SBC/sbcblk.F90
r12724 r13193 627 627 628 628 END SELECT 629 629 630 IF( iom_use('Cd_oce') ) CALL iom_put("Cd_oce", zcd_oce * tmask(:,:,1)) 631 IF( iom_use('Ce_oce') ) CALL iom_put("Ce_oce", zce_oce * tmask(:,:,1)) 632 IF( iom_use('Ch_oce') ) CALL iom_put("Ch_oce", zch_oce * tmask(:,:,1)) 633 !! LB: mainly here for debugging purpose: 634 IF( iom_use('theta_zt') ) CALL iom_put("theta_zt", (ztpot-rt0) * tmask(:,:,1)) ! potential temperature at z=zt 635 IF( iom_use('q_zt') ) CALL iom_put("q_zt", zqair * tmask(:,:,1)) ! specific humidity " 636 IF( iom_use('theta_zu') ) CALL iom_put("theta_zu", (t_zu -rt0) * tmask(:,:,1)) ! potential temperature at z=zu 637 IF( iom_use('q_zu') ) CALL iom_put("q_zu", q_zu * tmask(:,:,1)) ! specific humidity " 638 IF( iom_use('ssq') ) CALL iom_put("ssq", pssq * tmask(:,:,1)) ! saturation specific humidity at z=0 639 IF( iom_use('wspd_blk') ) CALL iom_put("wspd_blk", zU_zu * tmask(:,:,1)) ! bulk wind speed at z=zu 640 630 641 IF( ln_skin_cs .OR. ln_skin_wl ) THEN 631 642 !! ptsk and pssq have been updated!!! … … 643 654 644 655 IF( ln_abl ) THEN !== ABL formulation ==! multiplication by rho_air and turbulent fluxes computation done in ablstp 645 !! FL do we need this multiplication by tmask ... ???646 656 DO_2D_11_11 647 zztmp = zU_zu(ji,jj) !* tmask(ji,jj,1)657 zztmp = zU_zu(ji,jj) 648 658 wndm(ji,jj) = zztmp ! Store zU_zu in wndm to compute ustar2 in ablmod 649 659 pcd_du(ji,jj) = zztmp * zcd_oce(ji,jj) 650 660 psen(ji,jj) = zztmp * zch_oce(ji,jj) 651 661 pevp(ji,jj) = zztmp * zce_oce(ji,jj) 662 rhoa(ji,jj) = rho_air( ptair(ji,jj), phumi(ji,jj), pslp(ji,jj) ) 652 663 END_2D 653 664 ELSE !== BLK formulation ==! turbulent fluxes computation … … 673 684 ! ... utau, vtau at U- and V_points, resp. 674 685 ! Note the use of 0.5*(2-umask) in order to unmask the stress along coastlines 675 ! Note th e use of MAX(tmask(i,j),tmask(i+1,j) is to mask tau over ice shelves676 DO_2D_ 10_10686 ! Note that coastal wind stress is not used in the code... so this extra care has no effect 687 DO_2D_00_00 677 688 utau(ji,jj) = 0.5 * ( 2. - umask(ji,jj,1) ) * ( zwnd_i(ji,jj) + zwnd_i(ji+1,jj ) ) & 678 689 & * MAX(tmask(ji,jj,1),tmask(ji+1,jj,1)) … … 878 889 Ce_ice(:,:) = Ch_ice(:,:) ! sensible and latent heat transfer coef. are considered identical 879 890 ENDIF 880 881 !! IF ( iom_use("Cd_ice") ) CALL iom_put("Cd_ice", Cd_ice) ! output value of pure ice-atm. transfer coef. 882 !! IF ( iom_use("Ch_ice") ) CALL iom_put("Ch_ice", Ch_ice) ! output value of pure ice-atm. transfer coef. 883 891 892 IF( iom_use('Cd_ice') ) CALL iom_put("Cd_ice", Cd_ice) 893 IF( iom_use('Ce_ice') ) CALL iom_put("Ce_ice", Ce_ice) 894 IF( iom_use('Ch_ice') ) CALL iom_put("Ch_ice", Ch_ice) 895 884 896 ! local scalars ( place there for vector optimisation purposes) 885 !IF (ln_abl) rhoa (:,:) = rho_air( ptair(:,:), phumi(:,:), pslp(:,:) ) !!GS: rhoa must be (re)computed here with ABL to avoid division by zero after (TBI)886 897 zcd_dui(:,:) = wndm_ice(:,:) * Cd_ice(:,:) 887 898 888 899 IF( ln_blk ) THEN 889 ! ------------------------------------------------------------ ! 890 ! Wind stress relative to the moving ice ( U10m - U_ice ) ! 891 ! ------------------------------------------------------------ ! 892 ! C-grid ice dynamics : U & V-points (same as ocean) 893 DO_2D_00_00 894 putaui(ji,jj) = 0.5_wp * ( rhoa(ji+1,jj) * zcd_dui(ji+1,jj) & 895 & + rhoa(ji ,jj) * zcd_dui(ji ,jj) ) & 896 & * ( 0.5_wp * ( pwndi(ji+1,jj) + pwndi(ji,jj) ) - rn_vfac * puice(ji,jj) ) 897 pvtaui(ji,jj) = 0.5_wp * ( rhoa(ji,jj+1) * zcd_dui(ji,jj+1) & 898 & + rhoa(ji,jj ) * zcd_dui(ji,jj ) ) & 899 & * ( 0.5_wp * ( pwndj(ji,jj+1) + pwndj(ji,jj) ) - rn_vfac * pvice(ji,jj) ) 900 ! ------------------------------------------------------------- ! 901 ! Wind stress relative to the moving ice ( U10m - U_ice ) ! 902 ! ------------------------------------------------------------- ! 903 zztmp1 = rn_vfac * 0.5_wp 904 DO_2D_01_01 ! at T point 905 putaui(ji,jj) = rhoa(ji,jj) * zcd_dui(ji,jj) * ( pwndi(ji,jj) - zztmp1 * ( puice(ji-1,jj ) + puice(ji,jj) ) ) 906 pvtaui(ji,jj) = rhoa(ji,jj) * zcd_dui(ji,jj) * ( pwndj(ji,jj) - zztmp1 * ( pvice(ji ,jj-1) + pvice(ji,jj) ) ) 907 END_2D 908 ! 909 DO_2D_00_00 ! U & V-points (same as ocean). 910 ! take care of the land-sea mask to avoid "pollution" of coastal stress. p[uv]taui used in frazil and rheology 911 zztmp1 = 0.5_wp * ( 2. - umask(ji,jj,1) ) * MAX( tmask(ji,jj,1),tmask(ji+1,jj ,1) ) 912 zztmp2 = 0.5_wp * ( 2. - vmask(ji,jj,1) ) * MAX( tmask(ji,jj,1),tmask(ji ,jj+1,1) ) 913 putaui(ji,jj) = zztmp1 * ( putaui(ji,jj) + putaui(ji+1,jj ) ) 914 pvtaui(ji,jj) = zztmp2 * ( pvtaui(ji,jj) + pvtaui(ji ,jj+1) ) 900 915 END_2D 901 916 CALL lbc_lnk_multi( 'sbcblk', putaui, 'U', -1., pvtaui, 'V', -1. ) -
NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/src/OCE/SBC/sbcblk_phy.F90
r12724 r13193 31 31 REAL(wp), PARAMETER, PUBLIC :: R_vap = 461.495_wp !: Specific gas constant for water vapor [J/K/kg] 32 32 REAL(wp), PARAMETER, PUBLIC :: reps0 = R_dry/R_vap !: ratio of gas constant for dry air and water vapor => ~ 0.622 33 REAL(wp), PARAMETER, PUBLIC :: rctv0 = R_vap/R_dry !: for virtual temperature (== (1-eps)/eps) => ~ 0.60833 REAL(wp), PARAMETER, PUBLIC :: rctv0 = R_vap/R_dry - 1._wp !: for virtual temperature (== (1-eps)/eps) => ~ 0.608 34 34 REAL(wp), PARAMETER, PUBLIC :: rCp_air = 1000.5_wp !: specific heat of air (only used for ice fluxes now...) 35 35 REAL(wp), PARAMETER, PUBLIC :: rCd_ice = 1.4e-3_wp !: transfer coefficient over ice -
NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/src/OCE/SBC/sbccpl.F90
r12724 r13193 365 365 ! 366 366 ! Vectors: change of sign at north fold ONLY if on the local grid 367 IF( TRIM( sn_rcv_tau%cldes ) == 'oce only' .OR. TRIM(sn_rcv_tau%cldes ) == 'oce and ice') THEN ! avoid working with the atmospheric fields if they are not coupled 367 IF( TRIM( sn_rcv_tau%cldes ) == 'oce only' .OR. TRIM( sn_rcv_tau%cldes ) == 'oce and ice' & 368 .OR. TRIM( sn_rcv_tau%cldes ) == 'mixed oce-ice' ) THEN ! avoid working with the atmospheric fields if they are not coupled 369 368 370 IF( TRIM( sn_rcv_tau%clvor ) == 'local grid' ) srcv(jpr_otx1:jpr_itz2)%nsgn = -1. 369 371 … … 1482 1484 INTEGER :: ji, jj ! dummy loop indices 1483 1485 INTEGER :: itx ! index of taux over ice 1486 REAL(wp) :: zztmp1, zztmp2 1484 1487 REAL(wp), DIMENSION(jpi,jpj) :: ztx, zty 1485 1488 !!---------------------------------------------------------------------- … … 1545 1548 p_taui(:,:) = frcv(jpr_itx1)%z3(:,:,1) ! (U,V) ==> (U,V) 1546 1549 p_tauj(:,:) = frcv(jpr_ity1)%z3(:,:,1) 1547 CASE( 'F' )1548 DO_2D_00_001549 p_taui(ji,jj) = 0.5 * ( frcv(jpr_itx1)%z3(ji,jj,1) + frcv(jpr_itx1)%z3(ji ,jj-1,1) )1550 p_tauj(ji,jj) = 0.5 * ( frcv(jpr_ity1)%z3(ji,jj,1) + frcv(jpr_ity1)%z3(ji-1,jj ,1) )1551 END_2D1552 1550 CASE( 'T' ) 1553 1551 DO_2D_00_00 1554 p_taui(ji,jj) = 0.5 * ( frcv(jpr_itx1)%z3(ji+1,jj ,1) + frcv(jpr_itx1)%z3(ji,jj,1) ) 1555 p_tauj(ji,jj) = 0.5 * ( frcv(jpr_ity1)%z3(ji ,jj+1,1) + frcv(jpr_ity1)%z3(ji,jj,1) ) 1552 ! take care of the land-sea mask to avoid "pollution" of coastal stress. p[uv]taui used in frazil and rheology 1553 zztmp1 = 0.5_wp * ( 2. - umask(ji,jj,1) ) * MAX( tmask(ji,jj,1),tmask(ji+1,jj ,1) ) 1554 zztmp2 = 0.5_wp * ( 2. - vmask(ji,jj,1) ) * MAX( tmask(ji,jj,1),tmask(ji ,jj+1,1) ) 1555 p_taui(ji,jj) = zztmp1 * ( frcv(jpr_itx1)%z3(ji+1,jj ,1) + frcv(jpr_itx1)%z3(ji,jj,1) ) 1556 p_tauj(ji,jj) = zztmp2 * ( frcv(jpr_ity1)%z3(ji ,jj+1,1) + frcv(jpr_ity1)%z3(ji,jj,1) ) 1556 1557 END_2D 1557 CASE( 'I' ) 1558 DO_2D_00_00 1559 p_taui(ji,jj) = 0.5 * ( frcv(jpr_itx1)%z3(ji+1,jj+1,1) + frcv(jpr_itx1)%z3(ji+1,jj ,1) ) 1560 p_tauj(ji,jj) = 0.5 * ( frcv(jpr_ity1)%z3(ji+1,jj+1,1) + frcv(jpr_ity1)%z3(ji ,jj+1,1) ) 1561 END_2D 1558 CALL lbc_lnk_multi( 'sbccpl', p_taui, 'U', -1., p_tauj, 'V', -1. ) 1562 1559 END SELECT 1563 IF( srcv(jpr_itx1)%clgrid /= 'U' ) THEN1564 CALL lbc_lnk_multi( 'sbccpl', p_taui, 'U', -1., p_tauj, 'V', -1. )1565 ENDIF1566 1560 1567 1561 ENDIF … … 1792 1786 ENDDO 1793 1787 ELSE 1794 qns_tot(:,:) =qns_tot(:,:) + picefr(:,:) * frcv(jpr_qnsice)%z3(:,:,1)1788 zqns_tot(:,:) = zqns_tot(:,:) + picefr(:,:) * frcv(jpr_qnsice)%z3(:,:,1) 1795 1789 DO jl = 1, jpl 1796 zqns_tot(:,: ) = zqns_tot(:,:) + picefr(:,:) * frcv(jpr_qnsice)%z3(:,:,1)1797 1790 zqns_ice(:,:,jl) = frcv(jpr_qnsice)%z3(:,:,1) 1798 1791 END DO … … 1935 1928 END DO 1936 1929 ELSE 1937 qsr_tot(:,: ) =qsr_tot(:,:) + picefr(:,:) * frcv(jpr_qsrice)%z3(:,:,1)1930 zqsr_tot(:,:) = zqsr_tot(:,:) + picefr(:,:) * frcv(jpr_qsrice)%z3(:,:,1) 1938 1931 DO jl = 1, jpl 1939 zqsr_tot(:,: ) = zqsr_tot(:,:) + picefr(:,:) * frcv(jpr_qsrice)%z3(:,:,1)1940 1932 zqsr_ice(:,:,jl) = frcv(jpr_qsrice)%z3(:,:,1) 1941 1933 END DO -
NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/src/OCE/SBC/sbcmod.F90
r12724 r13193 120 120 ncom_fsbc = nn_fsbc ! make nn_fsbc available for lib_mpp 121 121 #endif 122 ! !* overwrite namelist parameter using CPP key information123 #if defined key_agrif124 IF( Agrif_Root() ) THEN ! AGRIF zoom (cf r1242: possibility to run without ice in fine grid)125 IF( lk_si3 ) nn_ice = 2126 IF( lk_cice ) nn_ice = 3127 ENDIF128 !!GS: TBD129 !#else130 ! IF( lk_si3 ) nn_ice = 2131 ! IF( lk_cice ) nn_ice = 3132 #endif133 122 ! 134 123 IF(lwp) THEN !* Control print -
NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/src/OCE/SBC/sbcwave.F90
r12622 r13193 212 212 END_3D 213 213 ! 214 #if defined key_agrif215 IF( .NOT. Agrif_Root() ) THEN216 IF( nbondi == -1 .OR. nbondi == 2 ) ze3divh( 2:nbghostcells+1,: ,:) = 0._wp ! west217 IF( nbondi == 1 .OR. nbondi == 2 ) ze3divh( nlci-nbghostcells:nlci-1,:,:) = 0._wp ! east218 IF( nbondj == -1 .OR. nbondj == 2 ) ze3divh( :,2:nbghostcells+1 ,:) = 0._wp ! south219 IF( nbondj == 1 .OR. nbondj == 2 ) ze3divh( :,nlcj-nbghostcells:nlcj-1,:) = 0._wp ! north220 ENDIF221 #endif222 !223 214 CALL lbc_lnk( 'sbcwave', ze3divh, 'T', 1. ) 224 215 ! -
NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/src/OCE/STO/stopar.F90
r12377 r13193 684 684 !! ** Purpose : read stochastic parameters from restart file 685 685 !!---------------------------------------------------------------------- 686 INTEGER :: jsto, jseed 686 INTEGER :: jsto, jseed 687 INTEGER :: idg ! number of digits 687 688 INTEGER(KIND=8) :: ziseed(4) ! RNG seeds in integer type 688 689 REAL(KIND=8) :: zrseed(4) ! RNG seeds in real type (with same bits to save in restart) 689 690 CHARACTER(LEN=9) :: clsto2d='sto2d_000' ! stochastic parameter variable name 690 691 CHARACTER(LEN=9) :: clsto3d='sto3d_000' ! stochastic parameter variable name 691 CHARACTER(LEN=10) :: clseed='seed0_0000' ! seed variable name 692 CHARACTER(LEN=15) :: clseed='seed0_0000' ! seed variable name 693 CHARACTER(LEN=6) :: clfmt ! writing format 692 694 !!---------------------------------------------------------------------- 693 695 … … 717 719 IF (ln_rstseed) THEN 718 720 ! Get saved state of the random number generator 721 idg = MAX( INT(LOG10(REAL(jpnij,wp))) + 1, 4 ) ! how many digits to we need to write? min=4, max=9 722 WRITE(clfmt, "('(i', i1, '.', i1, ')')") idg, idg ! "(ix.x)" 719 723 DO jseed = 1 , 4 720 WRITE(clseed(5:5) ,'(i1.1)') jseed721 WRITE(clseed(7: 10),'(i4.4)') narea722 CALL iom_get( numstor, clseed , zrseed(jseed) )724 WRITE(clseed(5:5) ,'(i1.1)') jseed 725 WRITE(clseed(7:7+idg-1), clfmt ) narea 726 CALL iom_get( numstor, clseed(1:7+idg-1) , zrseed(jseed) ) 723 727 END DO 724 728 ziseed = TRANSFER( zrseed , ziseed) … … 742 746 INTEGER, INTENT(in) :: kt ! ocean time-step 743 747 !! 744 INTEGER :: jsto, jseed 748 INTEGER :: jsto, jseed 749 INTEGER :: idg ! number of digits 745 750 INTEGER(KIND=8) :: ziseed(4) ! RNG seeds in integer type 746 751 REAL(KIND=8) :: zrseed(4) ! RNG seeds in real type (with same bits to save in restart) … … 749 754 CHARACTER(LEN=9) :: clsto2d='sto2d_000' ! stochastic parameter variable name 750 755 CHARACTER(LEN=9) :: clsto3d='sto3d_000' ! stochastic parameter variable name 751 CHARACTER(LEN=10) :: clseed='seed0_0000' ! seed variable name 756 CHARACTER(LEN=15) :: clseed='seed0_0000' ! seed variable name 757 CHARACTER(LEN=6) :: clfmt ! writing format 752 758 !!---------------------------------------------------------------------- 753 759 … … 771 777 CALL kiss_state( ziseed(1) , ziseed(2) , ziseed(3) , ziseed(4) ) 772 778 zrseed = TRANSFER( ziseed , zrseed) 779 idg = MAX( INT(LOG10(REAL(jpnij,wp))) + 1, 4 ) ! how many digits to we need to write? min=4, max=9 780 WRITE(clfmt, "('(i', i1, '.', i1, ')')") idg, idg ! "(ix.x)" 773 781 DO jseed = 1 , 4 774 WRITE(clseed(5:5) ,'(i1.1)') jseed775 WRITE(clseed(7: 10),'(i4.4)') narea776 CALL iom_rstput( kt, nitrst, numstow, clseed 782 WRITE(clseed(5:5) ,'(i1.1)') jseed 783 WRITE(clseed(7:7+idg-1), clfmt ) narea 784 CALL iom_rstput( kt, nitrst, numstow, clseed(1:7+idg-1), zrseed(jseed) ) 777 785 END DO 778 786 ! 2D stochastic parameters -
NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/src/OCE/TRD/trdtra.F90
r12724 r13193 83 83 REAL(wp), DIMENSION(jpi,jpj,jpk), INTENT(in), OPTIONAL :: ptra ! now tracer variable 84 84 ! 85 INTEGER :: jk ! loop indices 85 INTEGER :: jk ! loop indices 86 INTEGER :: i01 ! 0 or 1 86 87 REAL(wp), DIMENSION(jpi,jpj,jpk) :: ztrds ! 3D workspace 87 88 REAL(wp), ALLOCATABLE, DIMENSION(:,:,:) :: zwt, zws, ztrdt ! 3D workspace … … 91 92 IF( trd_tra_alloc() /= 0 ) CALL ctl_stop( 'STOP', 'trd_tra : unable to allocate arrays' ) 92 93 ENDIF 93 94 ! 95 i01 = COUNT( (/ PRESENT(pu) .OR. ( ktrd /= jptra_xad .AND. ktrd /= jptra_yad .AND. ktrd /= jptra_zad ) /) ) 96 ! 94 97 IF( ctype == 'TRA' .AND. ktra == jp_tem ) THEN !== Temperature trend ==! 95 98 ! 96 SELECT CASE( ktrd )99 SELECT CASE( ktrd*i01 ) 97 100 ! ! advection: transform the advective flux into a trend 98 101 CASE( jptra_xad ) ; CALL trd_tra_adv( ptrd, pu, ptra, 'X', trdtx, Kmm ) … … 113 116 IF( ctype == 'TRA' .AND. ktra == jp_sal ) THEN !== Salinity trends ==! 114 117 ! 115 SELECT CASE( ktrd )118 SELECT CASE( ktrd*i01 ) 116 119 ! ! advection: transform the advective flux into a trend 117 120 ! ! and send T & S trends to trd_tra_mng … … 168 171 IF( ctype == 'TRC' ) THEN !== passive tracer trend ==! 169 172 ! 170 SELECT CASE( ktrd )173 SELECT CASE( ktrd*i01 ) 171 174 ! ! advection: transform the advective flux into a masked trend 172 175 CASE( jptra_xad ) ; CALL trd_tra_adv( ptrd , pu , ptra, 'X', ztrds, Kmm ) -
NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/src/OCE/USR/usrdef_zgr.F90
r12377 r13193 202 202 CALL lbc_lnk( 'usrdef_zgr', z2d, 'T', 1. ) ! set surrounding land to zero (here jperio=0 ==>> closed) 203 203 ! 204 k_bot(:,:) = NINT( z2d(:,:) ) 204 k_bot(:,:) = NINT( z2d(:,:) ) ! =jpkm1 over the ocean point, =0 elsewhere 205 205 ! 206 206 k_top(:,:) = MIN( 1 , k_bot(:,:) ) ! = 1 over the ocean point, =0 elsewhere -
NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/src/OCE/ZDF/zdftke.F90
r12724 r13193 45 45 USE zdfdrg ! vertical physics: top/bottom drag coef. 46 46 USE zdfmxl ! vertical physics: mixed layer 47 #if defined key_si3 48 USE ice, ONLY: hm_i, h_i 49 #endif 50 #if defined key_cice 51 USE sbc_ice, ONLY: h_i 52 #endif 47 53 ! 48 54 USE in_out_manager ! I/O manager … … 64 70 INTEGER :: nn_mxl ! type of mixing length (=0/1/2/3) 65 71 REAL(wp) :: rn_mxl0 ! surface min value of mixing length (kappa*z_o=0.4*0.1 m) [m] 72 INTEGER :: nn_mxlice ! type of scaling under sea-ice 73 REAL(wp) :: rn_mxlice ! max constant ice thickness value when scaling under sea-ice ( nn_mxlice=1) 66 74 INTEGER :: nn_pdl ! Prandtl number or not (ratio avt/avm) (=0/1) 67 75 REAL(wp) :: rn_ediff ! coefficient for avt: avt=rn_ediff*mxl*sqrt(e) … … 245 253 zetop = - 0.001875_wp * rCdU_top(ji,jj) * SQRT( ( zmsku*( uu(ji,jj,mikt(ji,jj),Kbb)+uu(ji-1,jj,mikt(ji,jj),Kbb) ) )**2 & 246 254 & + ( zmskv*( vv(ji,jj,mikt(ji,jj),Kbb)+vv(ji,jj-1,mikt(ji,jj),Kbb) ) )**2 ) 247 en(ji,jj,mikt(ji,jj)) = en(ji,jj,1) * tmask(ji,jj,1) + MAX( zetop, rn_emin ) * (1._wp - tmask(ji,jj,1)) * ssmask(ji,jj) ! masked at ocean surface 255 ! (1._wp - tmask(ji,jj,1)) * ssmask(ji,jj) = 1 where ice shelves are present 256 en(ji,jj,mikt(ji,jj)) = en(ji,jj,1) * tmask(ji,jj,1) & 257 & + MAX( zetop, rn_emin ) * (1._wp - tmask(ji,jj,1)) * ssmask(ji,jj) 248 258 END_2D 249 259 ENDIF … … 424 434 REAL(wp) :: zrn2, zraug, zcoef, zav ! local scalars 425 435 REAL(wp) :: zdku, zdkv, zsqen ! - - 426 REAL(wp) :: zemxl, zemlm, zemlp ! - -436 REAL(wp) :: zemxl, zemlm, zemlp, zmaxice ! - - 427 437 REAL(wp), DIMENSION(jpi,jpj,jpk) :: zmxlm, zmxld ! 3D workspace 428 438 !!-------------------------------------------------------------------- … … 438 448 zmxld(:,:,:) = rmxl_min 439 449 ! 440 IF( ln_mxl0 ) THEN ! surface mixing length = F(stress) : l=vkarmn*2.e5*taum/(rho0*g) 450 IF( ln_mxl0 ) THEN ! surface mixing length = F(stress) : l=vkarmn*2.e5*taum/(rho0*g) 451 ! 441 452 zraug = vkarmn * 2.e5_wp / ( rho0 * grav ) 453 #if ! defined key_si3 && ! defined key_cice 442 454 DO_2D_00_00 443 zmxlm(ji,jj,1) = MAX( rn_mxl0, zraug * taum(ji,jj) * tmask(ji,jj,1))455 zmxlm(ji,jj,1) = zraug * taum(ji,jj) * tmask(ji,jj,1) 444 456 END_2D 445 ELSE 457 #else 458 SELECT CASE( nn_mxlice ) ! Type of scaling under sea-ice 459 ! 460 CASE( 0 ) ! No scaling under sea-ice 461 DO_2D_00_00 462 zmxlm(ji,jj,1) = zraug * taum(ji,jj) * tmask(ji,jj,1) 463 END_2D 464 ! 465 CASE( 1 ) ! scaling with constant sea-ice thickness 466 DO_2D_00_00 467 zmxlm(ji,jj,1) = ( ( 1. - fr_i(ji,jj) ) * zraug * taum(ji,jj) + fr_i(ji,jj) * rn_mxlice ) * tmask(ji,jj,1) 468 END_2D 469 ! 470 CASE( 2 ) ! scaling with mean sea-ice thickness 471 DO_2D_00_00 472 #if defined key_si3 473 zmxlm(ji,jj,1) = ( ( 1. - fr_i(ji,jj) ) * zraug * taum(ji,jj) + fr_i(ji,jj) * hm_i(ji,jj) * 2. ) * tmask(ji,jj,1) 474 #elif defined key_cice 475 zmaxice = MAXVAL( h_i(ji,jj,:) ) 476 zmxlm(ji,jj,1) = ( ( 1. - fr_i(ji,jj) ) * zraug * taum(ji,jj) + fr_i(ji,jj) * zmaxice ) * tmask(ji,jj,1) 477 #endif 478 END_2D 479 ! 480 CASE( 3 ) ! scaling with max sea-ice thickness 481 DO_2D_00_00 482 zmaxice = MAXVAL( h_i(ji,jj,:) ) 483 zmxlm(ji,jj,1) = ( ( 1. - fr_i(ji,jj) ) * zraug * taum(ji,jj) + fr_i(ji,jj) * zmaxice ) * tmask(ji,jj,1) 484 END_2D 485 ! 486 END SELECT 487 #endif 488 ! 489 DO_2D_00_00 490 zmxlm(ji,jj,1) = MAX( rn_mxl0, zmxlm(ji,jj,1) ) 491 END_2D 492 ! 493 ELSE 446 494 zmxlm(:,:,1) = rn_mxl0 447 495 ENDIF 496 448 497 ! 449 498 DO_3D_00_00( 2, jpkm1 ) … … 554 603 INTEGER :: ios 555 604 !! 556 NAMELIST/namzdf_tke/ rn_ediff, rn_ediss , rn_ebb , rn_emin , & 557 & rn_emin0, rn_bshear, nn_mxl , ln_mxl0 , & 558 & rn_mxl0 , nn_pdl , ln_drg , ln_lc , rn_lc, & 559 & nn_etau , nn_htau , rn_efr , rn_eice 605 NAMELIST/namzdf_tke/ rn_ediff, rn_ediss , rn_ebb , rn_emin , & 606 & rn_emin0, rn_bshear, nn_mxl , ln_mxl0 , & 607 & rn_mxl0 , nn_mxlice, rn_mxlice, & 608 & nn_pdl , ln_drg , ln_lc , rn_lc, & 609 & nn_etau , nn_htau , rn_efr , rn_eice 560 610 !!---------------------------------------------------------------------- 561 611 ! … … 583 633 WRITE(numout,*) ' mixing length type nn_mxl = ', nn_mxl 584 634 WRITE(numout,*) ' surface mixing length = F(stress) or not ln_mxl0 = ', ln_mxl0 635 IF( ln_mxl0 ) THEN 636 WRITE(numout,*) ' type of scaling under sea-ice nn_mxlice = ', nn_mxlice 637 IF( nn_mxlice == 1 ) & 638 WRITE(numout,*) ' ice thickness when scaling under sea-ice rn_mxlice = ', rn_mxlice 639 ENDIF 585 640 WRITE(numout,*) ' surface mixing length minimum value rn_mxl0 = ', rn_mxl0 586 641 WRITE(numout,*) ' top/bottom friction forcing flag ln_drg = ', ln_drg -
NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/src/OCE/nemogcm.F90
r12732 r13193 194 194 END DO 195 195 ! 196 IF( .NOT. Agrif_Root() ) THEN197 CALL Agrif_ParentGrid_To_ChildGrid()198 IF( ln_diaobs ) CALL dia_obs_wri199 IF( ln_timing ) CALL timing_finalize200 CALL Agrif_ChildGrid_To_ParentGrid()201 ENDIF202 !203 196 # else 204 197 ! … … 249 242 IF( nstop /= 0 .AND. lwp ) THEN ! error print 250 243 WRITE(ctmp1,*) ' ==>>> nemo_gcm: a total of ', nstop, ' errors have been found' 251 CALL ctl_stop( ctmp1 ) 244 IF( ngrdstop > 0 ) THEN 245 WRITE(ctmp9,'(i2)') ngrdstop 246 WRITE(ctmp2,*) ' E R R O R detected in Agrif grid '//TRIM(ctmp9) 247 WRITE(ctmp3,*) ' Look for "E R R O R" messages in all existing '//TRIM(ctmp9)//'_ocean_output* files' 248 CALL ctl_stop( ' ', ctmp1, ' ', ctmp2, ' ', ctmp3 ) 249 ELSE 250 WRITE(ctmp2,*) ' Look for "E R R O R" messages in all existing ocean_output* files' 251 CALL ctl_stop( ' ', ctmp1, ' ', ctmp2 ) 252 ENDIF 252 253 ENDIF 253 254 ! … … 261 262 #else 262 263 IF ( lk_oasis ) THEN ; CALL cpl_finalize ! end coupling and mpp communications with OASIS 263 ELSEIF( lk_mpp ) THEN ; CALL mppstop ! end mpp communications264 ELSEIF( lk_mpp ) THEN ; CALL mppstop ! end mpp communications 264 265 ENDIF 265 266 #endif … … 347 348 ! 348 349 ! finalize the definition of namctl variables 349 IF( sn_cfctl%l_allon ) THEN 350 ! Turn on all options. 351 CALL nemo_set_cfctl( sn_cfctl, .TRUE., .TRUE. ) 352 ! Ensure all processors are active 353 sn_cfctl%procmin = 0 ; sn_cfctl%procmax = 1000000 ; sn_cfctl%procincr = 1 354 ELSEIF( sn_cfctl%l_config ) THEN 355 ! Activate finer control of report outputs 356 ! optionally switch off output from selected areas (note this only 357 ! applies to output which does not involve global communications) 358 IF( ( narea < sn_cfctl%procmin .OR. narea > sn_cfctl%procmax ) .OR. & 359 & ( MOD( narea - sn_cfctl%procmin, sn_cfctl%procincr ) /= 0 ) ) & 360 & CALL nemo_set_cfctl( sn_cfctl, .FALSE., .FALSE. ) 361 ELSE 362 ! turn off all options. 363 CALL nemo_set_cfctl( sn_cfctl, .FALSE., .TRUE. ) 364 ENDIF 350 IF( narea < sn_cfctl%procmin .OR. narea > sn_cfctl%procmax .OR. MOD( narea - sn_cfctl%procmin, sn_cfctl%procincr ) /= 0 ) & 351 & CALL nemo_set_cfctl( sn_cfctl, .FALSE. ) 365 352 ! 366 353 lwp = (narea == 1) .OR. sn_cfctl%l_oceout ! control of all listing output print … … 546 533 WRITE(numout,*) '~~~~~~~~' 547 534 WRITE(numout,*) ' Namelist namctl' 548 WRITE(numout,*) ' sn_cfctl%l_glochk = ', sn_cfctl%l_glochk549 WRITE(numout,*) ' sn_cfctl%l_allon = ', sn_cfctl%l_allon550 WRITE(numout,*) ' finer control over o/p sn_cfctl%l_config = ', sn_cfctl%l_config551 535 WRITE(numout,*) ' sn_cfctl%l_runstat = ', sn_cfctl%l_runstat 552 536 WRITE(numout,*) ' sn_cfctl%l_trcstat = ', sn_cfctl%l_trcstat … … 696 680 697 681 698 SUBROUTINE nemo_set_cfctl(sn_cfctl, setto , for_all)682 SUBROUTINE nemo_set_cfctl(sn_cfctl, setto ) 699 683 !!---------------------------------------------------------------------- 700 684 !! *** ROUTINE nemo_set_cfctl *** 701 685 !! 702 686 !! ** Purpose : Set elements of the output control structure to setto. 703 !! for_all should be .false. unless all areas are to be704 !! treated identically.705 687 !! 706 688 !! ** Method : Note this routine can be used to switch on/off some 707 !! types of output for selected areas but any output types 708 !! that involve global communications (e.g. mpp_max, glob_sum) 709 !! should be protected from selective switching by the 710 !! for_all argument 711 !!---------------------------------------------------------------------- 712 LOGICAL :: setto, for_all 713 TYPE(sn_ctl) :: sn_cfctl 714 !!---------------------------------------------------------------------- 715 IF( for_all ) THEN 716 sn_cfctl%l_runstat = setto 717 sn_cfctl%l_trcstat = setto 718 ENDIF 689 !! types of output for selected areas. 690 !!---------------------------------------------------------------------- 691 TYPE(sn_ctl), INTENT(inout) :: sn_cfctl 692 LOGICAL , INTENT(in ) :: setto 693 !!---------------------------------------------------------------------- 694 sn_cfctl%l_runstat = setto 695 sn_cfctl%l_trcstat = setto 719 696 sn_cfctl%l_oceout = setto 720 697 sn_cfctl%l_layout = setto -
NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/src/OCE/step.F90
r12724 r13193 86 86 !!---------------------------------------------------------------------- 87 87 INTEGER :: ji, jj, jk ! dummy loop indice 88 INTEGER :: indic ! error indicator if < 089 88 !!gm kcall can be removed, I guess 90 89 INTEGER :: kcall ! optional integer argument (dom_vvl_sf_nxt) 91 90 !! --------------------------------------------------------------------- 92 91 #if defined key_agrif 93 IF( nstop > 0 ) return ! avoid to go further if an error was detected during previous time step92 IF( nstop > 0 ) RETURN ! avoid to go further if an error was detected during previous time step (child grid) 94 93 kstp = nit000 + Agrif_Nb_Step() 95 94 Kbb_a = Nbb; Kmm_a = Nnn; Krhs_a = Nrhs ! agrif_oce module copies of time level indices … … 119 118 ! update I/O and calendar 120 119 !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 121 indic = 0 ! reset to no error condition122 123 120 IF( kstp == nit000 ) THEN ! initialize IOM context (must be done after nemo_init for AGRIF+XIOS+OASIS) 124 CALL iom_init( cxios_context, ld_closedef=.FALSE. ) ! for model grid (including p assible AGRIF zoom)121 CALL iom_init( cxios_context, ld_closedef=.FALSE. ) ! for model grid (including possible AGRIF zoom) 125 122 IF( lk_diamlr ) CALL dia_mlr_iom_init ! with additional setup for multiple-linear-regression analysis 126 123 CALL iom_init_closedef … … 318 315 Kbb_a = Nbb; Kmm_a = Nnn; Krhs_a = Nrhs ! agrif_oce module copies of time level indices 319 316 CALL Agrif_Integrate_ChildGrids( stp ) ! allows to finish all the Child Grids before updating 317 320 318 #endif 321 319 !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 322 320 ! Control 323 321 !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 324 CALL stp_ctl ( kstp, Nbb, Nnn, indic ) 322 CALL stp_ctl ( kstp, Nnn ) 323 325 324 #if defined key_agrif 326 325 !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> … … 330 329 CALL Agrif_update_all( ) ! Update all components 331 330 ENDIF 332 #endif 333 IF( ln_diaobs ) CALL dia_obs ( kstp, Nnn ) ! obs-minus-model (assimilation) diagnostics (call after dynamics update) 331 332 #endif 333 IF( ln_diaobs .AND. nstop == 0 ) CALL dia_obs( kstp, Nnn ) ! obs-minus-model (assimilation) diags (after dynamics update) 334 334 335 335 !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> … … 345 345 ! Coupled mode 346 346 !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 347 !!gm why lk_oasis and not lk_cpl ???? 348 IF( lk_oasis ) CALL sbc_cpl_snd( kstp, Nbb, Nnn ) ! coupled mode : field exchanges 347 IF( lk_oasis .AND. nstop == 0 ) CALL sbc_cpl_snd( kstp, Nbb, Nnn ) ! coupled mode : field exchanges 349 348 ! 350 349 #if defined key_iomput … … 352 351 ! Finalize contextes if end of simulation or error detected 353 352 !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 354 IF( kstp == nitend .OR. indic <0 ) THEN353 IF( kstp == nitend .OR. nstop > 0 ) THEN 355 354 CALL iom_context_finalize( cxios_context ) ! needed for XIOS+AGRIF 356 355 IF( lrxios ) CALL iom_context_finalize( crxios_context ) -
NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/src/OCE/stpctl.F90
r12377 r13193 19 19 USE dom_oce ! ocean space and time domain variables 20 20 USE c1d ! 1D vertical configuration 21 USE zdf_oce , ONLY : ln_zad_Aimp ! ocean vertical physics variables 22 USE wet_dry, ONLY : ll_wd, ssh_ref ! reference depth for negative bathy 23 ! 21 24 USE diawri ! Standard run outputs (dia_wri_state routine) 22 !23 25 USE in_out_manager ! I/O manager 24 26 USE lbclnk ! ocean lateral boundary conditions (or mpp link) 25 27 USE lib_mpp ! distributed memory computing 26 USE zdf_oce , ONLY : ln_zad_Aimp ! ocean vertical physics variables 27 USE wet_dry, ONLY : ll_wd, ssh_ref ! reference depth for negative bathy 28 28 ! 29 29 USE netcdf ! NetCDF library 30 30 IMPLICIT NONE … … 33 33 PUBLIC stp_ctl ! routine called by step.F90 34 34 35 INTEGER :: idrun, idtime, idssh, idu, ids1, ids2, idt1, idt2, idc1, idw1, istatus36 LOGICAL :: lsomeoce35 INTEGER :: nrunid ! netcdf file id 36 INTEGER, DIMENSION(8) :: nvarid ! netcdf variable id 37 37 !!---------------------------------------------------------------------- 38 38 !! NEMO/OCE 4.0 , NEMO Consortium (2018) … … 42 42 CONTAINS 43 43 44 SUBROUTINE stp_ctl( kt, K bb, Kmm, kindic)44 SUBROUTINE stp_ctl( kt, Kmm ) 45 45 !!---------------------------------------------------------------------- 46 46 !! *** ROUTINE stp_ctl *** … … 50 50 !! ** Method : - Save the time step in numstp 51 51 !! - Print it each 50 time steps 52 !! - Stop the run IF problem encountered by setting indic=-352 !! - Stop the run IF problem encountered by setting nstop > 0 53 53 !! Problems checked: |ssh| maximum larger than 10 m 54 54 !! |U| maximum larger than 10 m/s … … 57 57 !! ** Actions : "time.step" file = last ocean time-step 58 58 !! "run.stat" file = run statistics 59 !! nstop indicator sheared among all local domain (lk_mpp=T)59 !! nstop indicator sheared among all local domain 60 60 !!---------------------------------------------------------------------- 61 61 INTEGER, INTENT(in ) :: kt ! ocean time-step index 62 INTEGER, INTENT(in ) :: Kbb, Kmm ! ocean time level index 63 INTEGER, INTENT(inout) :: kindic ! error indicator 64 !! 65 INTEGER :: ji, jj, jk ! dummy loop indices 66 INTEGER, DIMENSION(2) :: ih ! min/max loc indices 67 INTEGER, DIMENSION(3) :: iu, is1, is2 ! min/max loc indices 68 REAL(wp) :: zzz ! local real 69 REAL(wp), DIMENSION(9) :: zmax 70 LOGICAL :: ll_wrtstp, ll_colruns, ll_wrtruns 71 CHARACTER(len=20) :: clname 72 !!---------------------------------------------------------------------- 73 ! 74 ll_wrtstp = ( MOD( kt, sn_cfctl%ptimincr ) == 0 ) .OR. ( kt == nitend ) 75 ll_colruns = ll_wrtstp .AND. ( sn_cfctl%l_runstat ) 76 ll_wrtruns = ll_colruns .AND. lwm 77 IF( kt == nit000 .AND. lwp ) THEN 78 WRITE(numout,*) 79 WRITE(numout,*) 'stp_ctl : time-stepping control' 80 WRITE(numout,*) '~~~~~~~' 81 ! ! open time.step file 82 IF( lwm ) CALL ctl_opn( numstp, 'time.step', 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, numout, lwp, narea ) 83 ! ! open run.stat file(s) at start whatever 84 ! ! the value of sn_cfctl%ptimincr 85 IF( lwm .AND. ( sn_cfctl%l_runstat ) ) THEN 62 INTEGER, INTENT(in ) :: Kmm ! ocean time level index 63 !! 64 INTEGER :: ji ! dummy loop indices 65 INTEGER :: idtime, istatus 66 INTEGER , DIMENSION(9) :: iareasum, iareamin, iareamax 67 INTEGER , DIMENSION(3,4) :: iloc ! min/max loc indices 68 REAL(wp) :: zzz ! local real 69 REAL(wp), DIMENSION(9) :: zmax, zmaxlocal 70 LOGICAL :: ll_wrtstp, ll_colruns, ll_wrtruns 71 LOGICAL, DIMENSION(jpi,jpj,jpk) :: llmsk 72 CHARACTER(len=20) :: clname 73 !!---------------------------------------------------------------------- 74 IF( nstop > 0 .AND. ngrdstop > -1 ) RETURN ! stpctl was already called by a child grid 75 ! 76 ll_wrtstp = ( MOD( kt-nit000, sn_cfctl%ptimincr ) == 0 ) .OR. ( kt == nitend ) 77 ll_colruns = ll_wrtstp .AND. sn_cfctl%l_runstat .AND. jpnij > 1 78 ll_wrtruns = ( ll_colruns .OR. jpnij == 1 ) .AND. lwm 79 ! 80 IF( kt == nit000 ) THEN 81 ! 82 IF( lwp ) THEN 83 WRITE(numout,*) 84 WRITE(numout,*) 'stp_ctl : time-stepping control' 85 WRITE(numout,*) '~~~~~~~' 86 ENDIF 87 ! ! open time.step ascii file, done only by 1st subdomain 88 IF( lwm ) CALL ctl_opn( numstp, 'time.step', 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, numout, lwp, narea ) 89 ! 90 IF( ll_wrtruns ) THEN 91 ! ! open run.stat ascii file, done only by 1st subdomain 86 92 CALL ctl_opn( numrun, 'run.stat', 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, numout, lwp, narea ) 93 ! ! open run.stat.nc netcdf file, done only by 1st subdomain 87 94 clname = 'run.stat.nc' 88 95 IF( .NOT. Agrif_Root() ) clname = TRIM(Agrif_CFixed())//"_"//TRIM(clname) 89 istatus = NF90_CREATE( TRIM(clname), NF90_CLOBBER, idrun)90 istatus = NF90_DEF_DIM( idrun, 'time', NF90_UNLIMITED, idtime )91 istatus = NF90_DEF_VAR( idrun, 'abs_ssh_max', NF90_DOUBLE, (/ idtime /), idssh)92 istatus = NF90_DEF_VAR( idrun, 'abs_u_max', NF90_DOUBLE, (/ idtime /), idu)93 istatus = NF90_DEF_VAR( idrun, 's_min', NF90_DOUBLE, (/ idtime /), ids1)94 istatus = NF90_DEF_VAR( idrun, 's_max', NF90_DOUBLE, (/ idtime /), ids2)95 istatus = NF90_DEF_VAR( idrun, 't_min', NF90_DOUBLE, (/ idtime /), idt1)96 istatus = NF90_DEF_VAR( idrun, 't_max', NF90_DOUBLE, (/ idtime /), idt2)96 istatus = NF90_CREATE( TRIM(clname), NF90_CLOBBER, nrunid ) 97 istatus = NF90_DEF_DIM( nrunid, 'time', NF90_UNLIMITED, idtime ) 98 istatus = NF90_DEF_VAR( nrunid, 'abs_ssh_max', NF90_DOUBLE, (/ idtime /), nvarid(1) ) 99 istatus = NF90_DEF_VAR( nrunid, 'abs_u_max', NF90_DOUBLE, (/ idtime /), nvarid(2) ) 100 istatus = NF90_DEF_VAR( nrunid, 's_min', NF90_DOUBLE, (/ idtime /), nvarid(3) ) 101 istatus = NF90_DEF_VAR( nrunid, 's_max', NF90_DOUBLE, (/ idtime /), nvarid(4) ) 102 istatus = NF90_DEF_VAR( nrunid, 't_min', NF90_DOUBLE, (/ idtime /), nvarid(5) ) 103 istatus = NF90_DEF_VAR( nrunid, 't_max', NF90_DOUBLE, (/ idtime /), nvarid(6) ) 97 104 IF( ln_zad_Aimp ) THEN 98 istatus = NF90_DEF_VAR( idrun, 'abs_wi_max', NF90_DOUBLE, (/ idtime /), idw1)99 istatus = NF90_DEF_VAR( idrun, 'Cf_max', NF90_DOUBLE, (/ idtime /), idc1)105 istatus = NF90_DEF_VAR( nrunid, 'Cf_max', NF90_DOUBLE, (/ idtime /), nvarid(7) ) 106 istatus = NF90_DEF_VAR( nrunid,'abs_wi_max',NF90_DOUBLE, (/ idtime /), nvarid(8) ) 100 107 ENDIF 101 istatus = NF90_ENDDEF(idrun) 102 zmax(8:9) = 0._wp ! initialise to zero in case ln_zad_Aimp option is not in use 103 ENDIF 104 ENDIF 105 IF( kt == nit000 ) lsomeoce = COUNT( ssmask(:,:) == 1._wp ) > 0 106 ! 107 IF(lwm .AND. ll_wrtstp) THEN !== current time step ==! ("time.step" file) 108 istatus = NF90_ENDDEF(nrunid) 109 ENDIF 110 ! 111 ENDIF 112 ! 113 ! !== write current time step ==! 114 ! !== done only by 1st subdomain at writting timestep ==! 115 IF( lwm .AND. ll_wrtstp ) THEN 108 116 WRITE ( numstp, '(1x, i8)' ) kt 109 117 REWIND( numstp ) 110 118 ENDIF 111 ! 112 ! !== test of extrema ==! 113 IF( ll_wd ) THEN 114 zmax(1) = MAXVAL( ABS( ssh(:,:,Kmm) + ssh_ref*tmask(:,:,1) ) ) ! ssh max 115 ELSE 116 zmax(1) = MAXVAL( ABS( ssh(:,:,Kmm) ) ) ! ssh max 117 ENDIF 118 zmax(2) = MAXVAL( ABS( uu(:,:,:,Kmm) ) ) ! velocity max (zonal only) 119 zmax(3) = MAXVAL( -ts(:,:,:,jp_sal,Kmm) , mask = tmask(:,:,:) == 1._wp ) ! minus salinity max 120 zmax(4) = MAXVAL( ts(:,:,:,jp_sal,Kmm) , mask = tmask(:,:,:) == 1._wp ) ! salinity max 121 zmax(5) = MAXVAL( -ts(:,:,:,jp_tem,Kmm) , mask = tmask(:,:,:) == 1._wp ) ! minus temperature max 122 zmax(6) = MAXVAL( ts(:,:,:,jp_tem,Kmm) , mask = tmask(:,:,:) == 1._wp ) ! temperature max 123 zmax(7) = REAL( nstop , wp ) ! stop indicator 124 IF( ln_zad_Aimp ) THEN 125 zmax(8) = MAXVAL( ABS( wi(:,:,:) ) , mask = wmask(:,:,:) == 1._wp ) ! implicit vertical vel. max 126 zmax(9) = MAXVAL( Cu_adv(:,:,:) , mask = tmask(:,:,:) == 1._wp ) ! partitioning coeff. max 127 ENDIF 128 ! 119 ! !== test of local extrema ==! 120 ! !== done by all processes at every time step ==! 121 ! 122 ! define zmax default value. needed for land processors 123 IF( ll_colruns ) THEN ! default value: must not be kept when calling mpp_max -> must be as small as possible 124 zmax(:) = -HUGE(1._wp) 125 ELSE ! default value: must not give true for any of the tests bellow (-> avoid manipulating HUGE...) 126 zmax(:) = 0._wp 127 zmax(3) = -1._wp ! avoid salinity minimum at 0. 128 ENDIF 129 ! 130 llmsk(:,:,1) = ssmask(:,:) == 1._wp 131 IF( COUNT( llmsk(:,:,1) ) > 0 ) THEN ! avoid huge values sent back for land processors... 132 IF( ll_wd ) THEN 133 zmax(1) = MAXVAL( ABS( ssh(:,:,Kmm) + ssh_ref ), mask = llmsk(:,:,1) ) ! ssh max 134 ELSE 135 zmax(1) = MAXVAL( ABS( ssh(:,:,Kmm) ), mask = llmsk(:,:,1) ) ! ssh max 136 ENDIF 137 ENDIF 138 zmax(2) = MAXVAL( ABS( uu(:,:,:,Kmm) ) ) ! velocity max (zonal only) 139 llmsk(:,:,:) = tmask(:,:,:) == 1._wp 140 IF( COUNT( llmsk(:,:,:) ) > 0 ) THEN ! avoid huge values sent back for land processors... 141 zmax(3) = MAXVAL( -ts(:,:,:,jp_sal,Kmm), mask = llmsk ) ! minus salinity max 142 zmax(4) = MAXVAL( ts(:,:,:,jp_sal,Kmm), mask = llmsk ) ! salinity max 143 IF( ll_colruns .OR. jpnij == 1 ) THEN ! following variables are used only in the netcdf file 144 zmax(5) = MAXVAL( -ts(:,:,:,jp_tem,Kmm), mask = llmsk ) ! minus temperature max 145 zmax(6) = MAXVAL( ts(:,:,:,jp_tem,Kmm), mask = llmsk ) ! temperature max 146 IF( ln_zad_Aimp ) THEN 147 zmax(7) = MAXVAL( Cu_adv(:,:,:) , mask = llmsk ) ! partitioning coeff. max 148 llmsk(:,:,:) = wmask(:,:,:) == 1._wp 149 IF( COUNT( llmsk(:,:,:) ) > 0 ) THEN ! avoid huge values sent back for land processors... 150 zmax(8) = MAXVAL(ABS( wi(:,:,:) ), mask = llmsk ) ! implicit vertical vel. max 151 ENDIF 152 ENDIF 153 ENDIF 154 ENDIF 155 zmax(9) = REAL( nstop, wp ) ! stop indicator 156 ! !== get global extrema ==! 157 ! !== done by all processes if writting run.stat ==! 129 158 IF( ll_colruns ) THEN 159 zmaxlocal(:) = zmax(:) 130 160 CALL mpp_max( "stpctl", zmax ) ! max over the global domain 131 nstop = NINT( zmax(7) ) ! nstop indicator sheared among all local domains 132 ENDIF 133 ! !== run statistics ==! ("run.stat" files) 161 nstop = NINT( zmax(9) ) ! update nstop indicator (now sheared among all local domains) 162 ENDIF 163 ! !== write "run.stat" files ==! 164 ! !== done only by 1st subdomain at writting timestep ==! 134 165 IF( ll_wrtruns ) THEN 135 166 WRITE(numrun,9500) kt, zmax(1), zmax(2), -zmax(3), zmax(4) 136 istatus = NF90_PUT_VAR( idrun, idssh, (/ zmax(1)/), (/kt/), (/1/) )137 istatus = NF90_PUT_VAR( idrun, idu, (/ zmax(2)/), (/kt/), (/1/) )138 istatus = NF90_PUT_VAR( idrun, ids1, (/-zmax(3)/), (/kt/), (/1/) )139 istatus = NF90_PUT_VAR( idrun, ids2, (/ zmax(4)/), (/kt/), (/1/) )140 istatus = NF90_PUT_VAR( idrun, idt1, (/-zmax(5)/), (/kt/), (/1/) )141 istatus = NF90_PUT_VAR( idrun, idt2, (/ zmax(6)/), (/kt/), (/1/) )167 istatus = NF90_PUT_VAR( nrunid, nvarid(1), (/ zmax(1)/), (/kt/), (/1/) ) 168 istatus = NF90_PUT_VAR( nrunid, nvarid(2), (/ zmax(2)/), (/kt/), (/1/) ) 169 istatus = NF90_PUT_VAR( nrunid, nvarid(3), (/-zmax(3)/), (/kt/), (/1/) ) 170 istatus = NF90_PUT_VAR( nrunid, nvarid(4), (/ zmax(4)/), (/kt/), (/1/) ) 171 istatus = NF90_PUT_VAR( nrunid, nvarid(5), (/-zmax(5)/), (/kt/), (/1/) ) 172 istatus = NF90_PUT_VAR( nrunid, nvarid(6), (/ zmax(6)/), (/kt/), (/1/) ) 142 173 IF( ln_zad_Aimp ) THEN 143 istatus = NF90_PUT_VAR( idrun, idw1, (/ zmax(8)/), (/kt/), (/1/) ) 144 istatus = NF90_PUT_VAR( idrun, idc1, (/ zmax(9)/), (/kt/), (/1/) ) 145 ENDIF 146 IF( MOD( kt , 100 ) == 0 ) istatus = NF90_SYNC(idrun) 147 IF( kt == nitend ) istatus = NF90_CLOSE(idrun) 174 istatus = NF90_PUT_VAR( nrunid, nvarid(7), (/ zmax(7)/), (/kt/), (/1/) ) 175 istatus = NF90_PUT_VAR( nrunid, nvarid(8), (/ zmax(8)/), (/kt/), (/1/) ) 176 ENDIF 177 IF( kt == nitend ) istatus = NF90_CLOSE(nrunid) 148 178 END IF 149 ! !== error handling ==! 150 IF( ( sn_cfctl%l_glochk .OR. lsomeoce ) .AND. ( & ! domain contains some ocean points, check for sensible ranges 151 & zmax(1) > 20._wp .OR. & ! too large sea surface height ( > 20 m ) 152 & zmax(2) > 10._wp .OR. & ! too large velocity ( > 10 m/s) 153 & zmax(3) >= 0._wp .OR. & ! negative or zero sea surface salinity 154 & zmax(4) >= 100._wp .OR. & ! too large sea surface salinity ( > 100 ) 155 & zmax(4) < 0._wp .OR. & ! too large sea surface salinity (keep this line for sea-ice) 156 & ISNAN( zmax(1) + zmax(2) + zmax(3) ) ) ) THEN ! NaN encounter in the tests 157 IF( lk_mpp .AND. sn_cfctl%l_glochk ) THEN 158 ! have use mpp_max (because sn_cfctl%l_glochk=.T. and distributed) 159 CALL mpp_maxloc( 'stpctl', ABS(ssh(:,:,Kmm)) , ssmask(:,:) , zzz, ih ) 160 CALL mpp_maxloc( 'stpctl', ABS(uu(:,:,:,Kmm)) , umask (:,:,:), zzz, iu ) 161 CALL mpp_minloc( 'stpctl', ts(:,:,:,jp_sal,Kmm), tmask (:,:,:), zzz, is1 ) 162 CALL mpp_maxloc( 'stpctl', ts(:,:,:,jp_sal,Kmm), tmask (:,:,:), zzz, is2 ) 179 ! !== error handling ==! 180 ! !== done by all processes at every time step ==! 181 ! 182 IF( zmax(1) > 20._wp .OR. & ! too large sea surface height ( > 20 m ) 183 & zmax(2) > 10._wp .OR. & ! too large velocity ( > 10 m/s) 184 & zmax(3) >= 0._wp .OR. & ! negative or zero sea surface salinity 185 & zmax(4) >= 100._wp .OR. & ! too large sea surface salinity ( > 100 ) 186 & zmax(4) < 0._wp .OR. & ! too large sea surface salinity (keep this line for sea-ice) 187 & ISNAN( zmax(1) + zmax(2) + zmax(3) ) .OR. & ! NaN encounter in the tests 188 & ABS( zmax(1) + zmax(2) + zmax(3) ) > HUGE(1._wp) ) THEN ! Infinity encounter in the tests 189 ! 190 iloc(:,:) = 0 191 IF( ll_colruns ) THEN ! zmax is global, so it is the same on all subdomains -> no dead lock with mpp_maxloc 192 ! first: close the netcdf file, so we can read it 193 IF( lwm .AND. kt /= nitend ) istatus = NF90_CLOSE(nrunid) 194 ! get global loc on the min/max 195 CALL mpp_maxloc( 'stpctl', ABS(ssh(:,:, Kmm)), ssmask(:,: ), zzz, iloc(1:2,1) ) ! mpp_maxloc ok if mask = F 196 CALL mpp_maxloc( 'stpctl', ABS( uu(:,:,:, Kmm)), umask(:,:,:), zzz, iloc(1:3,2) ) 197 CALL mpp_minloc( 'stpctl', ts(:,:,:,jp_sal,Kmm) , tmask(:,:,:), zzz, iloc(1:3,3) ) 198 CALL mpp_maxloc( 'stpctl', ts(:,:,:,jp_sal,Kmm) , tmask(:,:,:), zzz, iloc(1:3,4) ) 199 ! find which subdomain has the max. 200 iareamin(:) = jpnij+1 ; iareamax(:) = 0 ; iareasum(:) = 0 201 DO ji = 1, 9 202 IF( zmaxlocal(ji) == zmax(ji) ) THEN 203 iareamin(ji) = narea ; iareamax(ji) = narea ; iareasum(ji) = 1 204 ENDIF 205 END DO 206 CALL mpp_min( "stpctl", iareamin ) ! min over the global domain 207 CALL mpp_max( "stpctl", iareamax ) ! max over the global domain 208 CALL mpp_sum( "stpctl", iareasum ) ! sum over the global domain 209 ELSE ! find local min and max locations: 210 ! if we are here, this means that the subdomain contains some oce points -> no need to test the mask used in maxloc 211 iloc(1:2,1) = MAXLOC( ABS( ssh(:,:, Kmm)), mask = ssmask(:,: ) == 1._wp ) + (/ nimpp - 1, njmpp - 1 /) 212 iloc(1:3,2) = MAXLOC( ABS( uu(:,:,:, Kmm)), mask = umask(:,:,:) == 1._wp ) + (/ nimpp - 1, njmpp - 1, 0 /) 213 iloc(1:3,3) = MINLOC( ts(:,:,:,jp_sal,Kmm) , mask = tmask(:,:,:) == 1._wp ) + (/ nimpp - 1, njmpp - 1, 0 /) 214 iloc(1:3,4) = MAXLOC( ts(:,:,:,jp_sal,Kmm) , mask = tmask(:,:,:) == 1._wp ) + (/ nimpp - 1, njmpp - 1, 0 /) 215 iareamin(:) = narea ; iareamax(:) = narea ; iareasum(:) = 1 ! this is local information 216 ENDIF 217 ! 218 WRITE(ctmp1,*) ' stp_ctl: |ssh| > 20 m or |U| > 10 m/s or S <= 0 or S >= 100 or NaN encounter in the tests' 219 CALL wrt_line( ctmp2, kt, '|ssh| max', zmax(1), iloc(:,1), iareasum(1), iareamin(1), iareamax(1) ) 220 CALL wrt_line( ctmp3, kt, '|U| max', zmax(2), iloc(:,2), iareasum(2), iareamin(2), iareamax(2) ) 221 CALL wrt_line( ctmp4, kt, 'Sal min', -zmax(3), iloc(:,3), iareasum(3), iareamin(3), iareamax(3) ) 222 CALL wrt_line( ctmp5, kt, 'Sal max', zmax(4), iloc(:,4), iareasum(4), iareamin(4), iareamax(4) ) 223 IF( Agrif_Root() ) THEN 224 WRITE(ctmp6,*) ' ===> output of last computed fields in output.abort* files' 163 225 ELSE 164 ! find local min and max locations 165 ih(:) = MAXLOC( ABS( ssh(:,:,Kmm) ) ) + (/ nimpp - 1, njmpp - 1 /) 166 iu(:) = MAXLOC( ABS( uu (:,:,:,Kmm) ) ) + (/ nimpp - 1, njmpp - 1, 0 /) 167 is1(:) = MINLOC( ts(:,:,:,jp_sal,Kmm), mask = tmask(:,:,:) == 1._wp ) + (/ nimpp - 1, njmpp - 1, 0 /) 168 is2(:) = MAXLOC( ts(:,:,:,jp_sal,Kmm), mask = tmask(:,:,:) == 1._wp ) + (/ nimpp - 1, njmpp - 1, 0 /) 169 ENDIF 170 171 WRITE(ctmp1,*) ' stp_ctl: |ssh| > 20 m or |U| > 10 m/s or S <= 0 or S >= 100 or NaN encounter in the tests' 172 WRITE(ctmp2,9100) kt, zmax(1), ih(1) , ih(2) 173 WRITE(ctmp3,9200) kt, zmax(2), iu(1) , iu(2) , iu(3) 174 WRITE(ctmp4,9300) kt, - zmax(3), is1(1), is1(2), is1(3) 175 WRITE(ctmp5,9400) kt, zmax(4), is2(1), is2(2), is2(3) 176 WRITE(ctmp6,*) ' ===> output of last computed fields in output.abort.nc file' 177 226 WRITE(ctmp6,*) ' ===> output of last computed fields in '//TRIM(Agrif_CFixed())//'_output.abort* files' 227 ENDIF 228 ! 178 229 CALL dia_wri_state( Kmm, 'output.abort' ) ! create an output.abort file 179 180 IF( .NOT. sn_cfctl%l_glochk ) THEN181 WRITE(ctmp8,*) 'E R R O R message from sub-domain: ', narea182 CALL ctl_stop( 'STOP', ctmp1, ' ', ctmp8, ' ', ctmp2, ctmp3, ctmp4, ctmp5, ctmp6)183 ELSE184 CALL ctl_stop( ctmp1, ' ', ctmp2, ctmp3, ctmp4, ctmp5, ' ', ctmp6, ' ' )185 ENDIF186 187 kindic = -3188 !189 ENDIF190 !191 9100 FORMAT (' kt=',i8,' |ssh| max: ',1pg11.4,', at i j : ',2i5) 192 9200 FORMAT (' kt=',i8,' |U| max: ',1pg11.4,', at i j k: ',3i5) 193 9300 FORMAT (' kt=',i8,' S min: ',1pg11.4,', at i j k: ',3i5) 194 9400 FORMAT (' kt=',i8,' S max: ',1pg11.4,', at i j k: ',3i5) 230 ! 231 IF( ll_colruns .or. jpnij == 1 ) THEN ! all processes synchronized -> use lwp to print in opened ocean.output files 232 IF(lwp) THEN ; CALL ctl_stop( ctmp1, ' ', ctmp2, ctmp3, ctmp4, ctmp5, ' ', ctmp6 ) 233 ELSE ; nstop = MAX(1, nstop) ! make sure nstop > 0 (automatically done when calling ctl_stop) 234 ENDIF 235 ELSE ! only mpi subdomains with errors are here -> STOP now 236 CALL ctl_stop( 'STOP', ctmp1, ' ', ctmp2, ctmp3, ctmp4, ctmp5, ' ', ctmp6 ) 237 ENDIF 238 ! 239 ENDIF 240 ! 241 IF( nstop > 0 ) THEN ! an error was detected and we did not abort yet... 242 ngrdstop = Agrif_Fixed() ! store which grid got this error 243 IF( .NOT. ll_colruns .AND. jpnij > 1 ) CALL ctl_stop( 'STOP' ) ! we must abort here to avoid MPI deadlock 244 ENDIF 245 ! 195 246 9500 FORMAT(' it :', i8, ' |ssh|_max: ', D23.16, ' |U|_max: ', D23.16,' S_min: ', D23.16,' S_max: ', D23.16) 196 247 ! 197 248 END SUBROUTINE stp_ctl 249 250 251 SUBROUTINE wrt_line( cdline, kt, cdprefix, pval, kloc, ksum, kmin, kmax ) 252 !!---------------------------------------------------------------------- 253 !! *** ROUTINE wrt_line *** 254 !! 255 !! ** Purpose : write information line 256 !! 257 !!---------------------------------------------------------------------- 258 CHARACTER(len=*), INTENT( out) :: cdline 259 CHARACTER(len=*), INTENT(in ) :: cdprefix 260 REAL(wp), INTENT(in ) :: pval 261 INTEGER, DIMENSION(3), INTENT(in ) :: kloc 262 INTEGER, INTENT(in ) :: kt, ksum, kmin, kmax 263 ! 264 CHARACTER(len=80) :: clsuff 265 CHARACTER(len=9 ) :: clkt, clsum, clmin, clmax 266 CHARACTER(len=9 ) :: cli, clj, clk 267 CHARACTER(len=1 ) :: clfmt 268 CHARACTER(len=4 ) :: cl4 ! needed to be able to compile with Agrif, I don't know why 269 INTEGER :: ifmtk 270 !!---------------------------------------------------------------------- 271 WRITE(clkt , '(i9)') kt 272 273 WRITE(clfmt, '(i1)') INT(LOG10(REAL(jpnij ,wp))) + 1 ! how many digits to we need to write ? (we decide max = 9) 274 !!! WRITE(clsum, '(i'//clfmt//')') ksum ! this is creating a compilation error with AGRIF 275 cl4 = '(i'//clfmt//')' ; WRITE(clsum, cl4) ksum 276 WRITE(clfmt, '(i1)') INT(LOG10(REAL(MAX(1,jpnij-1),wp))) + 1 ! how many digits to we need to write ? (we decide max = 9) 277 cl4 = '(i'//clfmt//')' ; WRITE(clmin, cl4) kmin-1 278 WRITE(clmax, cl4) kmax-1 279 ! 280 WRITE(clfmt, '(i1)') INT(LOG10(REAL(jpiglo,wp))) + 1 ! how many digits to we need to write jpiglo? (we decide max = 9) 281 cl4 = '(i'//clfmt//')' ; WRITE(cli, cl4) kloc(1) ! this is ok with AGRIF 282 WRITE(clfmt, '(i1)') INT(LOG10(REAL(jpjglo,wp))) + 1 ! how many digits to we need to write jpjglo? (we decide max = 9) 283 cl4 = '(i'//clfmt//')' ; WRITE(clj, cl4) kloc(2) ! this is ok with AGRIF 284 ! 285 IF( ksum == 1 ) THEN ; WRITE(clsuff,9100) TRIM(clmin) 286 ELSE ; WRITE(clsuff,9200) TRIM(clsum), TRIM(clmin), TRIM(clmax) 287 ENDIF 288 IF(kloc(3) == 0) THEN 289 ifmtk = INT(LOG10(REAL(jpk,wp))) + 1 ! how many digits to we need to write jpk? (we decide max = 9) 290 clk = REPEAT(' ', ifmtk) ! create the equivalent in blank string 291 WRITE(cdline,9300) TRIM(ADJUSTL(clkt)), TRIM(ADJUSTL(cdprefix)), pval, TRIM(cli), TRIM(clj), clk(1:ifmtk), TRIM(clsuff) 292 ELSE 293 WRITE(clfmt, '(i1)') INT(LOG10(REAL(jpk,wp))) + 1 ! how many digits to we need to write jpk? (we decide max = 9) 294 !!! WRITE(clk, '(i'//clfmt//')') kloc(3) ! this is creating a compilation error with AGRIF 295 cl4 = '(i'//clfmt//')' ; WRITE(clk, cl4) kloc(3) ! this is ok with AGRIF 296 WRITE(cdline,9400) TRIM(ADJUSTL(clkt)), TRIM(ADJUSTL(cdprefix)), pval, TRIM(cli), TRIM(clj), TRIM(clk), TRIM(clsuff) 297 ENDIF 298 ! 299 9100 FORMAT('MPI rank ', a) 300 9200 FORMAT('found in ', a, ' MPI tasks, spread out among ranks ', a, ' to ', a) 301 9300 FORMAT('kt ', a, ' ', a, ' ', 1pg11.4, ' at i j ', a, ' ', a, ' ', a, ' ', a) 302 9400 FORMAT('kt ', a, ' ', a, ' ', 1pg11.4, ' at i j k ', a, ' ', a, ' ', a, ' ', a) 303 ! 304 END SUBROUTINE wrt_line 305 198 306 199 307 !!====================================================================== -
NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/src/OFF/nemogcm.F90
r12779 r13193 31 31 USE domqco ! tools for scale factor (dom_qco_r3c routine) 32 32 #endif 33 USE bdy_oce, ONLY : ln_bdy 34 USE bdyini ! open boundary cond. setting (bdy_init routine) 33 35 ! ! ocean physics 34 36 USE ldftra ! lateral diffusivity setting (ldf_tra_init routine) … … 93 95 !! Madec, 2008, internal report, IPSL. 94 96 !!---------------------------------------------------------------------- 95 INTEGER :: istp , indic! time step index97 INTEGER :: istp ! time step index 96 98 !!---------------------------------------------------------------------- 97 99 … … 145 147 # endif 146 148 #endif 147 CALL stp_ctl ( istp , indic )! Time loop: control and print149 CALL stp_ctl ( istp ) ! Time loop: control and print 148 150 istp = istp + 1 149 151 END DO … … 160 162 IF( nstop /= 0 .AND. lwp ) THEN ! error print 161 163 WRITE(ctmp1,*) ' ==>>> nemo_gcm: a total of ', nstop, ' errors have been found' 162 CALL ctl_stop( ctmp1 ) 164 WRITE(ctmp2,*) ' Look for "E R R O R" messages in all existing ocean_output* files' 165 CALL ctl_stop( ' ', ctmp1, ' ', ctmp2 ) 163 166 ENDIF 164 167 ! … … 242 245 ! 243 246 ! finalize the definition of namctl variables 244 IF( sn_cfctl%l_allon ) THEN 245 ! Turn on all options. 246 CALL nemo_set_cfctl( sn_cfctl, .TRUE., .TRUE. ) 247 ! Ensure all processors are active 248 sn_cfctl%procmin = 0 ; sn_cfctl%procmax = 1000000 ; sn_cfctl%procincr = 1 249 ELSEIF( sn_cfctl%l_config ) THEN 250 ! Activate finer control of report outputs 251 ! optionally switch off output from selected areas (note this only 252 ! applies to output which does not involve global communications) 253 IF( ( narea < sn_cfctl%procmin .OR. narea > sn_cfctl%procmax ) .OR. & 254 & ( MOD( narea - sn_cfctl%procmin, sn_cfctl%procincr ) /= 0 ) ) & 255 & CALL nemo_set_cfctl( sn_cfctl, .FALSE., .FALSE. ) 256 ELSE 257 ! turn off all options. 258 CALL nemo_set_cfctl( sn_cfctl, .FALSE., .TRUE. ) 259 ENDIF 247 IF( narea < sn_cfctl%procmin .OR. narea > sn_cfctl%procmax .OR. MOD( narea - sn_cfctl%procmin, sn_cfctl%procincr ) /= 0 ) & 248 & CALL nemo_set_cfctl( sn_cfctl, .FALSE. ) 260 249 ! 261 250 lwp = (narea == 1) .OR. sn_cfctl%l_oceout ! control of all listing output print … … 322 311 ! Initialise time level indices 323 312 Nbb = 1; Nnn = 2; Naa = 3; Nrhs = Naa 324 325 313 326 314 ! !-------------------------------! … … 344 332 345 333 CALL sbc_init( Nbb, Nnn, Naa ) ! Forcings : surface module 334 CALL bdy_init ! Open boundaries initialisation 346 335 347 336 ! ! Tracer physics … … 386 375 WRITE(numout,*) '~~~~~~~~' 387 376 WRITE(numout,*) ' Namelist namctl' 388 WRITE(numout,*) ' sn_cfctl%l_glochk = ', sn_cfctl%l_glochk389 WRITE(numout,*) ' sn_cfctl%l_allon = ', sn_cfctl%l_allon390 WRITE(numout,*) ' finer control over o/p sn_cfctl%l_config = ', sn_cfctl%l_config391 377 WRITE(numout,*) ' sn_cfctl%l_runstat = ', sn_cfctl%l_runstat 392 378 WRITE(numout,*) ' sn_cfctl%l_trcstat = ', sn_cfctl%l_trcstat … … 507 493 USE zdf_oce, ONLY : zdf_oce_alloc 508 494 USE trc_oce, ONLY : trc_oce_alloc 495 USE bdy_oce, ONLY : bdy_oce_alloc 509 496 ! 510 497 INTEGER :: ierr … … 516 503 ierr = ierr + zdf_oce_alloc() ! ocean vertical physics 517 504 ierr = ierr + trc_oce_alloc() ! shared TRC / TRA arrays 505 ierr = ierr + bdy_oce_alloc() ! bdy masks (incl. initialization) 518 506 ! 519 507 CALL mpp_sum( 'nemogcm', ierr ) … … 522 510 END SUBROUTINE nemo_alloc 523 511 524 SUBROUTINE nemo_set_cfctl(sn_cfctl, setto , for_all)512 SUBROUTINE nemo_set_cfctl(sn_cfctl, setto ) 525 513 !!---------------------------------------------------------------------- 526 514 !! *** ROUTINE nemo_set_cfctl *** 527 515 !! 528 516 !! ** Purpose : Set elements of the output control structure to setto. 529 !! for_all should be .false. unless all areas are to be 530 !! treated identically. 531 !! 517 !! 532 518 !! ** Method : Note this routine can be used to switch on/off some 533 !! types of output for selected areas but any output types 534 !! that involve global communications (e.g. mpp_max, glob_sum) 535 !! should be protected from selective switching by the 536 !! for_all argument 537 !!---------------------------------------------------------------------- 538 LOGICAL :: setto, for_all 539 TYPE(sn_ctl) :: sn_cfctl 540 !!---------------------------------------------------------------------- 541 IF( for_all ) THEN 542 sn_cfctl%l_runstat = setto 543 sn_cfctl%l_trcstat = setto 544 ENDIF 519 !! types of output for selected areas. 520 !!---------------------------------------------------------------------- 521 TYPE(sn_ctl), INTENT(inout) :: sn_cfctl 522 LOGICAL , INTENT(in ) :: setto 523 !!---------------------------------------------------------------------- 524 sn_cfctl%l_runstat = setto 525 sn_cfctl%l_trcstat = setto 545 526 sn_cfctl%l_oceout = setto 546 527 sn_cfctl%l_layout = setto … … 572 553 573 554 574 SUBROUTINE stp_ctl( kt , kindic)555 SUBROUTINE stp_ctl( kt ) 575 556 !!---------------------------------------------------------------------- 576 557 !! *** ROUTINE stp_ctl *** … … 583 564 !!---------------------------------------------------------------------- 584 565 INTEGER, INTENT(in ) :: kt ! ocean time-step index 585 INTEGER, INTENT(inout) :: kindic ! indicator of solver convergence586 566 !!---------------------------------------------------------------------- 587 567 ! -
NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/src/SAO/nemogcm.F90
r12724 r13193 158 158 ! 159 159 ! finalize the definition of namctl variables 160 IF( sn_cfctl%l_allon ) THEN 161 ! Turn on all options. 162 CALL nemo_set_cfctl( sn_cfctl, .TRUE., .TRUE. ) 163 ! Ensure all processors are active 164 sn_cfctl%procmin = 0 ; sn_cfctl%procmax = 1000000 ; sn_cfctl%procincr = 1 165 ELSEIF( sn_cfctl%l_config ) THEN 166 ! Activate finer control of report outputs 167 ! optionally switch off output from selected areas (note this only 168 ! applies to output which does not involve global communications) 169 IF( ( narea < sn_cfctl%procmin .OR. narea > sn_cfctl%procmax ) .OR. & 170 & ( MOD( narea - sn_cfctl%procmin, sn_cfctl%procincr ) /= 0 ) ) & 171 & CALL nemo_set_cfctl( sn_cfctl, .FALSE., .FALSE. ) 172 ELSE 173 ! turn off all options. 174 CALL nemo_set_cfctl( sn_cfctl, .FALSE., .TRUE. ) 175 ENDIF 160 IF( narea < sn_cfctl%procmin .OR. narea > sn_cfctl%procmax .OR. MOD( narea - sn_cfctl%procmin, sn_cfctl%procincr ) /= 0 ) & 161 & CALL nemo_set_cfctl( sn_cfctl, .FALSE. ) 176 162 ! 177 163 lwp = (narea == 1) .OR. sn_cfctl%l_oceout ! control of all listing output print … … 270 256 WRITE(numout,*) '~~~~~~~~' 271 257 WRITE(numout,*) ' Namelist namctl' 272 WRITE(numout,*) ' sn_cfctl%l_glochk = ', sn_cfctl%l_glochk273 WRITE(numout,*) ' sn_cfctl%l_allon = ', sn_cfctl%l_allon274 WRITE(numout,*) ' finer control over o/p sn_cfctl%l_config = ', sn_cfctl%l_config275 258 WRITE(numout,*) ' sn_cfctl%l_runstat = ', sn_cfctl%l_runstat 276 259 WRITE(numout,*) ' sn_cfctl%l_trcstat = ', sn_cfctl%l_trcstat … … 410 393 END SUBROUTINE nemo_alloc 411 394 412 SUBROUTINE nemo_set_cfctl(sn_cfctl, setto , for_all)395 SUBROUTINE nemo_set_cfctl(sn_cfctl, setto ) 413 396 !!---------------------------------------------------------------------- 414 397 !! *** ROUTINE nemo_set_cfctl *** 415 398 !! 416 399 !! ** Purpose : Set elements of the output control structure to setto. 417 !! for_all should be .false. unless all areas are to be418 !! treated identically.419 400 !! 420 401 !! ** Method : Note this routine can be used to switch on/off some 421 !! types of output for selected areas but any output types 422 !! that involve global communications (e.g. mpp_max, glob_sum) 423 !! should be protected from selective switching by the 424 !! for_all argument 425 !!---------------------------------------------------------------------- 426 LOGICAL :: setto, for_all 427 TYPE(sn_ctl) :: sn_cfctl 428 !!---------------------------------------------------------------------- 429 IF( for_all ) THEN 430 sn_cfctl%l_runstat = setto 431 sn_cfctl%l_trcstat = setto 432 ENDIF 402 !! types of output for selected areas. 403 !!---------------------------------------------------------------------- 404 TYPE(sn_ctl), INTENT(inout) :: sn_cfctl 405 LOGICAL , INTENT(in ) :: setto 406 !!---------------------------------------------------------------------- 407 sn_cfctl%l_runstat = setto 408 sn_cfctl%l_trcstat = setto 433 409 sn_cfctl%l_oceout = setto 434 410 sn_cfctl%l_layout = setto -
NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/src/SAS/diawri.F90
r12724 r13193 138 138 !! Each nn_write time step, output the instantaneous or mean fields 139 139 !!---------------------------------------------------------------------- 140 !!141 140 INTEGER, INTENT( in ) :: kt ! ocean time-step index 142 INTEGER, INTENT( in ) :: Kmm ! ocean time level index141 INTEGER, INTENT( in ) :: Kmm ! ocean time level index 143 142 !! 144 143 LOGICAL :: ll_print = .FALSE. ! =T print and flush numout … … 462 461 CALL iom_close( inum ) 463 462 ENDIF 464 #endif 465 463 ! 464 #endif 466 465 END SUBROUTINE dia_wri_state 467 466 -
NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/src/SAS/nemogcm.F90
r12724 r13193 126 126 END DO 127 127 ! 128 IF( .NOT. Agrif_Root() ) THEN129 CALL Agrif_ParentGrid_To_ChildGrid()130 IF( ln_timing ) CALL timing_finalize131 CALL Agrif_ChildGrid_To_ParentGrid()132 ENDIF133 !134 128 #else 135 129 ! … … 166 160 IF( nstop /= 0 .AND. lwp ) THEN ! error print 167 161 WRITE(ctmp1,*) ' ==>>> nemo_gcm: a total of ', nstop, ' errors have been found' 168 CALL ctl_stop( ctmp1 ) 162 IF( ngrdstop > 0 ) THEN 163 WRITE(ctmp9,'(i2)') ngrdstop 164 WRITE(ctmp2,*) ' E R R O R detected in Agrif grid '//TRIM(ctmp9) 165 WRITE(ctmp3,*) ' Look for "E R R O R" messages in all existing '//TRIM(ctmp9)//'_ocean_output* files' 166 CALL ctl_stop( ' ', ctmp1, ' ', ctmp2, ' ', ctmp3 ) 167 ELSE 168 WRITE(ctmp2,*) ' Look for "E R R O R" messages in all existing ocean_output* files' 169 CALL ctl_stop( ' ', ctmp1, ' ', ctmp2 ) 170 ENDIF 169 171 ENDIF 170 172 ! … … 275 277 ! 276 278 ! finalize the definition of namctl variables 277 IF( sn_cfctl%l_allon ) THEN 278 ! Turn on all options. 279 CALL nemo_set_cfctl( sn_cfctl, .TRUE., .TRUE. ) 280 ! Ensure all processors are active 281 sn_cfctl%procmin = 0 ; sn_cfctl%procmax = 1000000 ; sn_cfctl%procincr = 1 282 ELSEIF( sn_cfctl%l_config ) THEN 283 ! Activate finer control of report outputs 284 ! optionally switch off output from selected areas (note this only 285 ! applies to output which does not involve global communications) 286 IF( ( narea < sn_cfctl%procmin .OR. narea > sn_cfctl%procmax ) .OR. & 287 & ( MOD( narea - sn_cfctl%procmin, sn_cfctl%procincr ) /= 0 ) ) & 288 & CALL nemo_set_cfctl( sn_cfctl, .FALSE., .FALSE. ) 289 ELSE 290 ! turn off all options. 291 CALL nemo_set_cfctl( sn_cfctl, .FALSE., .TRUE. ) 292 ENDIF 279 IF( narea < sn_cfctl%procmin .OR. narea > sn_cfctl%procmax .OR. MOD( narea - sn_cfctl%procmin, sn_cfctl%procincr ) /= 0 ) & 280 & CALL nemo_set_cfctl( sn_cfctl, .FALSE. ) 293 281 ! 294 282 lwp = (narea == 1) .OR. sn_cfctl%l_oceout ! control of all listing output print … … 408 396 WRITE(numout,*) '~~~~~~~~' 409 397 WRITE(numout,*) ' Namelist namctl' 410 WRITE(numout,*) ' sn_cfctl%l_glochk = ', sn_cfctl%l_glochk411 WRITE(numout,*) ' sn_cfctl%l_allon = ', sn_cfctl%l_allon412 WRITE(numout,*) ' finer control over o/p sn_cfctl%l_config = ', sn_cfctl%l_config413 398 WRITE(numout,*) ' sn_cfctl%l_runstat = ', sn_cfctl%l_runstat 414 399 WRITE(numout,*) ' sn_cfctl%l_trcstat = ', sn_cfctl%l_trcstat … … 552 537 END SUBROUTINE nemo_alloc 553 538 554 SUBROUTINE nemo_set_cfctl(sn_cfctl, setto , for_all)539 SUBROUTINE nemo_set_cfctl(sn_cfctl, setto ) 555 540 !!---------------------------------------------------------------------- 556 541 !! *** ROUTINE nemo_set_cfctl *** 557 542 !! 558 543 !! ** Purpose : Set elements of the output control structure to setto. 559 !! for_all should be .false. unless all areas are to be560 !! treated identically.561 544 !! 562 545 !! ** Method : Note this routine can be used to switch on/off some 563 !! types of output for selected areas but any output types 564 !! that involve global communications (e.g. mpp_max, glob_sum) 565 !! should be protected from selective switching by the 566 !! for_all argument 567 !!---------------------------------------------------------------------- 568 LOGICAL :: setto, for_all 569 TYPE(sn_ctl) :: sn_cfctl 570 !!---------------------------------------------------------------------- 571 IF( for_all ) THEN 572 sn_cfctl%l_runstat = setto 573 sn_cfctl%l_trcstat = setto 574 ENDIF 546 !! types of output for selected areas. 547 !!---------------------------------------------------------------------- 548 TYPE(sn_ctl), INTENT(inout) :: sn_cfctl 549 LOGICAL , INTENT(in ) :: setto 550 !!---------------------------------------------------------------------- 551 sn_cfctl%l_runstat = setto 552 sn_cfctl%l_trcstat = setto 575 553 sn_cfctl%l_oceout = setto 576 554 sn_cfctl%l_layout = setto -
NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/src/SAS/step.F90
r12724 r13193 74 74 !! -2- Outputs and diagnostics 75 75 !!---------------------------------------------------------------------- 76 INTEGER :: indic ! error indicator if < 077 !! ---------------------------------------------------------------------78 76 79 77 #if defined key_agrif 80 IF( nstop > 0 ) return ! avoid to go further if an error was detected during previous time step78 IF( nstop > 0 ) RETURN ! avoid to go further if an error was detected during previous time step (child grid) 81 79 kstp = nit000 + Agrif_Nb_Step() 82 80 Kbb_a = Nbb; Kmm_a = Nnn; Krhs_a = Nrhs ! agrif_oce module copies of time level indices 83 IF 84 IF ( Agrif_Root() .and. lwp) Write(*,*) '---'85 IF (lwp) Write(*,*) 'Grid Number',Agrif_Fixed(),' time step ',kstp, 'int tstep',Agrif_NbStepint()81 IF( lk_agrif_debug ) THEN 82 IF( Agrif_Root() .and. lwp) WRITE(*,*) '---' 83 IF(lwp) WRITE(*,*) 'Grid Number', Agrif_Fixed(),' time step ', kstp, 'int tstep', Agrif_NbStepint() 86 84 ENDIF 87 88 IF ( kstp == (nit000 + 1) ) lk_agrif_fstep = .FALSE. 89 85 IF( kstp == nit000 + 1 ) lk_agrif_fstep = .FALSE. 90 86 # if defined key_iomput 91 87 IF( Agrif_Nbstepint() == 0 ) CALL iom_swap( cxios_context ) 92 88 # endif 93 89 #endif 94 indic = 0 ! although indic is not changed in stp_ctl95 ! need to keep the same interface96 90 IF( kstp == nit000 ) CALL iom_init( cxios_context ) ! iom_put initialization (must be done after nemo_init for AGRIF+XIOS+OASIS) 97 91 IF( kstp /= nit000 ) CALL day( kstp ) ! Calendar (day was already called at nit000 in day_init) … … 112 106 ! AGRIF recursive integration 113 107 !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 114 CALL Agrif_Integrate_ChildGrids( stp ) 115 #endif 108 CALL Agrif_Integrate_ChildGrids( stp ) 116 109 110 #endif 117 111 !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 118 112 ! Control 119 113 !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 120 CALL stp_ctl( kstp, indic ) 121 IF( indic < 0 ) THEN 122 CALL ctl_stop( 'step: indic < 0' ) 123 CALL dia_wri_state( Nnn, 'output.abort' ) 124 ENDIF 114 CALL stp_ctl( kstp, Nnn ) 115 125 116 #if defined key_agrif 126 117 !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> … … 132 123 #endif 133 124 ENDIF 125 134 126 #endif 135 127 !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> … … 141 133 ! Coupled mode 142 134 !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 143 IF( lk_oasis ) CALL sbc_cpl_snd( kstp, Nbb, Nnn )! coupled mode : field exchanges if OASIS-coupled ice135 IF( lk_oasis .AND. nstop == 0 ) CALL sbc_cpl_snd( kstp, Nbb, Nnn ) ! coupled mode : field exchanges if OASIS-coupled ice 144 136 145 137 #if defined key_iomput … …