- Timestamp:
- 2020-03-23T22:16:19+01:00 (4 years ago)
- Location:
- NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/src/OCE/SBC
- Files:
-
- 2 edited
-
sbcrnf.F90 (modified) (7 diffs)
-
sbcssm.F90 (modified) (6 diffs)
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/src/OCE/SBC/sbcrnf.F90
r12377 r12590 34 34 PUBLIC sbc_rnf_alloc ! called in sbcmod module 35 35 PUBLIC sbc_rnf_init ! called in sbcmod module 36 36 37 37 ! !!* namsbc_rnf namelist * 38 38 CHARACTER(len=100) :: cn_dir !: Root directory for location of rnf files … … 58 58 LOGICAL , PUBLIC :: l_rnfcpl = .false. !: runoffs recieved from oasis 59 59 INTEGER , PUBLIC :: nkrnf = 0 !: nb of levels over which Kz is increased at river mouths 60 60 61 61 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: rnfmsk !: river mouth mask (hori.) 62 62 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) :: rnfmsk_z !: river mouth mask (vert.) 63 63 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: h_rnf !: depth of runoff in m 64 64 INTEGER, PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: nk_rnf !: depth of runoff in model levels 65 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: rnf_tsc_b, rnf_tsc !: before and now T & S runoff contents [K.m/s & PSU.m/s] 65 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: rnf_tsc_b, rnf_tsc !: before and now T & S runoff contents [K.m/s & PSU.m/s] 66 66 67 67 TYPE(FLD), ALLOCATABLE, DIMENSION(:) :: sf_rnf ! structure: river runoff (file information, fields read) 68 68 TYPE(FLD), ALLOCATABLE, DIMENSION(:) :: sf_i_rnf ! structure: iceberg flux (file information, fields read) 69 TYPE(FLD), ALLOCATABLE, DIMENSION(:) :: sf_s_rnf ! structure: river runoff salinity (file information, fields read) 70 TYPE(FLD), ALLOCATABLE, DIMENSION(:) :: sf_t_rnf ! structure: river runoff temperature (file information, fields read) 71 69 TYPE(FLD), ALLOCATABLE, DIMENSION(:) :: sf_s_rnf ! structure: river runoff salinity (file information, fields read) 70 TYPE(FLD), ALLOCATABLE, DIMENSION(:) :: sf_t_rnf ! structure: river runoff temperature (file information, fields read) 71 72 72 !! * Substitutions 73 73 # include "do_loop_substitute.h90" 74 # include "domzgr_substitute.h90" 74 75 !!---------------------------------------------------------------------- 75 76 !! NEMO/OCE 4.0 , NEMO Consortium (2018) … … 227 228 ELSE !== runoff put only at the surface ==! 228 229 h_rnf (:,:) = e3t (:,:,1,Kmm) ! update h_rnf to be depth of top box 229 phdivn(:,:,1) = phdivn(:,:,1) - ( rnf(:,:) +rnf_b(:,:) ) * zfact * r1_rau0 / e3t(:,:,1,Kmm)230 phdivn(:,:,1) = phdivn(:,:,1) - ( rnf(:,:)+rnf_b(:,:) ) * zfact * r1_rau0 / e3t(:,:,1,Kmm) 230 231 ENDIF 231 232 ! … … 249 250 INTEGER :: ios ! Local integer output status for namelist read 250 251 INTEGER :: nbrec ! temporary integer 251 REAL(wp) :: zacoef 252 REAL(wp), DIMENSION(jpi,jpj,2) :: zrnfcl 252 REAL(wp) :: zacoef 253 REAL(wp), DIMENSION(jpi,jpj,2) :: zrnfcl 253 254 !! 254 255 NAMELIST/namsbc_rnf/ cn_dir , ln_rnf_depth, ln_rnf_tem, ln_rnf_sal, ln_rnf_icb, & … … 261 262 IF( sbc_rnf_alloc() /= 0 ) CALL ctl_stop( 'STOP', 'sbc_rnf_alloc : unable to allocate arrays' ) 262 263 ! 263 IF( .NOT. ln_rnf ) THEN ! no specific treatment in vicinity of river mouths 264 IF( .NOT. ln_rnf ) THEN ! no specific treatment in vicinity of river mouths 264 265 ln_rnf_mouth = .FALSE. ! default definition needed for example by sbc_ssr or by tra_adv_muscl 265 266 nkrnf = 0 … … 297 298 ! ! ================== 298 299 ! 299 IF( .NOT. l_rnfcpl ) THEN 300 IF( .NOT. l_rnfcpl ) THEN 300 301 ALLOCATE( sf_rnf(1), STAT=ierror ) ! Create sf_rnf structure (runoff inflow) 301 302 IF(lwp) WRITE(numout,*) … … 352 353 IF(lwp) WRITE(numout,*) ' ==>>> runoffs depth read in a file' 353 354 rn_dep_file = TRIM( cn_dir )//TRIM( sn_dep_rnf%clname ) 354 IF( .NOT. sn_dep_rnf%ln_clim ) THEN ; WRITE(rn_dep_file, '(a,"_y",i4)' ) TRIM( rn_dep_file ), nyear ! add year 355 IF( sn_dep_rnf%cltype == 'monthly' ) WRITE(rn_dep_file, '(a,"m",i2)' ) TRIM( rn_dep_file ), nmonth ! add month 355 IF( .NOT. sn_dep_rnf%ln_clim ) THEN ; WRITE(rn_dep_file, '(a,"_y",i4)' ) TRIM( rn_dep_file ), nyear ! add year 356 IF( sn_dep_rnf%cltype == 'monthly' ) WRITE(rn_dep_file, '(a,"m",i2)' ) TRIM( rn_dep_file ), nmonth ! add month 356 357 ENDIF 357 358 CALL iom_open ( rn_dep_file, inum ) ! open file -
NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/src/OCE/SBC/sbcssm.F90
r12377 r12590 10 10 11 11 !!---------------------------------------------------------------------- 12 !! sbc_ssm : calculate sea surface mean currents, temperature, 12 !! sbc_ssm : calculate sea surface mean currents, temperature, 13 13 !! and salinity over nn_fsbc time-step 14 14 !!---------------------------------------------------------------------- … … 31 31 32 32 LOGICAL, SAVE :: l_ssm_mean = .FALSE. ! keep track of whether means have been read from restart file 33 33 34 # include "domzgr_substitute.h90" 34 35 !!---------------------------------------------------------------------- 35 36 !! NEMO/OCE 4.0 , NEMO Consortium (2018) … … 42 43 !!--------------------------------------------------------------------- 43 44 !! *** ROUTINE sbc_oce *** 44 !! 45 !! 45 46 !! ** Purpose : provide ocean surface variable to sea-surface boundary 46 !! condition computation 47 !! 48 !! ** Method : compute mean surface velocity (2 components at U and 47 !! condition computation 48 !! 49 !! ** Method : compute mean surface velocity (2 components at U and 49 50 !! V-points) [m/s], temperature [Celsius] and salinity [psu] over 50 51 !! the periode (kt - nn_fsbc) to kt … … 200 201 ! 201 202 ELSE 202 ! 203 ! 203 204 IF(lwp) WRITE(numout,*) 204 205 IF(lwp) WRITE(numout,*) 'sbc_ssm_init : sea surface mean fields' … … 222 223 ! 223 224 IF( zf_sbc /= REAL( nn_fsbc, wp ) ) THEN ! nn_fsbc has changed between 2 runs 224 IF(lwp) WRITE(numout,*) ' restart with a change in the frequency of mean from ', zf_sbc, ' to ', nn_fsbc 225 zcoef = REAL( nn_fsbc - 1, wp ) / zf_sbc 226 ssu_m(:,:) = zcoef * ssu_m(:,:) 225 IF(lwp) WRITE(numout,*) ' restart with a change in the frequency of mean from ', zf_sbc, ' to ', nn_fsbc 226 zcoef = REAL( nn_fsbc - 1, wp ) / zf_sbc 227 ssu_m(:,:) = zcoef * ssu_m(:,:) 227 228 ssv_m(:,:) = zcoef * ssv_m(:,:) 228 229 sst_m(:,:) = zcoef * sst_m(:,:) … … 252 253 ENDIF 253 254 ! 254 IF( .NOT. ln_traqsr ) fraqsr_1lev(:,:) = 1._wp ! default definition: qsr 100% in the fisrt level 255 IF( .NOT. ln_traqsr ) fraqsr_1lev(:,:) = 1._wp ! default definition: qsr 100% in the fisrt level 255 256 ! 256 257 IF( lwxios.AND.nn_fsbc > 1 ) THEN
Note: See TracChangeset
for help on using the changeset viewer.
