Changeset 4187 for dynamico_lmdz


Ignore:
Timestamp:
Dec 19, 2019, 10:26:38 PM (5 years ago)
Author:
dubos
Message:

simple_physics : cleanup

Location:
dynamico_lmdz/simple_physics/phyparam
Files:
3 edited

Legend:

Unmodified
Added
Removed

  • dynamico_lmdz/simple_physics/phyparam/param/vdif.F

    r4186 r4187  
    77     $                lwrite)
    88      USE phys_const
    9       USE planet, ONLY : karman
    109      USE vdif_mod
    1110      IMPLICIT NONE

  • dynamico_lmdz/simple_physics/phyparam/physics/planet.F90

    r4184 r4187  
    77       lmixmin=100., emin_turb=1e-8
    88 
    9   REAL, PARAMETER :: karman=0.4
    10 
    119END MODULE planet

  • dynamico_lmdz/simple_physics/phyparam/physics/vdif_mod.F90

    r4186 r4187  
    22  IMPLICIT NONE
    33  SAVE
     4  PRIVATE
    45
     6  REAL, PARAMETER :: karman=0.4
     7
     8  PUBLIC :: vdif_cd, vdif_k
     9 
    510CONTAINS
    611
    712  SUBROUTINE vdif_cd( ngrid,pz0,pg,pz,pu,pv,pts,ph,pcdv,pcdh)
    8     USE planet, ONLY : karman
    913    !=======================================================================
    1014    !
    … …  
    5256    !   ------
    5357   
     58    REAL, PARAMETER :: b=5., c=5., d=5., umin2=1e-12, &
     59         c2b=2.*b, c3bc=3.*b*c, c3b=3.*b
     60 
    5461    INTEGER ig
    5562   
    5663    REAL zu2,z1,zri,zcd0,zz
    5764   
    58     REAL b,c,d,c2b,c3bc,c3b,z0,umin2
    59     LOGICAL firstcal
    60     DATA b,c,d,umin2/5.,5.,5.,1.e-12/
    61     DATA firstcal/.true./
    62     SAVE b,c,d,c2b,c3bc,c3b,firstcal,umin2
    6365   
    6466    !-----------------------------------------------------------------------
    … …  
    7375    !     RETURN
    7476   
    75     IF (firstcal) THEN
    76        c2b=2.*b
    77        c3bc=3.*b*c
    78        c3b=3.*b
    79        firstcal=.false.
    80     ENDIF
    8177   
    8278!!!! WARNING, verifier la formule originale de Louis!
    … …  
    106102  SUBROUTINE vdif_k(ngrid,nlay,   &
    107103       ptimestep,pg,pzlev,pzlay,pz0,pu,pv,ph,pcdv,pkv,pkh)
    108     ! FIXME : pkh :: pkv
     104    ! FIXME : pkh := pkv
    109105    USE planet
    110106    INTEGER ngrid,nlay
Note: See TracChangeset for help on using the changeset viewer.