- Timestamp:
- 2020-09-14T17:40:34+02:00 (4 years ago)
- Location:
- NEMO/branches/2019/dev_r11351_fldread_with_XIOS
- Files:
-
- 3 edited
-
. (modified) (1 prop)
-
cfgs/C1D_PAPA/MY_SRC/usrdef_nam.F90 (modified) (4 diffs)
-
cfgs/C1D_PAPA/MY_SRC/usrdef_zgr.F90 (modified) (2 diffs)
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/dev_r11351_fldread_with_XIOS
- Property svn:externals
-
old new 3 3 ^/utils/build/mk@HEAD mk 4 4 ^/utils/tools@HEAD tools 5 ^/vendors/AGRIF/dev @HEADext/AGRIF5 ^/vendors/AGRIF/dev_r12970_AGRIF_CMEMS ext/AGRIF 6 6 ^/vendors/FCM@HEAD ext/FCM 7 7 ^/vendors/IOIPSL@HEAD ext/IOIPSL 8 9 # SETTE 10 ^/utils/CI/sette@13382 sette
-
- Property svn:externals
-
NEMO/branches/2019/dev_r11351_fldread_with_XIOS/cfgs/C1D_PAPA/MY_SRC/usrdef_nam.F90
r10072 r13463 39 39 CONTAINS 40 40 41 SUBROUTINE usr_def_nam( ldtxt, ldnam,cd_cfg, kk_cfg, kpi, kpj, kpk, kperio )41 SUBROUTINE usr_def_nam( cd_cfg, kk_cfg, kpi, kpj, kpk, kperio ) 42 42 !!---------------------------------------------------------------------- 43 43 !! *** ROUTINE dom_nam *** … … 51 51 !! ** input : - namusr_def namelist found in namelist_cfg 52 52 !!---------------------------------------------------------------------- 53 CHARACTER(len=*), DIMENSION(:), INTENT(out) :: ldtxt, ldnam ! stored print information54 53 CHARACTER(len=*) , INTENT(out) :: cd_cfg ! configuration name 55 54 INTEGER , INTENT(out) :: kk_cfg ! configuration resolution … … 57 56 INTEGER , INTENT(out) :: kperio ! lateral global domain b.c. 58 57 ! 59 INTEGER :: ios , ii! Local integer58 INTEGER :: ios ! Local integer 60 59 !! 61 60 NAMELIST/namusr_def/ rn_bathy 62 61 !!---------------------------------------------------------------------- 63 62 ! 64 ii = 1 63 READ ( numnam_cfg, namusr_def, IOSTAT = ios, ERR = 902 ) 64 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namusr_def in configuration namelist' ) 65 65 ! 66 REWIND( numnam_cfg ) ! Namelist namusr_def (exist in namelist_cfg only) 67 READ ( numnam_cfg, namusr_def, IOSTAT = ios, ERR = 902 ) 68 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namusr_def in configuration namelist', .TRUE. ) 69 ! 70 WRITE( ldnam(:), namusr_def ) 66 IF(lwm) WRITE( numond, namusr_def ) 71 67 ! 72 68 cd_cfg = 'C1D' ! name & resolution (not used) … … 77 73 kpj = 3 78 74 kpk = 75 79 !80 ! ! control print81 WRITE(ldtxt(ii),*) ' ' ; ii = ii + 182 WRITE(ldtxt(ii),*) 'usr_def_nam : read the user defined namelist (namusr_def) in namelist_cfg' ; ii = ii + 183 WRITE(ldtxt(ii),*) '~~~~~~~~~~~ ' ; ii = ii + 184 WRITE(ldtxt(ii),*) ' Namelist namusr_def : C1 case' ; ii = ii + 185 WRITE(ldtxt(ii),*) ' type of vertical coordinate : ' ; ii = ii + 186 WRITE(ldtxt(ii),*) ' z-coordinate flag ln_zco = ', ln_zco ; ii = ii + 187 WRITE(ldtxt(ii),*) ' z-partial-step coordinate flag ln_zps = ', ln_zps ; ii = ii + 188 WRITE(ldtxt(ii),*) ' s-coordinate flag ln_sco = ', ln_sco ; ii = ii + 189 WRITE(ldtxt(ii),*) ' C1D domain = 3 x 3 x75 grid-points ' ; ii = ii + 190 WRITE(ldtxt(ii),*) ' resulting global domain size : jpiglo = ', kpi ; ii = ii + 191 WRITE(ldtxt(ii),*) ' jpjglo = ', kpj ; ii = ii + 192 WRITE(ldtxt(ii),*) ' jpkglo = ', kpk ; ii = ii + 193 94 !95 75 ! ! Set the lateral boundary condition of the global domain 96 76 kperio = 7 ! C1D configuration : 3x3 basin with cyclic Est-West and Norht-South condition 97 77 ! 98 WRITE(ldtxt(ii),*) ' Lateral boundary condition of the global domain' ; ii = ii + 1 99 WRITE(ldtxt(ii),*) ' C1D : closed basin jperio = ', kperio ; ii = ii + 1 78 ! ! control print 79 IF(lwp) THEN 80 WRITE(numout,*) ' ' 81 WRITE(numout,*) 'usr_def_nam : read the user defined namelist (namusr_def) in namelist_cfg' 82 WRITE(numout,*) '~~~~~~~~~~~ ' 83 WRITE(numout,*) ' Namelist namusr_def : C1 case' 84 WRITE(numout,*) ' type of vertical coordinate : ' 85 WRITE(numout,*) ' z-coordinate flag ln_zco = ', ln_zco 86 WRITE(numout,*) ' z-partial-step coordinate flag ln_zps = ', ln_zps 87 WRITE(numout,*) ' s-coordinate flag ln_sco = ', ln_sco 88 WRITE(numout,*) ' C1D domain = 3 x 3 x75 grid-points ' 89 WRITE(numout,*) ' resulting global domain size : jpiglo = ', kpi 90 WRITE(numout,*) ' jpjglo = ', kpj 91 WRITE(numout,*) ' jpkglo = ', kpk 92 WRITE(numout,*) ' Lateral boundary condition of the global domain' 93 WRITE(numout,*) ' C1D : closed basin jperio = ', kperio 94 ENDIF 100 95 ! 101 96 END SUBROUTINE usr_def_nam -
NEMO/branches/2019/dev_r11351_fldread_with_XIOS/cfgs/C1D_PAPA/MY_SRC/usrdef_zgr.F90
r10072 r13463 31 31 32 32 !! * Substitutions 33 # include " vectopt_loop_substitute.h90"33 # include "do_loop_substitute.h90" 34 34 !!---------------------------------------------------------------------- 35 35 !! NEMO/OCE 4.0 , NEMO Consortium (2018) … … 159 159 pe3vw(:,:,jk) = pe3w_1d (jk) 160 160 END DO 161 DO jj = 1, jpj ! bottom scale factors and depth at T- and W-points 162 DO ji = 1, jpi 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 DO 172 END DO 161 ! bottom scale factors and depth at T- and W-points 162 DO_2D( 1, 1, 1, 1 ) 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 173 172 ! ! bottom scale factors and depth at U-, V-, UW and VW-points 174 173 ! ! usually Computed as the minimum of neighbooring scale factors
Note: See TracChangeset
for help on using the changeset viewer.
