[888] | 1 | MODULE sbcrnf |
---|
| 2 | !!====================================================================== |
---|
| 3 | !! *** MODULE sbcrnf *** |
---|
| 4 | !! Ocean forcing: river runoff |
---|
| 5 | !!===================================================================== |
---|
[1482] | 6 | !! History : OPA ! 2000-11 (R. Hordoir, E. Durand) NetCDF FORMAT |
---|
| 7 | !! NEMO 1.0 ! 2002-09 (G. Madec) F90: Free form and module |
---|
| 8 | !! 3.0 ! 2006-07 (G. Madec) Surface module |
---|
| 9 | !! 3.2 ! 2009-04 (B. Lemaire) Introduce iom_put |
---|
[888] | 10 | !!---------------------------------------------------------------------- |
---|
| 11 | |
---|
| 12 | !!---------------------------------------------------------------------- |
---|
| 13 | !! sbc_rnf : monthly runoffs read in a NetCDF file |
---|
[1116] | 14 | !! sbc_rnf_init : runoffs initialisation |
---|
[888] | 15 | !! rnf_mouth : set river mouth mask |
---|
| 16 | !!---------------------------------------------------------------------- |
---|
| 17 | USE dom_oce ! ocean space and time domain |
---|
| 18 | USE phycst ! physical constants |
---|
| 19 | USE sbc_oce ! surface boundary condition variables |
---|
| 20 | USE fldread ! ??? |
---|
| 21 | USE in_out_manager ! I/O manager |
---|
| 22 | USE iom ! I/O module |
---|
| 23 | |
---|
| 24 | IMPLICIT NONE |
---|
| 25 | PRIVATE |
---|
| 26 | |
---|
| 27 | PUBLIC sbc_rnf ! routine call in step module |
---|
| 28 | |
---|
[1482] | 29 | ! !!* namsbc_rnf namelist * |
---|
[1116] | 30 | CHARACTER(len=100), PUBLIC :: cn_dir = './' !: Root directory for location of ssr files |
---|
| 31 | LOGICAL , PUBLIC :: ln_rnf_emp = .false. !: runoffs into a file to be read or already into precipitation |
---|
| 32 | TYPE(FLD_N) , PUBLIC :: sn_rnf !: information about the runoff file to be read |
---|
| 33 | TYPE(FLD_N) , PUBLIC :: sn_cnf !: information about the runoff mouth file to be read |
---|
[1938] | 34 | TYPE(FLD_N) :: sn_sal_rnf !: information about the salinities of runoff file to be read |
---|
| 35 | TYPE(FLD_N) :: sn_tmp_rnf !: information about the temperatures of runoff file to be read |
---|
| 36 | TYPE(FLD_N) :: sn_dep_rnf !: information about the depth which river inflow affects |
---|
[1116] | 37 | LOGICAL , PUBLIC :: ln_rnf_mouth = .false. !: specific treatment in mouths vicinity |
---|
| 38 | REAL(wp) , PUBLIC :: rn_hrnf = 0.e0 !: runoffs, depth over which enhanced vertical mixing is used |
---|
| 39 | REAL(wp) , PUBLIC :: rn_avt_rnf = 0.e0 !: runoffs, value of the additional vertical mixing coef. [m2/s] |
---|
[1998] | 40 | LOGICAL , PUBLIC :: ln_rnf_att = .false. !: river runoffs attributes (temp, sal & depth) are specified in a file |
---|
[1601] | 41 | REAL(wp) , PUBLIC :: rn_rfact = 1.e0 !: multiplicative factor for runoff |
---|
[888] | 42 | |
---|
| 43 | INTEGER , PUBLIC :: nkrnf = 0 !: number of levels over which Kz is increased at river mouths |
---|
| 44 | REAL(wp), PUBLIC, DIMENSION(jpi,jpj) :: rnfmsk !: river mouth mask (hori.) |
---|
| 45 | REAL(wp), PUBLIC, DIMENSION(jpk) :: rnfmsk_z !: river mouth mask (vert.) |
---|
| 46 | |
---|
[1998] | 47 | TYPE(FLD), ALLOCATABLE, DIMENSION(:) :: sf_rnf !: structure of input river runoff (file information, fields read) |
---|
[888] | 48 | |
---|
[1998] | 49 | TYPE(FLD), ALLOCATABLE, DIMENSION(:) :: sf_sal_rnf !: structure of input river runoff salinity (file information, fields read) |
---|
| 50 | TYPE(FLD), ALLOCATABLE, DIMENSION(:) :: sf_tmp_rnf !: structure of input river runoff temperature (file information, fields read) |
---|
[1938] | 51 | |
---|
| 52 | REAL, PUBLIC, DIMENSION(jpi,jpj) :: rnf_dep !: depth of runoff in m |
---|
| 53 | INTEGER, PUBLIC, DIMENSION(jpi,jpj) :: rnf_mod_dep !: depth of runoff in model levels |
---|
| 54 | REAL, PUBLIC, DIMENSION(jpi,jpj) :: rnf_sal !: salinity of river runoff |
---|
| 55 | REAL, PUBLIC, DIMENSION(jpi,jpj) :: rnf_tmp !: temperature of river runoff |
---|
| 56 | |
---|
| 57 | INTEGER :: ji, jj ,jk ! dummy loop indices |
---|
| 58 | INTEGER :: inum ! temporary logical unit |
---|
| 59 | |
---|
| 60 | !! * Substitutions |
---|
| 61 | # include "domzgr_substitute.h90" |
---|
| 62 | |
---|
[888] | 63 | !!---------------------------------------------------------------------- |
---|
[1482] | 64 | !! NEMO/OPA 3.2 , LOCEAN-IPSL (2009) |
---|
[1146] | 65 | !! $Id$ |
---|
[888] | 66 | !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt) |
---|
| 67 | !!---------------------------------------------------------------------- |
---|
| 68 | |
---|
| 69 | CONTAINS |
---|
| 70 | |
---|
| 71 | SUBROUTINE sbc_rnf( kt ) |
---|
| 72 | !!---------------------------------------------------------------------- |
---|
| 73 | !! *** ROUTINE sbc_rnf *** |
---|
| 74 | !! |
---|
| 75 | !! ** Purpose : Introduce a climatological run off forcing |
---|
| 76 | !! |
---|
| 77 | !! ** Method : Set each river mouth with a monthly climatology |
---|
| 78 | !! provided from different data. |
---|
| 79 | !! CAUTION : upward water flux, runoff forced to be < 0 |
---|
| 80 | !! |
---|
| 81 | !! ** Action : runoff updated runoff field at time-step kt |
---|
| 82 | !!---------------------------------------------------------------------- |
---|
| 83 | INTEGER, INTENT(in) :: kt ! ocean time step |
---|
[1482] | 84 | !! |
---|
[888] | 85 | INTEGER :: ji, jj ! dummy loop indices |
---|
[1242] | 86 | INTEGER :: ierror ! temporary integer |
---|
[888] | 87 | !!---------------------------------------------------------------------- |
---|
[1116] | 88 | ! |
---|
[1242] | 89 | IF( kt == nit000 ) THEN |
---|
[1998] | 90 | CALL sbc_rnf_init ! Read namelist and allocate structures |
---|
[1242] | 91 | ENDIF |
---|
[888] | 92 | |
---|
| 93 | ! !-------------------! |
---|
[1116] | 94 | IF( .NOT. ln_rnf_emp ) THEN ! Update runoff ! |
---|
[888] | 95 | ! !-------------------! |
---|
| 96 | ! |
---|
| 97 | CALL fld_read( kt, nn_fsbc, sf_rnf ) ! Read Runoffs data and provides it |
---|
| 98 | ! ! at the current time-step |
---|
[1938] | 99 | IF ( ln_rnf_att ) THEN |
---|
| 100 | CALL fld_read ( kt, nn_fsbc, sf_sal_rnf ) |
---|
| 101 | CALL fld_read ( kt, nn_fsbc, sf_tmp_rnf ) |
---|
| 102 | ENDIF |
---|
[888] | 103 | |
---|
| 104 | ! Runoff reduction only associated to the ORCA2_LIM configuration |
---|
| 105 | ! when reading the NetCDF file runoff_1m_nomask.nc |
---|
| 106 | IF( cp_cfg == 'orca' .AND. jp_cfg == 2 ) THEN |
---|
| 107 | DO jj = 1, jpj |
---|
| 108 | DO ji = 1, jpi |
---|
| 109 | IF( gphit(ji,jj) > 40 .AND. gphit(ji,jj) < 65 ) sf_rnf(1)%fnow(ji,jj) = 0.85 * sf_rnf(1)%fnow(ji,jj) |
---|
| 110 | END DO |
---|
| 111 | END DO |
---|
| 112 | ENDIF |
---|
| 113 | |
---|
| 114 | ! C a u t i o n : runoff is negative and in kg/m2/s |
---|
| 115 | |
---|
[1303] | 116 | IF( MOD( kt-1, nn_fsbc ) == 0 ) THEN |
---|
[1938] | 117 | rnf(:,:) = rn_rfact * ( sf_rnf(1)%fnow(:,:) ) |
---|
| 118 | IF ( ln_rnf_att ) THEN |
---|
| 119 | rnf_sal(:,:) = ( sf_sal_rnf(1)%fnow(:,:) ) |
---|
| 120 | rnf_tmp(:,:) = ( sf_tmp_rnf(1)%fnow(:,:) ) |
---|
| 121 | ELSE |
---|
| 122 | rnf_sal(:,:) = 0 |
---|
| 123 | rnf_tmp(:,:) = -999 |
---|
| 124 | ENDIF |
---|
[1998] | 125 | CALL iom_put( "runoffs", rnf ) ! runoffs |
---|
[1303] | 126 | ENDIF |
---|
[888] | 127 | ! |
---|
| 128 | ENDIF |
---|
| 129 | ! |
---|
| 130 | END SUBROUTINE sbc_rnf |
---|
| 131 | |
---|
| 132 | |
---|
[1998] | 133 | SUBROUTINE sbc_rnf_init |
---|
[1116] | 134 | !!---------------------------------------------------------------------- |
---|
| 135 | !! *** ROUTINE sbc_rnf_init *** |
---|
| 136 | !! |
---|
| 137 | !! ** Purpose : Initialisation of the runoffs if (ln_rnf=T) |
---|
| 138 | !! |
---|
| 139 | !! ** Method : - read the runoff namsbc_rnf namelist |
---|
| 140 | !! |
---|
| 141 | !! ** Action : - read parameters |
---|
| 142 | !!---------------------------------------------------------------------- |
---|
[1938] | 143 | CHARACTER(len=32) :: rn_dep_file ! runoff file name |
---|
| 144 | !! |
---|
| 145 | NAMELIST/namsbc_rnf/ cn_dir, ln_rnf_emp, sn_rnf, sn_cnf, sn_sal_rnf, sn_tmp_rnf, sn_dep_rnf, & |
---|
| 146 | & ln_rnf_mouth, ln_rnf_att, rn_hrnf, rn_avt_rnf, rn_rfact |
---|
[1116] | 147 | !!---------------------------------------------------------------------- |
---|
| 148 | |
---|
| 149 | ! ! ============ |
---|
| 150 | ! ! Namelist |
---|
| 151 | ! ! ============ |
---|
[1133] | 152 | ! (NB: frequency positive => hours, negative => months) |
---|
[1275] | 153 | ! ! file ! frequency ! variable ! time intep ! clim ! 'yearly' or ! weights ! rotation ! |
---|
| 154 | ! ! name ! (hours) ! name ! (T/F) ! (T/F) ! 'monthly' ! filename ! pairs ! |
---|
[1730] | 155 | sn_rnf = FLD_N( 'runoffs', -1 , 'sorunoff' , .TRUE. , .true. , 'yearly' , '' , '' ) |
---|
| 156 | sn_cnf = FLD_N( 'runoffs', 0 , 'sorunoff' , .FALSE. , .true. , 'yearly' , '' , '' ) |
---|
[1116] | 157 | |
---|
[1938] | 158 | sn_sal_rnf = FLD_N( 'runoffs', 24. , 'rosaline' , .TRUE. , .true. , 'yearly' , '' , '' ) |
---|
| 159 | sn_tmp_rnf = FLD_N( 'runoffs', 24. , 'rotemper' , .TRUE. , .true. , 'yearly' , '' , '' ) |
---|
| 160 | sn_dep_rnf = FLD_N( 'runoffs', 0. , 'rodepth' , .FALSE. , .true. , 'yearly' , '' , '' ) |
---|
[1116] | 161 | ! |
---|
| 162 | REWIND ( numnam ) ! Read Namelist namsbc_rnf |
---|
| 163 | READ ( numnam, namsbc_rnf ) |
---|
| 164 | |
---|
| 165 | ! ! Control print |
---|
| 166 | IF(lwp) THEN |
---|
| 167 | WRITE(numout,*) |
---|
| 168 | WRITE(numout,*) 'sbc_rnf : runoff ' |
---|
| 169 | WRITE(numout,*) '~~~~~~~ ' |
---|
| 170 | WRITE(numout,*) ' Namelist namsbc_rnf' |
---|
| 171 | WRITE(numout,*) ' runoff in a file to be read ln_rnf_emp = ', ln_rnf_emp |
---|
| 172 | WRITE(numout,*) ' specific river mouths treatment ln_rnf_mouth = ', ln_rnf_mouth |
---|
| 173 | WRITE(numout,*) ' river mouth additional Kz rn_avt_rnf = ', rn_avt_rnf |
---|
| 174 | WRITE(numout,*) ' depth of river mouth additional mixing rn_hrnf = ', rn_hrnf |
---|
[1601] | 175 | WRITE(numout,*) ' multiplicative factor for runoff rn_rfact = ', rn_rfact |
---|
[1116] | 176 | ENDIF |
---|
| 177 | |
---|
| 178 | ! ! ================== |
---|
| 179 | ! ! Type of runoff |
---|
| 180 | ! ! ================== |
---|
| 181 | ! |
---|
| 182 | IF( ln_rnf_emp ) THEN ! runoffs directly provided in the precipitations |
---|
| 183 | IF(lwp) WRITE(numout,*) |
---|
| 184 | IF(lwp) WRITE(numout,*) ' runoffs directly provided in the precipitations' |
---|
[1998] | 185 | IF ( ln_rnf_att ) THEN |
---|
| 186 | CALL ctl_warn( 'runoffs already included in precipitations & so runoff attributes will not be used' ) |
---|
| 187 | ln_rnf_att = .FALSE. |
---|
| 188 | ENDIF |
---|
[1116] | 189 | ! |
---|
| 190 | ELSE ! runoffs read in a file : set sf_rnf structure |
---|
| 191 | ! |
---|
[1998] | 192 | ! Allocate sf_rnf structure and (if required) sf_sal_rnf and sf_tmp_rnf structures |
---|
| 193 | ALLOCATE( sf_rnf(1), STAT=ierror ) |
---|
| 194 | IF( ierror > 0 ) THEN |
---|
| 195 | CALL ctl_stop( 'sbc_rnf: unable to allocate sf_rnf structure' ) ; RETURN |
---|
| 196 | ENDIF |
---|
| 197 | ALLOCATE( sf_rnf(1)%fnow(jpi,jpj) ) |
---|
| 198 | ALLOCATE( sf_rnf(1)%fdta(jpi,jpj,2) ) |
---|
| 199 | |
---|
| 200 | IF( ln_rnf_att ) THEN |
---|
| 201 | ALLOCATE( sf_sal_rnf(1), STAT=ierror ) |
---|
| 202 | IF( ierror > 0 ) THEN |
---|
| 203 | CALL ctl_stop( 'sbc_sal_rnf: unable to allocate sf_sal_rnf structure' ) ; RETURN |
---|
| 204 | ENDIF |
---|
| 205 | ALLOCATE( sf_sal_rnf(1)%fnow(jpi,jpj) ) |
---|
| 206 | ALLOCATE( sf_sal_rnf(1)%fdta(jpi,jpj,2) ) |
---|
| 207 | |
---|
| 208 | ALLOCATE( sf_tmp_rnf(1), STAT=ierror ) |
---|
| 209 | IF( ierror > 0 ) THEN |
---|
| 210 | CALL ctl_stop( 'sbc_tmp_rnf: unable to allocate sf_tmp_rnf structure' ) ; RETURN |
---|
| 211 | ENDIF |
---|
| 212 | ALLOCATE( sf_tmp_rnf(1)%fnow(jpi,jpj) ) |
---|
| 213 | ALLOCATE( sf_tmp_rnf(1)%fdta(jpi,jpj,2) ) |
---|
| 214 | ENDIF |
---|
[1133] | 215 | ! fill sf_rnf with sn_rnf and control print |
---|
| 216 | CALL fld_fill( sf_rnf, (/ sn_rnf /), cn_dir, 'sbc_rnf_init', 'read runoffs data', 'namsbc_rnf' ) |
---|
[1938] | 217 | |
---|
| 218 | IF ( ln_rnf_att ) THEN |
---|
| 219 | CALL fld_fill (sf_sal_rnf, (/ sn_sal_rnf /), cn_dir, 'sbc_rnf_init', 'read runoff salinity data', 'namsbc_rnf' ) |
---|
| 220 | CALL fld_fill (sf_tmp_rnf, (/ sn_tmp_rnf /), cn_dir, 'sbc_rnf_init', 'read runoff temperature data', 'namsbc_rnf' ) |
---|
| 221 | |
---|
| 222 | rn_dep_file = TRIM( cn_dir )//TRIM( sn_dep_rnf%clname ) |
---|
| 223 | CALL iom_open ( rn_dep_file, inum ) ! open file |
---|
| 224 | CALL iom_get ( inum, jpdom_data, sn_dep_rnf%clvar, rnf_dep ) ! read the river mouth array |
---|
| 225 | CALL iom_close( inum ) ! close file |
---|
| 226 | |
---|
| 227 | rnf_mod_dep(:,:)=0 |
---|
| 228 | DO jj=1,jpj |
---|
| 229 | DO ji=1,jpi |
---|
| 230 | IF ( rnf_dep(ji,jj) > 0.e0 ) THEN |
---|
| 231 | jk=2 |
---|
[2037] | 232 | DO WHILE ( jk/=(mbathy(ji,jj)-1) .AND. fsdept(ji,jj,jk) < rnf_dep(ji,jj) ); jk=jk+1; ENDDO |
---|
[1938] | 233 | rnf_mod_dep(ji,jj)=jk |
---|
| 234 | ELSE IF ( rnf_dep(ji,jj) .eq. -1 ) THEN |
---|
| 235 | rnf_mod_dep(ji,jj)=1 |
---|
| 236 | ELSE IF ( rnf_dep(ji,jj) .eq. -999 ) THEN |
---|
[2037] | 237 | rnf_mod_dep(ji,jj)=mbathy(ji,jj)-1 |
---|
[1938] | 238 | ELSE IF ( rnf_dep(ji,jj) /= 0 ) THEN |
---|
| 239 | CALL ctl_stop( 'runoff depth not positive, and not -999 or -1, rnf value in file fort.999' ) |
---|
| 240 | WRITE(999,*) 'ji, jj, rnf(ji,jj) :', ji, jj, rnf(ji,jj) |
---|
| 241 | ENDIF |
---|
| 242 | ENDDO |
---|
| 243 | ENDDO |
---|
| 244 | ELSE |
---|
| 245 | rnf_mod_dep(:,:)=1 |
---|
| 246 | ENDIF |
---|
| 247 | ! |
---|
[1116] | 248 | ENDIF |
---|
| 249 | |
---|
| 250 | ! ! ======================== |
---|
| 251 | ! ! River mouth vicinity |
---|
| 252 | ! ! ======================== |
---|
| 253 | ! |
---|
| 254 | IF( ln_rnf_mouth ) THEN ! Specific treatment in vicinity of river mouths : |
---|
| 255 | ! ! - Increase Kz in surface layers ( rn_hrnf > 0 ) |
---|
| 256 | ! ! - set to zero SSS damping (ln_ssr=T) |
---|
| 257 | ! ! - mixed upstream-centered (ln_traadv_cen2=T) |
---|
| 258 | ! |
---|
| 259 | ! ! Number of level over which Kz increase |
---|
[1938] | 260 | IF ( ln_rnf_att ) & |
---|
| 261 | & CALL ctl_warn( 'increased mixing turned on but effects may already be spread through depth by ln_rnf_att' ) |
---|
[1116] | 262 | nkrnf = 0 |
---|
| 263 | IF( rn_hrnf > 0.e0 ) THEN |
---|
| 264 | nkrnf = 2 |
---|
[1164] | 265 | DO WHILE( nkrnf /= jpkm1 .AND. gdepw_0(nkrnf+1) < rn_hrnf ) ; nkrnf = nkrnf + 1 ; END DO |
---|
[1116] | 266 | IF( ln_sco ) & |
---|
| 267 | CALL ctl_warn( 'sbc_rnf: number of levels over which Kz is increased is computed for zco...' ) |
---|
| 268 | ENDIF |
---|
| 269 | IF(lwp) WRITE(numout,*) |
---|
| 270 | IF(lwp) WRITE(numout,*) ' Specific treatment used in vicinity of river mouths :' |
---|
| 271 | IF(lwp) WRITE(numout,*) ' - Increase Kz in surface layers (if rn_hrnf > 0 )' |
---|
| 272 | IF(lwp) WRITE(numout,*) ' by ', rn_avt_rnf,' m2/s over ', nkrnf, ' w-levels' |
---|
| 273 | IF(lwp) WRITE(numout,*) ' - set to zero SSS damping (if ln_ssr=T)' |
---|
| 274 | IF(lwp) WRITE(numout,*) ' - mixed upstream-centered (if ln_traadv_cen2=T)' |
---|
| 275 | ! |
---|
| 276 | CALL rnf_mouth ! set river mouth mask |
---|
| 277 | ! |
---|
| 278 | ELSE ! No treatment at river mouths |
---|
| 279 | IF(lwp) WRITE(numout,*) |
---|
| 280 | IF(lwp) WRITE(numout,*) ' No specific treatment at river mouths' |
---|
| 281 | rnfmsk (:,:) = 0.e0 |
---|
| 282 | rnfmsk_z(:) = 0.e0 |
---|
| 283 | nkrnf = 0 |
---|
| 284 | ENDIF |
---|
| 285 | |
---|
| 286 | END SUBROUTINE sbc_rnf_init |
---|
| 287 | |
---|
| 288 | |
---|
[888] | 289 | SUBROUTINE rnf_mouth |
---|
| 290 | !!---------------------------------------------------------------------- |
---|
| 291 | !! *** ROUTINE rnf_mouth *** |
---|
| 292 | !! |
---|
| 293 | !! ** Purpose : define the river mouths mask |
---|
| 294 | !! |
---|
| 295 | !! ** Method : read the river mouth mask (=0/1) in the river runoff |
---|
| 296 | !! climatological file. Defined a given vertical structure. |
---|
| 297 | !! CAUTION, the vertical structure is hard coded on the |
---|
| 298 | !! first 5 levels. |
---|
| 299 | !! This fields can be used to: |
---|
| 300 | !! - set an upstream advection scheme |
---|
[1116] | 301 | !! (ln_rnf_mouth=T and ln_traadv_cen2=T) |
---|
[888] | 302 | !! - increase vertical on the top nn_krnf vertical levels |
---|
| 303 | !! at river runoff input grid point (nn_krnf>=2, see step.F90) |
---|
| 304 | !! - set to zero SSS restoring flux at river mouth grid points |
---|
| 305 | !! |
---|
| 306 | !! ** Action : rnfmsk set to 1 at river runoff input, 0 elsewhere |
---|
| 307 | !! rnfmsk_z vertical structure |
---|
| 308 | !!---------------------------------------------------------------------- |
---|
[1601] | 309 | USE closea, ONLY : clo_rnf ! rnfmsk update routine |
---|
[888] | 310 | ! |
---|
| 311 | INTEGER :: inum ! temporary integers |
---|
| 312 | CHARACTER(len=32) :: cl_rnfile ! runoff file name |
---|
| 313 | !!---------------------------------------------------------------------- |
---|
| 314 | ! |
---|
| 315 | IF(lwp) WRITE(numout,*) |
---|
| 316 | IF(lwp) WRITE(numout,*) 'rnf_mouth : river mouth mask' |
---|
| 317 | IF(lwp) WRITE(numout,*) '~~~~~~~~~ ' |
---|
| 318 | |
---|
[1133] | 319 | cl_rnfile = TRIM( cn_dir )//TRIM( sn_cnf%clname ) |
---|
| 320 | IF( .NOT. sn_cnf%ln_clim ) THEN ; WRITE(cl_rnfile, '(a,"_y",i4)' ) TRIM( cl_rnfile ), nyear ! add year |
---|
| 321 | IF( sn_cnf%cltype == 'monthly' ) WRITE(cl_rnfile, '(a,"m",i2)' ) TRIM( cl_rnfile ), nmonth ! add month |
---|
| 322 | ENDIF |
---|
| 323 | |
---|
[888] | 324 | ! horizontal mask (read in NetCDF file) |
---|
| 325 | CALL iom_open ( cl_rnfile, inum ) ! open file |
---|
| 326 | CALL iom_get ( inum, jpdom_data, sn_cnf%clvar, rnfmsk ) ! read the river mouth array |
---|
| 327 | CALL iom_close( inum ) ! close file |
---|
| 328 | |
---|
| 329 | IF( nclosea == 1 ) CALL clo_rnf( rnfmsk ) ! closed sea inflow set as ruver mouth |
---|
| 330 | |
---|
| 331 | rnfmsk_z(:) = 0.e0 ! vertical structure |
---|
| 332 | rnfmsk_z(1) = 1.0 |
---|
| 333 | rnfmsk_z(2) = 1.0 ! ********** |
---|
| 334 | rnfmsk_z(3) = 0.5 ! HARD CODED on the 5 first levels |
---|
| 335 | rnfmsk_z(4) = 0.25 ! ********** |
---|
| 336 | rnfmsk_z(5) = 0.125 |
---|
| 337 | ! |
---|
| 338 | END SUBROUTINE rnf_mouth |
---|
| 339 | |
---|
| 340 | !!====================================================================== |
---|
| 341 | END MODULE sbcrnf |
---|