Changeset 2128 for branches/devukmo2010/NEMO/OPA_SRC/LDF/ldfdyn_c1d.h90
- Timestamp:
- 2010-09-28T14:29:51+02:00 (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/devukmo2010/NEMO/OPA_SRC/LDF/ldfdyn_c1d.h90
r1152 r2128 17 17 !! ** Method : 1D eddy viscosity coefficients ( depth ) 18 18 !! ahm3, ahm4 never used 19 !! biharmonic or harmonic operator : ahm1=ahm2 defined at T-level 19 !! harmonic operator : ahm1 defined at T-level 20 !! biharmonic operator : ahm2 defined at T-level 20 21 !! isopycnal or geopotential harmonic operator 21 22 !! : ahm1 defined at T-level … … 27 28 28 29 !! * Local variables 29 REAL(wp) :: zdam, zwam, zm00, zm01, zmhf, zmhs 30 REAL(wp) :: zdam, zwam, zm00, zm01, zmhf, zmhs 31 REAL(wp) :: zdam2, zwam2, zm200, zm201, zmh2f, zmh2s 30 32 REAL(wp) :: zahmf, zahms 31 33 !!---------------------------------------------------------------------- … … 36 38 IF(lwp) WRITE(numout,*) 37 39 38 ! Set ahm1 =ahm2(always at t-level)40 ! Set ahm1 for laplacian (always at t-level) 39 41 ! ============= 40 42 ! (USER: modify ahm1 following your desiderata) … … 53 55 zmhs = zahms-zmhf * zm00 54 56 55 ! set ahm1=ahm2 at T-level 57 ! Set ahm2 for bilaplacian (always at t-level) 58 ! ============= 59 ! (USER: modify ahm2 following your desiderata) 60 61 ! initialization of the profile 62 ! ahms, ahmf: surface and bottom values 63 zahm2s = ahm0_blp 64 zahm2f = ahm0_blp/4. 65 ! zdam, zwam: depth of the inflection pt and width of inflection 66 zdam2 = -300. 67 zwam2 = 300. 68 ! computation coefficients 69 zm200 = TANH( (0-zdam2)/zwam2 ) 70 zm201 = TANH( (-fsdept(1,1,jpk)-zdam2)/zwam2 ) 71 zmh2f = (zahm2s-zahm2f)/(zm200-zm201) 72 zmh2s = zahm2s-zmh2f * zm00 73 74 75 ! set ahm1 and ahm2 at T-level 56 76 DO jk = 1, jpk 57 ahm1(jk) = zmhs + zmhf * TANH( (-fsdept(1,1,jk)-zdam) / zwam)58 ahm2(jk) = ahm1(jk)77 ahm1(jk) = zmhs + zmhf * TANH( (-fsdept(1,1,jk)-zdam ) / zwam ) 78 ahm2(jk) = zmh2s + zmh2f * TANH( (-fsdept(1,1,jk)-zdam2) / zwam2 ) 59 79 END DO 60 80 … … 62 82 IF(lwp .AND. ld_print ) THEN 63 83 WRITE(numout,*) 64 WRITE(numout,*) ' ahm profile : '84 WRITE(numout,*) ' ahm profile (laplacian): ' 65 85 WRITE(numout,*) 66 86 WRITE(numout,9100) 67 87 DO jk = 1, jpk 68 88 WRITE(numout,9110) jk, ahm1(jk), fsdept(1,1,jk) 89 END DO 90 WRITE(numout,*) 91 WRITE(numout,*) ' ahm profile (bilaplacian): ' 92 WRITE(numout,*) 93 WRITE(numout,9100) 94 DO jk = 1, jpk 95 WRITE(numout,9110) jk, ahm2(jk), fsdept(1,1,jk) 69 96 END DO 70 97 ENDIF
Note: See TracChangeset
for help on using the changeset viewer.