Changeset 2714 for trunk/LMDZ.GENERIC


Ignore:
Timestamp:
Jun 21, 2022, 11:05:41 AM (3 years ago)
Author:
aslmd
Message:

Print tendencies in diagfi.nc and small correction

Location:
trunk/LMDZ.GENERIC/libf/phystd
Files:
5 edited

Legend:

Unmodified
Added
Removed
  • trunk/LMDZ.GENERIC/libf/phystd/condensation_generic_mod.F90

    r2711 r2714  
    2121!     Authors
    2222!     -------
    23 !     Adapted from largescale.F90 by Lucas Teinturier (2022)
     23!     Adapted from largescale.F90 by Lucas Teinturier & Noé Clément (2022)
    2424!     largescale.F90 adapted from the LMDTERRE code by R. Wordsworth (2009)
    2525!     Original author Z. X. Li (1993)
  • trunk/LMDZ.GENERIC/libf/phystd/dyn1d/rcm1d.F

    r2713 r2714  
    687687                  q(:,iq)=0. !kg/kg of atmosphere
    688688                  qsurf(iq) = 0. !kg/kg of atmosphere
    689                   ! Look for a "profile_Fe_vap" input file*
     689                  ! Look for a "profile_...._vap" input file
    690690                  tracer_profile_file_name=""
    691691                  tracer_profile_file_name='profile_'//txt
  • trunk/LMDZ.GENERIC/libf/phystd/generic_cloud_common_h.F90

    r2710 r2714  
    11module generic_cloud_common_h
     2
     3    !============================================================================
     4    ! Purpose
     5    ! -------
     6    ! Set up relevant constants and parameters for generic condensable tracers cycles,
     7    ! and generic condensable tracers cloud properties
     8    !
     9    ! Authors
     10    ! -------
     11    ! Adapted from watercommon_h.F90
     12    ! by Lucas Teinturier & Noé Clément (2022)
     13    !============================================================================
     14
    215    use comcstfi_mod, only: r, cpp, mugaz
    316    implicit none
  • trunk/LMDZ.GENERIC/libf/phystd/phys_state_var_mod.F90

    r2550 r2714  
    115115!$OMP THREADPRIVATE(madjdE,madjdEz,lscaledE,lscaledEz)
    116116     
     117      real,dimension(:),allocatable,save :: genericconddE ! Heat from generic condensation (W.m-2)
     118!$OMP THREADPRIVATE(genericconddE)
     119
    117120      real,dimension(:),allocatable,save :: H2Omaxcol ! Maximum possible H2O column amount (at 100% saturation) (kg.m-2).
    118121!$OMP THREADPRIVATE(H2Omaxcol)
     
    202205        ALLOCATE(lscaledE(klon))
    203206        ALLOCATE(lscaledEz(klon,klev))
     207        ALLOCATE(genericconddE(klon))
    204208        ALLOCATE(H2Omaxcol(klon))
    205209        ! This is defined in comsaison_h
     
    293297        DEALLOCATE(lscaledE)
    294298        DEALLOCATE(lscaledEz)
     299        DEALLOCATE(genericconddE)
    295300        DEALLOCATE(H2Omaxcol)       
    296301        DEALLOCATE(mu0)
  • trunk/LMDZ.GENERIC/libf/phystd/physiq_mod.F90

    r2706 r2714  
    15561556            pdq(1:ngrid,1:nlayer,1:nq) = pdq(1:ngrid,1:nlayer,1:nq)+dqvaplscale_generic(1:ngrid,1:nlayer,1:nq)
    15571557            pdq(1:ngrid,1:nlayer,1:nq) = pdq(1:ngrid,1:nlayer,1:nq)+dqcldlscale_generic(1:ngrid,1:nlayer,1:nq)
     1558
     1559            if(enertest)then
     1560               do ig=1,ngrid
     1561                  genericconddE(ig) = cpp*SUM(mass(:,:)*dtlscale_generic(:,:))
     1562               enddo
     1563
     1564               call planetwide_sumval(cpp*massarea(:,:)*dtlscale_generic(:,:)/totarea_planet,dEtot)
     1565
     1566               if (is_master) print*,'In generic condensation atmospheric energy change =',dEtot,' W m-2'
     1567            end if
    15581568         endif !generic_condensation
    15591569
     
    16051615
    16061616         end if ! end of 'sedimentation'
     1617
     1618         write(*,*) 'end of sedimentation pdt', pdt
    16071619
    16081620
     
    22662278!             call writediagfi(ngrid,"h2o_max_col","maximum H2O column amount","kg.m^-2",2,H2Omaxcol)
    22672279
     2280         endif
     2281
     2282         if (generic_condensation) then
     2283
     2284            call writediagfi(ngrid,"genericconddE","heat from generic condensation","W m-2",2,genericconddE)
     2285         
    22682286         endif
    22692287         
Note: See TracChangeset for help on using the changeset viewer.