Ignore:
Timestamp:
2020-03-23T18:21:59+01:00 (7 months ago)
Author:
gsamson
Message:

revised ABL model version including:

  • albmod cleaning
  • new abl mixing length option (nn_amxl = 3)
  • use rho_air function from aerobulk everywhere
  • remove mxl_abl (replaced by master (mxlm_abl) and dissipative (mxlm_abl) mixing lengths)
  • temporary flag "ln_tpot" to disable potential temperature computation in sbcblk
File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2020/dev_r12563_ASINTER-06_ABL_improvement/src/ABL/abl.F90

    r12489 r12588  
    2929   REAL(wp), PUBLIC, ALLOCATABLE, DIMENSION(:,:,:)     ::   avm_abl      !: turbulent viscosity   [m2/s] 
    3030   REAL(wp), PUBLIC, ALLOCATABLE, DIMENSION(:,:,:)     ::   avt_abl      !: turbulent diffusivity [m2/s] 
    31    REAL(wp), PUBLIC, ALLOCATABLE, DIMENSION(:,:,:)     ::   mxl_abl      !: mixing length         [m] 
     31   REAL(wp), PUBLIC, ALLOCATABLE, DIMENSION(:,:,:)     ::   mxld_abl     !: dissipative mixing length    [m] 
     32   REAL(wp), PUBLIC, ALLOCATABLE, DIMENSION(:,:,:)     ::   mxlm_abl     !: master mixing length         [m] 
    3233   REAL(wp), PUBLIC, ALLOCATABLE, DIMENSION(:,:,:,:)   ::   tke_abl      !: turbulent kinetic energy [m2/s2] 
    3334   REAL(wp), PUBLIC, ALLOCATABLE, DIMENSION(:,:)       ::   fft_abl      !: Coriolis parameter    [1/s] 
     
    3839   REAL(wp), PUBLIC, ALLOCATABLE, DIMENSION(:,:)       ::   msk_abl 
    3940   REAL(wp), PUBLIC, ALLOCATABLE, DIMENSION(:,:)       ::   rest_eq 
     41    
     42   REAL(wp), PUBLIC, ALLOCATABLE, DIMENSION(:,:,:)     ::    cft_abl 
     43   REAL(wp), PUBLIC, ALLOCATABLE, DIMENSION(:,:  )     ::   taux_abl 
     44   REAL(wp), PUBLIC, ALLOCATABLE, DIMENSION(:,:  )     ::   tauy_abl   
    4045   ! 
    4146   INTEGER , PUBLIC :: nt_n, nt_a       !: now / after indices (equal 1 or 2) 
     
    5560      !!---------------------------------------------------------------------- 
    5661      ! 
    57       ALLOCATE( u_abl  (1:jpi,1:jpj,1:jpka,jptime), & 
    58          &      v_abl  (1:jpi,1:jpj,1:jpka,jptime), & 
    59          &      tq_abl (1:jpi,1:jpj,1:jpka,jptime,jptq), & 
    60          &      avm_abl(1:jpi,1:jpj,1:jpka), & 
    61          &      avt_abl(1:jpi,1:jpj,1:jpka), & 
    62          &      mxl_abl(1:jpi,1:jpj,1:jpka), &          
    63          &      tke_abl(1:jpi,1:jpj,1:jpka,jptime), & 
    64          &      fft_abl(1:jpi,1:jpj), & 
    65          &      pblh   (1:jpi,1:jpj), &          
    66          &      msk_abl(1:jpi,1:jpj), & 
    67          &      rest_eq(1:jpi,1:jpj), &          
    68          &      e3t_abl(1:jpka), e3w_abl(1:jpka), ght_abl(1:jpka), ghw_abl(1:jpka),                 STAT=ierr ) 
     62      ALLOCATE( u_abl   (1:jpi,1:jpj,1:jpka,jptime     ), & 
     63         &      v_abl   (1:jpi,1:jpj,1:jpka,jptime     ), & 
     64         &      tq_abl  (1:jpi,1:jpj,1:jpka,jptime,jptq), & 
     65         &      tke_abl (1:jpi,1:jpj,1:jpka,jptime     ), & 
     66         &      avm_abl (1:jpi,1:jpj,1:jpka            ), & 
     67         &      avt_abl (1:jpi,1:jpj,1:jpka            ), & 
     68         &      mxld_abl(1:jpi,1:jpj,1:jpka            ), &          
     69         &      mxlm_abl(1:jpi,1:jpj,1:jpka            ), &  
     70         &      cft_abl (1:jpi,1:jpj,1:jpka            ), &  
     71         &      fft_abl (1:jpi,1:jpj                   ), & 
     72         &      pblh    (1:jpi,1:jpj                   ), &          
     73         &      taux_abl(1:jpi,1:jpj                   ), &    
     74         &      tauy_abl(1:jpi,1:jpj                   ), &         
     75         &      msk_abl (1:jpi,1:jpj                   ), & 
     76         &      rest_eq (1:jpi,1:jpj                   ), &          
     77         &      e3t_abl (1:jpka), e3w_abl(1:jpka)       , & 
     78         &      ght_abl (1:jpka), ghw_abl(1:jpka)       , STAT=ierr ) 
    6979         ! 
    7080      abl_alloc = ierr 
Note: See TracChangeset for help on using the changeset viewer.