Changeset 12489 for NEMO/trunk/tests/VORTEX/MY_SRC/usrdef_istate.F90
- Timestamp:
- 2020-02-28T16:55:11+01:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk/tests/VORTEX/MY_SRC/usrdef_istate.F90
r10425 r12489 69 69 zH = 0.5_wp * 5000._wp 70 70 ! 71 zP0 = r au0 * zf0 * zumax * zlambda * SQRT(EXP(1._wp)/2._wp)71 zP0 = rho0 * zf0 * zumax * zlambda * SQRT(EXP(1._wp)/2._wp) 72 72 ! 73 73 ! Sea level: … … 77 77 zx = glamt(ji,jj) * 1.e3 78 78 zy = gphit(ji,jj) * 1.e3 79 zrho1 = r au0 + za * EXP(-(zx**2+zy**2)/zlambda**2)79 zrho1 = rho0 + za * EXP(-(zx**2+zy**2)/zlambda**2) 80 80 pssh(ji,jj) = zP0 * EXP(-(zx**2+zy**2)/zlambda**2)/(zrho1*grav) * ptmask(ji,jj,1) 81 81 END DO … … 89 89 DO jk=1,jpk 90 90 zdt = pdept(ji,jj,jk) 91 zrho1 = r au0 * (1._wp + zn2*zdt/grav)91 zrho1 = rho0 * (1._wp + zn2*zdt/grav) 92 92 IF (zdt < zH) THEN 93 93 zrho1 = zrho1 - zP0 * (1._wp-EXP(zdt-zH)) & 94 94 & * EXP(-(zx**2+zy**2)/zlambda**2) / (grav*(zH -1._wp + exp(-zH))); 95 95 ENDIF 96 pts(ji,jj,jk,jp_tem) = (20._wp + (r au0-zrho1) / 0.28_wp) * ptmask(ji,jj,jk)96 pts(ji,jj,jk,jp_tem) = (20._wp + (rho0-zrho1) / 0.28_wp) * ptmask(ji,jj,jk) 97 97 END DO 98 98 END DO … … 103 103 ! 104 104 ! velocities: 105 za = 2._wp * zP0 / (zf0 * r au0 * zlambda**2)105 za = 2._wp * zP0 / (zf0 * rho0 * zlambda**2) 106 106 DO ji=1, jpim1 107 107 DO jj=1, jpj
Note: See TracChangeset
for help on using the changeset viewer.