Changeset 3437 for trunk/LMDZ.GENERIC


Ignore:
Timestamp:
Sep 25, 2024, 1:50:44 PM (3 months ago)
Author:
emillour
Message:

Generic PCM Integrate update from newton branch by LT:

  • remove call to old newtrelax routine and replace it with a new one, newtcool_MOCHA, tuned to the MOCHA intercomparison project
Location:
trunk/LMDZ.GENERIC/libf/phystd
Files:
1 added
2 edited

Legend:

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

    r3309 r3437  
    2727  use ioipsl_getin_p_mod, only : getin_p
    2828  use mod_phys_lmdz_para, only : is_parallel, is_master, bcast
     29  use newton_cooling_hotJ, only: planetary_suffix
    2930
    3031!=======================================================================
     
    406407     tau_relax=tau_relax*24*3600 ! convert Earth days --> seconds
    407408
     409     if (is_master) write(*,*) trim(rname)//&
     410     ": planetary suffix for loading file in Newtonian cooling HotJ ?"
     411     planetary_suffix='default_value' ! default value
     412     call getin_p("planetary_suffix",planetary_suffix)
     413     if (is_master) write(*,*) trim(rname)//": planetary_suffix",planetary_suffix
     414     
    408415     if (is_master) write(*,*)trim(rname)//&
    409416       ": call thermal conduction in the soil ?"
  • trunk/LMDZ.GENERIC/libf/phystd/physiq_mod.F90

    r3436 r3437  
    8989      use condensation_generic_mod, only: condensation_generic
    9090      use datafile_mod, only: datadir
     91      use newton_cooling_hotJ, only: newtcool_MOCHA ! LT, adding for MOCHA protocol
     92     
    9193#ifndef MESOSCALE
    9294      use vertical_layers_mod, only: presnivs, pseudoalt
     
    11491151! II.b Call Newtonian cooling scheme
    11501152! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
    1151                call newtrelax(ngrid,nlayer,mu0,sinlat,zpopsk,pt,pplay,pplev,dtrad,firstcall)
     1153               ! call newtrelax(ngrid,nlayer,mu0,sinlat,zpopsk,pt,pplay,pplev,dtrad,firstcall)
     1154               call newtcool_MOCHA(ngrid,nlayer,coslon,coslat,pt,pplay,firstcall,lastcall,dtrad)
    11521155
    11531156               zdtsurf(1:ngrid) = +(pt(1:ngrid,1)-tsurf(1:ngrid))/ptimestep
    11541157               ! e.g. surface becomes proxy for 1st atmospheric layer ?
     1158               fluxrad_sky(1:ngrid) = 0. ! otherwise, not initialized!
    11551159
    11561160            else
     
    11811185         ! Transformation of the radiative tendencies
    11821186         ! ------------------------------------------
     1187
    11831188         zplanck(1:ngrid)=tsurf(1:ngrid)*tsurf(1:ngrid)
    11841189         zplanck(1:ngrid)=emis(1:ngrid)*sigma*zplanck(1:ngrid)*zplanck(1:ngrid)
Note: See TracChangeset for help on using the changeset viewer.