Ignore:
Timestamp:
2019-12-11T12:38:43+01:00 (20 months ago)
Author:
davestorkey
Message:

2019/dev_r11943_MERGE_2019: Merge in dev_ASINTER-01-05_merge.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2019/dev_r11943_MERGE_2019/src/SAS/diawri.F90

    r11949 r12182  
    2424   !!---------------------------------------------------------------------- 
    2525   USE oce             ! ocean dynamics and tracers  
     26   USE abl            ! abl variables in case ln_abl = .true. 
    2627   USE dom_oce         ! ocean space and time domain 
    2728   USE zdf_oce         ! ocean vertical physics 
     
    5152   PUBLIC   dia_wri_state 
    5253   PUBLIC   dia_wri_alloc           ! Called by nemogcm module 
    53  
     54#if ! defined key_iomput    
     55   PUBLIC   dia_wri_alloc_abl       ! Called by sbcabl  module (if ln_abl = .true.) 
     56#endif 
    5457   INTEGER ::   nid_T, nz_T, nh_T, ndim_T, ndim_hT   ! grid_T file 
    5558   INTEGER ::   nid_U, nz_U, nh_U, ndim_U, ndim_hU   ! grid_U file 
    5659   INTEGER ::   nid_V, nz_V, nh_V, ndim_V, ndim_hV   ! grid_V file 
     60   INTEGER ::   ndim_A, ndim_hA                      ! ABL file    
     61   INTEGER ::   nid_A, nz_A, nh_A                    ! grid_ABL file    
    5762   INTEGER ::   ndex(1)                              ! ??? 
    5863   INTEGER, SAVE, ALLOCATABLE, DIMENSION(:) :: ndex_hT, ndex_hU, ndex_hV 
     64   INTEGER, SAVE, ALLOCATABLE, DIMENSION(:) :: ndex_hA, ndex_A ! ABL 
    5965 
    6066   !! * Substitutions 
     
    115121   END FUNCTION dia_wri_alloc 
    116122    
     123   INTEGER FUNCTION dia_wri_alloc_abl() 
     124      !!---------------------------------------------------------------------- 
     125     ALLOCATE(   ndex_hA(jpi*jpj), ndex_A (jpi*jpj*jpkam1), STAT=dia_wri_alloc_abl) 
     126      CALL mpp_sum( 'diawri', dia_wri_alloc_abl ) 
     127      ! 
     128   END FUNCTION dia_wri_alloc_abl 
    117129   
    118130   SUBROUTINE dia_wri( kt ) 
     
    137149      INTEGER  ::   ierr                                     ! error code return from allocation 
    138150      INTEGER  ::   iimi, iima, ipk, it, itmod, ijmi, ijma   ! local integers 
     151      INTEGER  ::   ipka                                     ! ABL 
    139152      REAL(wp) ::   zsto, zout, zmax, zjulian                ! local scalars 
     153      REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: zw3d_abl   ! ABL 3D workspace 
    140154      !!---------------------------------------------------------------------- 
    141155      ! 
     
    175189      ijmi = 1      ;      ijma = jpj 
    176190      ipk = jpk 
     191     IF(ln_abl) ipka = jpkam1 
    177192 
    178193      ! define time axis 
     
    242257 
    243258         ! No W grid FILE 
     259         IF( ln_abl ) THEN  
     260         ! Define the ABL grid FILE ( nid_A ) 
     261            CALL dia_nam( clhstnam, nwrite, 'grid_ABL' ) 
     262            IF(lwp) WRITE(numout,*) " Name of NETCDF file ", clhstnam    ! filename 
     263            CALL histbeg( clhstnam, jpi, glamt, jpj, gphit,           &  ! Horizontal grid: glamt and gphit 
     264               &          iimi, iima-iimi+1, ijmi, ijma-ijmi+1,       & 
     265               &          nit000-1, zjulian, rdt, nh_A, nid_A, domain_id=nidom, snc4chunks=snc4set ) 
     266            CALL histvert( nid_A, "ght_abl", "Vertical T levels",      &  ! Vertical grid: gdept 
     267               &           "m", ipka, ght_abl(2:jpka), nz_A, "up" ) 
     268            !                                                            ! Index of ocean points 
     269         ALLOCATE( zw3d_abl(jpi,jpj,ipka) )  
     270         zw3d_abl(:,:,:) = 1._wp  
     271         CALL wheneq( jpi*jpj*ipka, zw3d_abl, 1, 1., ndex_A , ndim_A  )      ! volume 
     272            CALL wheneq( jpi*jpj     , zw3d_abl, 1, 1., ndex_hA, ndim_hA )      ! surface 
     273         DEALLOCATE(zw3d_abl) 
     274         ENDIF 
    244275 
    245276         ! Declare all the output fields as NETCDF variables 
     
    262293         CALL histdef( nid_T, "sowindsp", "wind speed at 10m"                  , "m/s"    ,   &  ! wndm 
    263294            &          jpi, jpj, nh_T, 1  , 1, 1  , -99 , 32, clop, zsto, zout ) 
     295! 
     296         IF( ln_abl ) THEN 
     297         !                                                                                      !!! nid_A : 3D 
     298         CALL histdef( nid_A, "t_abl", "Potential Temperature"     , "K"        ,       &  ! t_abl 
     299               &          jpi, jpj, nh_A, ipka, 1, ipka, nz_A, 32, clop, zsto, zout ) 
     300            CALL histdef( nid_A, "q_abl", "Humidity"                  , "kg/kg"    ,       &  ! q_abl 
     301               &          jpi, jpj, nh_A, ipka, 1, ipka, nz_A, 32, clop, zsto, zout )  
     302            CALL histdef( nid_A, "u_abl", "Atmospheric U-wind   "     , "m/s"        ,     &  ! u_abl 
     303               &          jpi, jpj, nh_A, ipka, 1, ipka, nz_A, 32, clop, zsto, zout ) 
     304            CALL histdef( nid_A, "v_abl", "Atmospheric V-wind   "     , "m/s"    ,         &  ! v_abl 
     305               &          jpi, jpj, nh_A, ipka, 1, ipka, nz_A, 32, clop, zsto, zout )  
     306            CALL histdef( nid_A, "tke_abl", "Atmospheric TKE   "     , "m2/s2"    ,        &  ! tke_abl 
     307               &          jpi, jpj, nh_A, ipka, 1, ipka, nz_A, 32, clop, zsto, zout )  
     308            CALL histdef( nid_A, "avm_abl", "Atmospheric turbulent viscosity", "m2/s"   ,  &  ! avm_abl 
     309               &          jpi, jpj, nh_A, ipka, 1, ipka, nz_A, 32, clop, zsto, zout )  
     310            CALL histdef( nid_A, "avt_abl", "Atmospheric turbulent diffusivity", "m2/s2",  &  ! avt_abl 
     311               &          jpi, jpj, nh_A, ipka, 1, ipka, nz_A, 32, clop, zsto, zout )  
     312            CALL histdef( nid_A, "pblh", "Atmospheric boundary layer height "  , "m",      &  ! pblh 
     313               &          jpi, jpj, nh_A,  1  , 1, 1   , -99 , 32, clop, zsto, zout )                  
     314#if defined key_si3 
     315            CALL histdef( nid_A, "oce_frac", "Fraction of open ocean"  , " ",      &  ! ato_i 
     316               &          jpi, jpj, nh_A,  1  , 1, 1   , -99 , 32, clop, zsto, zout ) 
     317#endif 
     318          CALL histend( nid_A, snc4chunks=snc4set ) 
     319       ! 
     320       ENDIF 
     321! 
    264322 
    265323         CALL histend( nid_T, snc4chunks=snc4set ) 
     
    311369      CALL histwrite( nid_T, "soicecov", it, fr_i          , ndim_hT, ndex_hT )   ! ice fraction    
    312370      CALL histwrite( nid_T, "sowindsp", it, wndm          , ndim_hT, ndex_hT )   ! wind speed    
     371! 
     372      IF( ln_abl ) THEN  
     373        ALLOCATE( zw3d_abl(jpi,jpj,jpka) ) 
     374        IF( ln_mskland )   THEN  
     375          DO jk=1,jpka 
     376             zw3d_abl(:,:,jk) = tmask(:,:,1) 
     377            END DO 
     378       ELSE 
     379            zw3d_abl(:,:,:) = 1._wp      
     380         ENDIF        
     381       CALL histwrite( nid_A,  "pblh"   , it, pblh(:,:)                  *zw3d_abl(:,:,1     ), ndim_hA, ndex_hA )   ! pblh  
     382        CALL histwrite( nid_A,  "u_abl"  , it, u_abl   (:,:,2:jpka,nt_n  )*zw3d_abl(:,:,2:jpka), ndim_A , ndex_A  )   ! u_abl 
     383        CALL histwrite( nid_A,  "v_abl"  , it, v_abl   (:,:,2:jpka,nt_n  )*zw3d_abl(:,:,2:jpka), ndim_A , ndex_A  )   ! v_abl 
     384        CALL histwrite( nid_A,  "t_abl"  , it, tq_abl  (:,:,2:jpka,nt_n,1)*zw3d_abl(:,:,2:jpka), ndim_A , ndex_A  )   ! t_abl 
     385        CALL histwrite( nid_A,  "q_abl"  , it, tq_abl  (:,:,2:jpka,nt_n,2)*zw3d_abl(:,:,2:jpka), ndim_A , ndex_A  )   ! q_abl      
     386        CALL histwrite( nid_A,  "tke_abl", it, tke_abl (:,:,2:jpka,nt_n  )*zw3d_abl(:,:,2:jpka), ndim_A , ndex_A  )   ! tke_abl 
     387        CALL histwrite( nid_A,  "avm_abl", it, avm_abl (:,:,2:jpka       )*zw3d_abl(:,:,2:jpka), ndim_A , ndex_A  )   ! avm_abl 
     388        CALL histwrite( nid_A,  "avt_abl", it, avt_abl (:,:,2:jpka       )*zw3d_abl(:,:,2:jpka), ndim_A , ndex_A  )   ! avt_abl   
     389#if defined key_si3 
     390         CALL histwrite( nid_A,  "oce_frac"   , it, ato_i(:,:)                                  , ndim_hA, ndex_hA )   ! ato_i 
     391#endif 
     392       DEALLOCATE(zw3d_abl) 
     393     ENDIF 
     394! 
    313395 
    314396         ! Write fields on U grid 
     
    326408         CALL histclo( nid_U ) 
    327409         CALL histclo( nid_V ) 
     410         IF(ln_abl) CALL histclo( nid_A ) 
    328411      ENDIF 
    329412      ! 
Note: See TracChangeset for help on using the changeset viewer.