Ignore:
Timestamp:
Nov 19, 2013, 12:52:22 PM (11 years ago)
Author:
emillour
Message:

Common dynamics: Updates and modifications to enable running Mars physics with

LMDZ.COMMON dynamics:

  • For compilation: adapted makelmdz, create_make_gcm and makelmdz_fcm, bld.cfg to compile aeronomy routines in "aerono$physique" if it exists, and added "-P -traditional" preprocessing flags in "arch-linux-ifort*"
  • Added function "cbrt.F" (cubic root) in 'bibio'
  • Adapted the reading/writing of dynamics (re)start.nc files for Mars. The main issue is that different information (on time, reference and current) is stored and used differently, hence a few if (planet_type =="mars") here and there. Moreover in the martian case there is the possibility to store fields over multiple times. Some Mars-specific variables (ecritphy,ecritstart,timestart) added in control_mod.F and (hour_ini) in temps.h

EM

File:
1 edited

Legend:

Unmodified
Added
Removed

  • trunk/LMDZ.COMMON/libf/dyn3d/calfis.F

    r1076 r1107  
    728728     .             tracerdyn)      !! tracerdyn <-- utilite ???
    729729
    730       else  ! a moduler pour Mars !!
    731 
    732          CALL physiq (ngridmx,
    733      .             llm,
    734      .             nqtot,
    735      .             debut_split,
    736      .             lafin_split,
    737      .             jD_cur,
    738      .             jH_cur_split,
    739      .             zdt_split,
    740      .             zplev,
    741      .             zplay,
    742      .             zpk,
    743      .             zphi,
    744      .             zphis,
    745      .             presnivs,
    746      .             zufi,
    747      .             zvfi,
    748      .             ztfi,
    749      .             zqfi,
    750      .             flxwfi,
    751      .             zdufi,
    752      .             zdvfi,
    753      .             zdtfi,
    754      .             zdqfi,
    755      .             zdpsrf)
     730      else if ( planet_type=="mars" ) then
     731
     732        CALL physiq (ngridmx,    ! ngrid
     733     .             llm,          ! nlayer
     734     .             nqtot,        ! nq
     735     .             debut_split,  ! firstcall
     736     .             lafin_split,  ! lastcall
     737     .             jD_cur,       ! pday
     738     .             jH_cur_split, ! ptime
     739     .             zdt_split,    ! ptimestep
     740     .             zplev,        ! pplev
     741     .             zplay,        ! pplay
     742     .             zphi,         ! pphi
     743     .             zufi,         ! pu
     744     .             zvfi,         ! pv
     745     .             ztfi,         ! pt
     746     .             zqfi,         ! pq
     747     .             flxwfi,       ! pw
     748     .             zdufi,        ! pdu
     749     .             zdvfi,        ! pdv
     750     .             zdtfi,        ! pdt
     751     .             zdqfi,        ! pdq
     752     .             zdpsrf,       ! pdpsrf
     753     .             tracerdyn)    ! tracerdyn (somewhat obsolete)
     754
     755      else ! unknown "planet_type"
     756
     757        write(lunout,*) "calfis_p: error, unknown planet_type: ",
     758     &                  trim(planet_type)
     759        stop
    756760
    757761      endif ! planet_type
Note: See TracChangeset for help on using the changeset viewer.