- Timestamp:
- 2019-10-29T18:14:49+01:00 (4 years ago)
- Location:
- NEMO/branches/2019/dev_r11085_ASINTER-05_Brodeau_Advanced_Bulk/tests/ISOMIP
- Files:
-
- 3 edited
-
EXPREF/context_nemo.xml (modified) (2 diffs)
-
EXPREF/namelist_cfg (modified) (2 diffs)
-
MY_SRC/usrdef_nam.F90 (modified) (4 diffs)
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/dev_r11085_ASINTER-05_Brodeau_Advanced_Bulk/tests/ISOMIP/EXPREF/context_nemo.xml
r9930 r11831 6 6 <context id="nemo"> 7 7 <!-- $id$ --> 8 <variable_definition> 9 <!-- Year of time origin for NetCDF files; defaults to 1800 --> 10 <variable id="ref_year" type="int" > 1800 </variable> 11 <variable id="rau0" type="float" > 1026.0 </variable> 12 <variable id="cpocean" type="float" > 3991.86795711963 </variable> 13 <variable id="convSpsu" type="float" > 0.99530670233846 </variable> 14 <variable id="rhoic" type="float" > 917.0 </variable> 15 <variable id="rhosn" type="float" > 330.0 </variable> 16 <variable id="missval" type="float" > 1.e20 </variable> 17 </variable_definition> 8 18 <!-- Fields definition --> 9 19 <field_definition src="./field_def_nemo-oce.xml"/> <!-- NEMO ocean dynamics --> … … 18 28 19 29 <axis_definition> 20 <axis id="deptht" long_name="Vertical T levels" unit="m" positive="down" />21 <axis id="depthu" long_name="Vertical U levels" unit="m" positive="down" />22 <axis id="depthv" long_name="Vertical V levels" unit="m" positive="down" />23 <axis id="depthw" long_name="Vertical W levels" unit="m" positive="down" />24 <axis id="nfloat" long_name="Float number" unit="-" />25 <axis id="icbcla" long_name="Iceberg class" unit="1"/>26 <axis id="ncatice" long_name="Ice category" unit="1"/>27 <axis id="iax_20C" long_name="20 degC isotherm" unit="degC"/>28 <axis id="iax_28C" long_name="28 degC isotherm" unit="degC"/>30 <axis id="deptht" long_name="Vertical T levels" unit="m" positive="down" /> 31 <axis id="depthu" long_name="Vertical U levels" unit="m" positive="down" /> 32 <axis id="depthv" long_name="Vertical V levels" unit="m" positive="down" /> 33 <axis id="depthw" long_name="Vertical W levels" unit="m" positive="down" /> 34 <axis id="nfloat" long_name="Float number" unit="-" /> 35 <axis id="icbcla" long_name="Iceberg class" unit="1" /> 36 <axis id="ncatice" long_name="Ice category" unit="1" /> 37 <axis id="iax_20C" long_name="20 degC isotherm" unit="degC" /> 38 <axis id="iax_28C" long_name="28 degC isotherm" unit="degC" /> 29 39 </axis_definition> 30 40 -
NEMO/branches/2019/dev_r11085_ASINTER-05_Brodeau_Advanced_Bulk/tests/ISOMIP/EXPREF/namelist_cfg
r10075 r11831 78 78 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! filename ! 79 79 ! nn_isf == 4 80 sn_fwfisf = 'rnfisf' , -12 80 sn_fwfisf = 'rnfisf' , -12. ,'sowflisf', .false. , .true. , 'yearly' , '' , '' , '' 81 81 ! nn_isf == 3 82 sn_rnfisf = 'rnfisf' , -12 82 sn_rnfisf = 'rnfisf' , -12. ,'sofwfisf', .false. , .true. , 'yearly' , '' , '' , '' 83 83 ! nn_isf == 2 and 3 84 sn_depmax_isf = 'rnfisf' , -12 85 sn_depmin_isf = 'rnfisf' , -12 84 sn_depmax_isf = 'rnfisf' , -12. ,'sozisfmax' , .false. , .true. , 'yearly' , '' , '' , '' 85 sn_depmin_isf = 'rnfisf' , -12. ,'sozisfmin' , .false. , .true. , 'yearly' , '' , '' , '' 86 86 ! nn_isf == 2 87 sn_Leff_isf = 'rnfisf' , -12 87 sn_Leff_isf = 'rnfisf' , -12. ,'Leff' , .false. , .true. , 'yearly' , '' , '' , '' 88 88 ! for all case 89 89 nn_isf = 1 ! ice shelf melting/freezing … … 272 272 !! namdiu Cool skin and warm layer models (default: OFF) 273 273 !! namdiu Cool skin and warm layer models (default: OFF) 274 !! namflo float parameters ( "key_float")275 !! nam_diaharm Harmonic analysis of tidal constituents ( "key_diaharm")276 !! nam dct transports through some sections ("key_diadct")274 !! namflo float parameters (default: OFF) 275 !! nam_diaharm Harmonic analysis of tidal constituents (default: OFF) 276 !! nam_diadct transports through some sections (default: OFF) 277 277 !! nam_diatmb Top Middle Bottom Output (default: OFF) 278 278 !! nam_dia25h 25h Mean Output (default: OFF) -
NEMO/branches/2019/dev_r11085_ASINTER-05_Brodeau_Advanced_Bulk/tests/ISOMIP/MY_SRC/usrdef_nam.F90
r10074 r11831 42 42 CONTAINS 43 43 44 SUBROUTINE usr_def_nam( ldtxt, ldnam,cd_cfg, kk_cfg, kpi, kpj, kpk, kperio )44 SUBROUTINE usr_def_nam( cd_cfg, kk_cfg, kpi, kpj, kpk, kperio ) 45 45 !!---------------------------------------------------------------------- 46 46 !! *** ROUTINE dom_nam *** … … 54 54 !! ** input : - namusr_def namelist found in namelist_cfg 55 55 !!---------------------------------------------------------------------- 56 CHARACTER(len=*), DIMENSION(:), INTENT(out) :: ldtxt, ldnam ! stored print information57 56 CHARACTER(len=*) , INTENT(out) :: cd_cfg ! configuration name 58 57 INTEGER , INTENT(out) :: kk_cfg ! configuration resolution … … 60 59 INTEGER , INTENT(out) :: kperio ! lateral global domain b.c. 61 60 ! 62 INTEGER :: ios , ii! Local integer61 INTEGER :: ios ! Local integer 63 62 !! 64 63 NAMELIST/namusr_def/ ln_zco, ln_zps, ln_sco, rn_e1deg, rn_e2deg, rn_e3 65 64 !!---------------------------------------------------------------------- 66 65 ! 67 ii = 168 !69 66 REWIND( numnam_cfg ) ! Namelist namusr_def (exist in namelist_cfg only) 70 67 READ ( numnam_cfg, namusr_def, IOSTAT = ios, ERR = 902 ) 71 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namusr_def in configuration namelist' , .TRUE.)68 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namusr_def in configuration namelist' ) 72 69 ! 73 WRITE( ldnam(:), namusr_def )70 IF(lwm) WRITE( numond, namusr_def ) 74 71 ! 75 72 cd_cfg = 'ISOMIP' ! name & resolution (not used) … … 81 78 kpk = INT( rbathy / rn_e3 ) + 1 ! add 1 for t-point in the seafloor 82 79 ! 83 ! ! control print84 WRITE(ldtxt(ii),*) ' ' ; ii = ii + 185 WRITE(ldtxt(ii),*) 'usr_def_nam : read the user defined namelist (namusr_def) in namelist_cfg' ; ii = ii + 186 WRITE(ldtxt(ii),*) '~~~~~~~~~~~ ' ; ii = ii + 187 WRITE(ldtxt(ii),*) ' Namelist namusr_def : ISOMIP test case' ; ii = ii + 188 WRITE(ldtxt(ii),*) ' type of vertical coordinate : ' ; ii = ii + 189 WRITE(ldtxt(ii),*) ' z-coordinate flag ln_zco = ', ln_zco ; ii = ii + 190 WRITE(ldtxt(ii),*) ' z-partial-step coordinate flag ln_zps = ', ln_zps ; ii = ii + 191 WRITE(ldtxt(ii),*) ' s-coordinate flag ln_sco = ', ln_sco ; ii = ii + 192 WRITE(ldtxt(ii),*) ' resolution' ; ii = ii + 193 WRITE(ldtxt(ii),*) ' zonal resolution rn_e1deg = ', rn_e1deg, ' degrees' ; ii = ii + 194 WRITE(ldtxt(ii),*) ' meridional resolution rn_e1deg = ', rn_e1deg, ' degrees' ; ii = ii + 195 WRITE(ldtxt(ii),*) ' vertical resolution rn_e3 = ', rn_e3 , ' meters' ; ii = ii + 196 WRITE(ldtxt(ii),*) ' ISOMIP domain = 15° x 10° x 900 m' ; ii = ii + 197 WRITE(ldtxt(ii),*) ' resulting global domain size : jpiglo = ', kpi ; ii = ii + 198 WRITE(ldtxt(ii),*) ' jpjglo = ', kpj ; ii = ii + 199 WRITE(ldtxt(ii),*) ' jpkglo = ', kpk ; ii = ii + 1100 !101 !102 80 ! ! Set the lateral boundary condition of the global domain 103 81 kperio = 0 ! ISOMIP configuration : close basin 104 82 ! 105 WRITE(ldtxt(ii),*) ' ' ; ii = ii + 1 106 WRITE(ldtxt(ii),*) ' Lateral boundary condition of the global domain' ; ii = ii + 1 107 WRITE(ldtxt(ii),*) ' ISOMIP : closed basin jperio = ', kperio ; ii = ii + 1 83 ! ! control print 84 IF(lwp) THEN 85 WRITE(numout,*) ' ' 86 WRITE(numout,*) 'usr_def_nam : read the user defined namelist (namusr_def) in namelist_cfg' 87 WRITE(numout,*) '~~~~~~~~~~~ ' 88 WRITE(numout,*) ' Namelist namusr_def : ISOMIP test case' 89 WRITE(numout,*) ' type of vertical coordinate : ' 90 WRITE(numout,*) ' z-coordinate flag ln_zco = ', ln_zco 91 WRITE(numout,*) ' z-partial-step coordinate flag ln_zps = ', ln_zps 92 WRITE(numout,*) ' s-coordinate flag ln_sco = ', ln_sco 93 WRITE(numout,*) ' resolution' 94 WRITE(numout,*) ' zonal resolution rn_e1deg = ', rn_e1deg, ' degrees' 95 WRITE(numout,*) ' meridional resolution rn_e1deg = ', rn_e1deg, ' degrees' 96 WRITE(numout,*) ' vertical resolution rn_e3 = ', rn_e3 , ' meters' 97 WRITE(numout,*) ' ISOMIP domain = 15° x 10° x 900 m' 98 WRITE(numout,*) ' resulting global domain size : jpiglo = ', kpi 99 WRITE(numout,*) ' jpjglo = ', kpj 100 WRITE(numout,*) ' jpkglo = ', kpk 101 WRITE(numout,*) ' ' 102 WRITE(numout,*) ' Lateral boundary condition of the global domain' 103 WRITE(numout,*) ' ISOMIP : closed basin jperio = ', kperio 104 ENDIF 108 105 ! 109 106 END SUBROUTINE usr_def_nam
Note: See TracChangeset
for help on using the changeset viewer.
