Ignore:
Timestamp:
May 5, 2014, 11:38:51 AM (11 years ago)
Author:
aslmd
Message:

MESOSCALE. A necessary complement commit to r1234 where a upgraded interface making use of modules was proposed. Completed the new formulation for module_lmd_driver for newphys with improved interface with both ini/bdy conditions and physical parameterizations. Changed the Registry accordingly. Finished changes about I/O with the LMD physics (see LMDZ.MARS/README). Made all those changes compatible for old interface, and LES runs (checked with test cases), as well as old input files. Changed makemeso to account for full flexibility on changin nx ny ntracers nproc with newphys. Cleaned the now obsolete bits of code used in LMD physics shared with the GCM. ----- Everything is now ready to properly code both restart runs and nesting for mesoscale runs with new physics.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/LMDZ.MARS/libf/phymars/vdif_cd.F

    r1226 r1236  
    11      SUBROUTINE vdif_cd(ngrid,nlay,pz0,
    2      &           pg,pz,pu,pv,wstar,pts,ph,pcdv,pcdh
    3 #ifdef MESOSCALE
    4      &                ,flag_LES
    5 #endif
    6      &                  )
     2     &           pg,pz,pu,pv,wstar,pts,ph,pcdv,pcdh)
    73      USE comcstfi_h
     4      use turb_mod, only: turb_resolved
    85      IMPLICIT NONE
    96c=======================================================================
     
    5249      REAL, INTENT(IN) :: wstar(ngrid)
    5350      REAL, INTENT(OUT) :: pcdv(ngrid),pcdh(ngrid) ! momentum and heat drag coefficient
    54 
    55 #ifdef MESOSCALE
    56       LOGICAL, INTENT(IN) :: flag_LES     !! pour LES avec isfflx!=0
    57 #endif
    5851
    5952c   Local:
     
    160153          zu2(ig)=pu(ig,1)*pu(ig,1) + pv(ig,1)*pv(ig,1)
    161154     &     + (log(1.+0.7*wstar(ig) + 2.3*wstar(ig)**2))**2
    162 #ifdef MESOSCALE
    163           if(flag_LES) then
     155          if(turb_resolved) then
    164156             zu2(ig)=MAX(zu2(ig),1.)
    165157          endif
    166 #endif
    167158!       zu2(ig)=pu(ig,1)*pu(ig,1) + pv(ig,1)*pv(ig,1) + (0.5*wstar(ig))**2
    168159
Note: See TracChangeset for help on using the changeset viewer.