Skip to content
Snippets Groups Projects
turb.f90 82.5 KiB
Newer Older
  • Learn to ignore specific revisions
  •   ZLM(:,:) = ZCOEF_AMPL(:,:)*ZLM_CLOUD(:,:)
    
    !$mnh_end_expand_where(JIJ=IIJB:IIJE,JK=1:IKT)
    
    !
    ! Cloud mixing length in the clouds at the points which do not verified the CEI
    !
    
    !$mnh_expand_where(JIJ=IIJB:IIJE,JK=1:IKT)
    
    WHERE (PCEI(:,:) == -1.)
      ZLM(:,:) = ZLM_CLOUD(:,:)
    
    !$mnh_end_expand_where(JIJ=IIJB:IIJE,JK=1:IKT)
    
    IF ( TURBN%LTURB_DIAG .AND. TPFILE%LOPENED ) THEN
    
      TZFIELD = TFIELDMETADATA(         &
        CMNHNAME   = 'COEF_AMPL',       &
        CSTDNAME   = '',                &
        CLONGNAME  = 'COEF_AMPL',       &
        CUNITS     = '1',               &
        CDIR       = 'XY',              &
        CCOMMENT   = 'X_Y_Z_COEF AMPL', &
        NGRID      = 1,                 &
        NTYPE      = TYPEREAL,          &
        NDIMS      = 3,                 &
        LTIMEDEP   = .TRUE.             )
    
      CALL IO_FIELD_WRITE_PHY(D,TPFILE,TZFIELD,ZCOEF_AMPL)
    
      TZFIELD = TFIELDMETADATA(        &
        CMNHNAME   = 'LM_CLOUD',       &
        CSTDNAME   = '',               &
        CLONGNAME  = 'LM_CLOUD',       &
        CUNITS     = 'm',              &
        CDIR       = 'XY',             &
        CCOMMENT   = 'X_Y_Z_LM CLOUD', &
        NGRID      = 1,                &
        NTYPE      = TYPEREAL,         &
        NDIMS      = 3,                &
        LTIMEDEP   = .TRUE.            )
    
      CALL IO_FIELD_WRITE_PHY(D,TPFILE,TZFIELD,ZLM_CLOUD)
    
    IF (LHOOK) CALL DR_HOOK('TURB:CLOUD_MODIF_LM',1,ZHOOK_HANDLE2)
    
    END SUBROUTINE TURB