- Timestamp:
- Apr 14, 2015, 9:46:47 AM (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/LMDZ.MARS/libf/dynlonlat_phylonlat/phymars/newstart.F
r1403 r1415 16 16 17 17 use ioipsl_getincom, only: getin 18 use infotrac, only: in iadvtrac, nqtot, tname18 use infotrac, only: infotrac_init, nqtot, tname 19 19 use tracer_mod, only: noms, mmol, 20 20 & igcm_dust_number, igcm_dust_mass, … … 26 26 & co2ice, emis 27 27 use comsoil_h, only: inertiedat, layer, mlayer, nsoilmx, tsoil 28 use control_mod, only: day_step, iphysiq, anneeref 28 use control_mod, only: day_step, iphysiq, anneeref, planet_type 29 29 use phyredem, only: physdem0, physdem1 30 30 use iostart, only: open_startphy … … 179 179 preff = 610. ! for Mars, instead of 101325. (Earth) 180 180 pa= 20 ! for Mars, instead of 500 (Earth) 181 planet_type="mars" 181 182 182 183 ! initialize "serial/parallel" related stuff … … 185 186 186 187 ! Load tracer number and names: 187 call iniadvtrac(nqtot,numvanle) 188 ! call iniadvtrac(nqtot,numvanle) 189 call infotrac_init 188 190 ! allocate arrays 189 191 allocate(q(iip1,jjp1,llm,nqtot)) … … 1628 1630 phis(iip1,:) = phis(1,:) 1629 1631 1630 CALL inidissip ( lstardis, nitergdiv, nitergrot, niterh,1631 * tetagdiv, tetagrot , tetatemp )1632 c CALL inidissip ( lstardis, nitergdiv, nitergrot, niterh, 1633 c * tetagdiv, tetagrot , tetatemp ) 1632 1634 itau=0 1633 1635 if (choix_1.eq.0) then … … 1644 1646 c $ phi,conser,du,dv,dteta,dp,w, pbaru,pbarv, day_ini ) 1645 1647 1646 CALL dynredem0("restart.nc",day_ini, anneeref,phis,nqtot)1648 CALL dynredem0("restart.nc",day_ini,phis) 1647 1649 CALL dynredem1("restart.nc",hour_ini,vcov,ucov,teta,q, 1648 . nqtot,masse,ps)1650 . masse,ps) 1649 1651 C 1650 1652 C Ecriture etat initial physique
Note: See TracChangeset
for help on using the changeset viewer.