Changeset 13998 for NEMO/branches/2020/dev_r13327_KERNEL-06_2_techene_e3/tests/CANAL/MY_SRC/usrdef_istate.F90
- Timestamp:
- 2020-12-02T14:55:21+01:00 (3 years ago)
- Location:
- NEMO/branches/2020/dev_r13327_KERNEL-06_2_techene_e3
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2020/dev_r13327_KERNEL-06_2_techene_e3
- Property svn:externals
-
old new 8 8 9 9 # SETTE 10 ^/utils/CI/sette@13 292sette10 ^/utils/CI/sette@13559 sette
-
- Property svn:externals
-
NEMO/branches/2020/dev_r13327_KERNEL-06_2_techene_e3/tests/CANAL/MY_SRC/usrdef_istate.F90
r13874 r13998 66 66 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~~~ ' 67 67 ! 68 IF (ln_sshnoise) CALL RANDOM_NUMBER(zrandom)69 68 zjetx = ABS(rn_ujetszx)/2. 70 69 zjety = ABS(rn_ujetszy)/2. 71 70 ! 71 zf0 = 2._wp * omega * SIN( rad * rn_ppgphi0 ) 72 ! 72 73 SELECT CASE(nn_initcase) 74 75 CASE(-1) ! stratif at rest 76 77 ! sea level: 78 pssh(:,:) = 0. 79 ! temperature: 80 pts(:,:,1,jp_tem) = 25. !!30._wp 81 pts(:,:,2:jpk,jp_tem) = 22. !!24._wp 82 ! salinity: 83 pts(:,:,:,jp_sal) = 35._wp 84 ! velocities: 85 pu(:,:,:) = 0. 86 pv(:,:,:) = 0. 87 73 88 CASE(0) ! rest 74 89 ! … … 88 103 pts(:,:,jpk,jp_sal) = 0. 89 104 DO jk=1, jpkm1 90 pts(:,:,jk,jp_sal) = gphit(:,:) 105 WHERE( ABS(gphit) <= zjety ) 106 !!$ WHERE( ABS(gphit) <= zjety*0.5 .AND. ABS(glamt) <= zjety*0.5 ) ! for a square of salt 107 pts(:,:,jk,jp_sal) = 35. 108 ELSEWHERE 109 pts(:,:,jk,jp_sal) = 30. 110 END WHERE 91 111 END DO 92 112 ! velocities: … … 102 122 pts(:,:,:,jp_tem) = 10._wp 103 123 ! salinity: 104 pts(:,:,:,jp_sal) = 2.105 DO jk=1, jpkm1 106 WHERE( ABS(gphiv) <= zjety ) pts(:,:,jk,jp_sal) = 2. + SIGN(1.,gphiv(:,:))124 pts(:,:,:,jp_sal) = 30. 125 DO jk=1, jpkm1 126 WHERE( ABS(gphiv) <= zjety ) pts(:,:,jk,jp_sal) = 30. + SIGN(1.,gphiv(:,:)) 107 127 END DO 108 128 ! velocities: … … 125 145 ! salinity: 126 146 DO jk=1, jpkm1 127 pts(:,:,jk,jp_sal) = gphit(:,:)147 pts(:,:,jk,jp_sal) = pssh(:,:) 128 148 END DO 129 149 ! velocities: … … 157 177 zf0 = 2._wp * omega * SIN( rad * rn_ppgphi0 ) 158 178 zumax = rn_vtxmax * SIGN(1._wp, zf0) ! Here Anticyclonic: set zumax=-1 for cyclonic 159 zlambda = SQRT(2._wp)*rn_lambda ! Horizontal scale in meters179 zlambda = SQRT(2._wp)*rn_lambda*1.e3 ! Horizontal scale in meters 160 180 zn2 = 3.e-3**2 161 181 zH = 0.5_wp * 5000._wp … … 352 372 ! !== add noise ==! 353 373 IF (ln_sshnoise) THEN 374 CALL RANDOM_SEED() 354 375 CALL RANDOM_NUMBER(zrandom) 355 376 pssh(:,:) = pssh(:,:) + ( 0.1 * zrandom(:,:) - 0.05 )
Note: See TracChangeset
for help on using the changeset viewer.