Changeset 479 for trunk/NEMO/OPA_SRC
- Timestamp:
- 2006-06-19T15:19:23+02:00 (18 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/OPA_SRC/istate.F90
r473 r479 238 238 INTEGER :: inum ! temporary logical unit 239 239 INTEGER :: ji, jj, jk ! dummy loop indices 240 INTEGER :: ijloc 240 241 REAL(wp) :: & 241 zh1, zh2, zslope, zcst ! temporary scalars 242 zh1, zh2, zslope, zcst, & ! temporary scalars 243 zfcor 242 244 REAL(wp) :: & 243 245 zt1 = 12._wp, & ! surface temperature value (EEL R5) 244 246 zt2 = 2._wp, & ! bottom temperature value (EEL R5) 245 zsal = 35.5_wp ! constant salinity (EEL R2, R5 and R6) 247 zsal = 35.5_wp, & ! constant salinity (EEL R2, R5 and R6) 248 zueel = 0.1_wp ! constant uniform zonal velocity (EEL R5) 246 249 # if ! defined key_dynspg_rl 247 250 REAL(wp), DIMENSION(jpiglo,jpjglo) :: & … … 291 294 ! Start EEL5 configuration with barotropic geostrophic velocities 292 295 ! according the sshb and sshn SSH imposed. 293 ub(:,:,:) = 0.1 * umask(:,:,:) 296 ! we assume a uniform grid (hence the use of e1t(1,1) for delta_y) 297 ! we use the Coriolis frequency at mid-channel. 298 299 ub(:,:,:) = zueel * umask(:,:,:) 294 300 un(:,:,:) = ub(:,:,:) 301 ijloc = mj0(INT(jpjglo-1)/2) 302 zfcor = ff(1,ijloc) 295 303 296 304 DO jj = 1, jpjglo 297 zssh(:,jj) = ( .22 - ( FLOAT(jj-3) * (0.44) ) / 99. )305 zssh(:,jj) = - (FLOAT(jj)- FLOAT(jpjglo-1)/2.)*zueel*e1t(1,1)*zfcor/grav 298 306 END DO 307 308 IF(lwp) THEN 309 WRITE(numout,*) ' Uniform zonal velocity for EEL R5:',zueel 310 WRITE(numout,*) ' Geostrophic SSH profile as a function of y:' 311 WRITE(numout,'(12(1x,f6.2))') zssh(1,:) 312 ENDIF 313 299 314 DO jj = 1, nlcj 300 315 DO ji = 1, nlci
Note: See TracChangeset
for help on using the changeset viewer.