Ignore:
Timestamp:
2010-06-24T17:52:01+02:00 (10 years ago)
Author:
acc
Message:

ticket #684 step 4: Add in changes between the head of the dev_r1821_mixed_ldfdyn branch and the trunk@1821

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/DEV_r1784_mid_year_merge_2010/NEMO/OPA_SRC/LDF/ldfdyn_c1d.h90

    r1152 r1954  
    1717      !! ** Method  :   1D eddy viscosity coefficients ( depth ) 
    1818      !!      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 
    2021      !!      isopycnal or geopotential harmonic operator 
    2122      !!                                      : ahm1 defined at T-level 
     
    2728 
    2829      !! * 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 
    3032      REAL(wp) ::   zahmf, zahms 
    3133      !!---------------------------------------------------------------------- 
     
    3638      IF(lwp) WRITE(numout,*) 
    3739 
    38       ! Set ahm1=ahm2   (always at t-level) 
     40      ! Set ahm1 for laplacian     (always at t-level) 
    3941      ! ============= 
    4042      ! (USER: modify ahm1 following your desiderata) 
     
    5355      zmhs = zahms-zmhf * zm00 
    5456 
    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 
    5676      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 ) 
    5979      END DO 
    6080 
     
    6282      IF(lwp .AND. ld_print ) THEN 
    6383         WRITE(numout,*) 
    64          WRITE(numout,*) '         ahm profile : ' 
     84         WRITE(numout,*) '         ahm profile (laplacian): ' 
    6585         WRITE(numout,*) 
    6686         WRITE(numout,9100) 
    6787         DO jk = 1, jpk 
    6888            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) 
    6996         END DO 
    7097      ENDIF 
Note: See TracChangeset for help on using the changeset viewer.