Changeset 8086


Ignore:
Timestamp:
2017-05-30T09:24:22+02:00 (3 years ago)
Author:
cbricaud
Message:

add missing zdf_mxl_alloc in CRS branch

Location:
branches/2015/dev_r5003_MERCATOR6_CRS/NEMOGCM/NEMO
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • branches/2015/dev_r5003_MERCATOR6_CRS/NEMOGCM/NEMO/OPA_SRC/LDF/ldftra_substitute.h90

    r7311 r8086  
    6262#      define   fsahtu_crs(i,j,k)   rldf * aht0 * r_fact_lap_crs(i,j,k) 
    6363#      define   fsahtv_crs(i,j,k)   rldf * aht0 
    64 #      define   fsahtw(_crsi,j,k)   rldf * aht0 
     64#      define   fsahtw_crs(i,j,k)   rldf * aht0 
     65 
    6566#endif 
    6667#endif 
  • branches/2015/dev_r5003_MERCATOR6_CRS/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfmxl.F90

    r7332 r8086  
    9797         IF(lwp) WRITE(numout,*) '~~~~~~~ ' 
    9898         !                             ! allocate zdfmxl arrays 
     99         IF( zdf_mxl_alloc() /= 0 )   CALL ctl_stop( 'STOP', 'zdf_mxl : unable to allocate arrays' ) 
    99100      ENDIF 
    100101 
  • branches/2015/dev_r5003_MERCATOR6_CRS/NEMOGCM/NEMO/TOP_SRC/TRP/trcldf_crs.F90

    r6101 r8086  
    1818   USE trc             ! ocean passive tracers variables 
    1919   USE trcnam_trp      ! passive tracers transport namelist variables 
    20    USE ldftra_oce      ! lateral diffusion coefficient on tracers 
     20   !USE ldftra_oce      ! lateral diffusion coefficient on tracers 
     21   USE ldftra_oce,ONLY: ln_traldf_grif,rn_aht_0,rn_ahtb_0,lk_traldf_eiv     ! lateral diffusion coefficient on tracers 
    2122   USE ldfslp          ! ??? 
    2223   USE traldf_bilapg   ! lateral mixing            (tra_ldf_bilapg routine) 
     
    5859      INTEGER, INTENT( in ) ::   kt   ! ocean time-step index 
    5960      !! 
    60       INTEGER            :: jn 
     61      INTEGER            :: ji,jj,jk,jn 
     62      REAL(wp)           :: zdep 
    6163      CHARACTER (len=22) :: charout 
    6264      REAL(wp), POINTER, DIMENSION(:,:,:,:) ::   ztrtrd 
     
    6870 
    6971      rldf = rldf_rat 
     72      ! 
     73      r_fact_lap(:,:,:) = 1. 
     74      DO jk= 1, jpk 
     75         DO jj = 1, jpj 
     76            DO ji = 1, jpi 
     77               IF( fsdept(ji,jj,jk) > 200. .AND. gphit(ji,jj) < 5. .AND. gphit(ji,jj) > -5. ) THEN 
     78                  zdep = MAX( fsdept(ji,jj,jk) - 1000., 0. ) / 1000. 
     79                  r_fact_lap(ji,jj,jk) = MAX( 1., rn_fact_lap * EXP( -zdep ) ) 
     80               ENDIF 
     81            END DO 
     82         END DO 
     83      END DO 
    7084 
    7185      IF( l_trdtrc )  THEN 
Note: See TracChangeset for help on using the changeset viewer.