Changeset 4043 for trunk


Ignore:
Timestamp:
Feb 5, 2026, 2:47:26 PM (6 weeks ago)
Author:
jbclement
Message:

Mars PCM:
Addition of 'pa' and 'preff' in the information provided by "start1D.txt".
JBC

Location:
trunk/LMDZ.MARS
Files:
5 edited

Legend:

Unmodified
Added
Removed

  • trunk/LMDZ.MARS/changelog.txt

    r4042 r4043  
    50355035Remove obsolete "user_manual.pdf", reference documentation is the MediaWiki
    50365036https://lmdz-forge.lmd.jussieu.fr/mediawiki/Planets/index.php/Overview_of_the_Mars_PCM
     5037
     5038== 05/02/2026 == JBC
     5039Addition of 'pa' and 'preff' in the information provided by "start1D.txt".

  • trunk/LMDZ.MARS/deftank/start1D.txt

    r3179 r4043  
    1  ps   610.0000000000
    2  co2                              0.0000000000   0.9500000000   0.9500000000   0.9500000000   0.9500000000   0.9500000000   0.9500000000   0.9500000000   0.9500000000   0.9500000000   0.9500000000   0.9500000000   0.9500000000   0.9500000000   0.9500000000   0.9500000000   0.9500000000   0.9500000000   0.9500000000   0.9500000000   0.9500000000   0.9500000000   0.9500000000   0.9500000000   0.9500000000   0.9500000000   0.9500000000   0.9500000000   0.9500000000   0.9500000000
    3  h2o_vap                          0.0000000000   0.0005000000   0.0005000000   0.0005000000   0.0005000000   0.0005000000   0.0005000000   0.0005000000   0.0005000000   0.0005000000   0.0005000000   0.0005000000   0.0005000000   0.0005000000   0.0005000000   0.0005000000   0.0005000000   0.0005000000   0.0005000000   0.0005000000   0.0005000000   0.0005000000   0.0005000000   0.0005000000   0.0005000000   0.0005000000   0.0005000000   0.0005000000   0.0005000000   0.0005000000
    4  h2o_ice                          0.0000000000   0.0000000000   0.0000000000   0.0000000000   0.0000000000   0.0000000000   0.0000000000   0.0000000000   0.0000000000   0.0000000000   0.0000000000   0.0000000000   0.0000000000   0.0000000000   0.0000000000   0.0000000000   0.0000000000   0.0000000000   0.0000000000   0.0000000000   0.0000000000   0.0000000000   0.0000000000   0.0000000000   0.0000000000   0.0000000000   0.0000000000   0.0000000000   0.0000000000   0.0000000000
    5  ccn_number                       0.0021943932   2.1938724518   2.5796332359   3.2571134567   3.4433057308   4.5493345261   4.7790579796   4.3691287041   4.7706074715   4.8893594742   4.4145894051   4.3622183800   5.3330893517   4.3492245674   1.7570885420   0.7927528620   2.8052117825   4.4238777161   5.5149497986   6.0550684929   5.3024525642   3.9112398624   4.6486749649   12.2077741623   22.6576042175   28.4508666992   41.3113174438   11.7739505768   2.6029777527   6.1708893776
    6  ccn_mass                         0.0021943932   2.1938724518   2.5796332359   3.2571134567   3.4433057308   4.5493345261   4.7790579796   4.3691287041   4.7706074715   4.8893594742   4.4145894051   4.3622183800   5.3330893517   4.3492245674   1.7570885420   0.7927528620   2.8052117825   4.4238777161   5.5149497986   6.0550684929   5.3024525642   3.9112398624   4.6486749649   12.2077741623   22.6576042175   28.4508666992   41.3113174438   11.7739505768   2.6029777527   6.1708893776
    7  dust_number                      6746394505248768.0000000000   171839504384.0000000000   161225474048.0000000000   143071084544.0000000000   137828204544.0000000000   144401481728.0000000000   148473561088.0000000000   152183554048.0000000000   149332672512.0000000000   135274577920.0000000000   129245749248.0000000000   144588570624.0000000000   172679544832.0000000000   180114243584.0000000000   156852289536.0000000000   111112560640.0000000000   57837268992.0000000000   19903340544.0000000000   3035456000.0000000000   5885964.0000000000   3328.6264648438   924.7664184570   925.3353881836   1096.6247558594   4722.9399414062   12446.5595703125   28204.7578125000   18738.2050781250   3464.9948730469   887.0166015625
    8  dust_mass                        125.4589462280   0.0026608715   0.0024969140   0.0022266677   0.0021634556   0.0022090040   0.0021991606   0.0021728165   0.0020881307   0.0019683379   0.0018992852   0.0018258542   0.0015930140   0.0011702793   0.0005762299   0.0002119728   0.0000371221   0.0000007424   0.0000000118   0.0000000000   0.0000000000   0.0000000000   0.0000000000   0.0000000000   0.0000000000   0.0000000000   0.0000000000   0.0000000000   0.0000000000   0.0000000000
    9  u   7.5996042233788010   7.6607046721354672   7.8182459568095410   8.1291326500216989   8.6003563912308749   9.1420699610852250   9.5971521808874876   9.8621303873952506   9.9658631538211928   9.9937508498778946   9.9991371629304844   9.9999090223895077   9.9999925678432433   9.9999995189646960   9.9999999746999979   9.9999999988992716   9.9999999999604992   9.9999999999988844   9.9999999999999840   10.000000000000004   9.9999999999999964   9.9999999999999982   10.000000000000002   9.9999999999999982   9.9999999999999964   10.000000000000004   10.000000000000002   10.000000000000004   10.000000000000004
    10  v   0.0000000000000000   0.0000000000000000   0.0000000000000000   0.0000000000000000   0.0000000000000000   0.0000000000000000   0.0000000000000000   0.0000000000000000   0.0000000000000000   0.0000000000000000   0.0000000000000000   0.0000000000000000   0.0000000000000000   0.0000000000000000   0.0000000000000000   0.0000000000000000   0.0000000000000000   0.0000000000000000   0.0000000000000000   0.0000000000000000   0.0000000000000000   0.0000000000000000   0.0000000000000000   0.0000000000000000   0.0000000000000000   0.0000000000000000   0.0000000000000000   0.0000000000000000   0.0000000000000000
    11  teta   255.36959122864869   262.94217129296442   262.86430823135720   262.65847466328324   262.21226857023794   261.37519598684537   259.94715098753250   257.67169898854451   254.34198389155611   250.07191983143383   245.22819265605887   239.99356815622977   234.45448519851908   228.69553782978173   222.71987927059794   216.40633765983554   209.50964495825343   201.82994000256119   193.66575283846379   185.52006334300367   172.11180547112119   154.22292138123350   152.26555642130651   151.86053274326906   153.28065584964054   154.47278979490088   153.89566817514981   151.07020093026608   144.95314172553861   137.21496977935399
     1 ps   610.00000000000000        30.000000000000000        610.00000000000000     
     2 co2                              207.28974396245630       0.95000000000143014       0.95000000000143014       0.95000000000143014       0.95000000000143014       0.95000000000142992       0.95000000000142903       0.95000000000142748       0.95000000000140117     
     3 h2o_ice                          2.0159175267859371        0.0000000000000000        0.0000000000000000        0.0000000000000000        0.0000000000000000        0.0000000000000000        0.0000000000000000        0.0000000000000000        0.0000000000000000     
     4 h2o_vap                          2.0005511063113608        1.0000000000000000E-004   1.0000000000000000E-004   1.0000000000000000E-004   1.0000000000000000E-004   1.0000000000000000E-004   1.0000000000000000E-004   1.0000000000000000E-004   1.0000000000000000E-004
     5 dust_number                      3683127184588657.0        908741027699.70618        908741039954.43994        908740972528.71692        908740463802.42969        908738139867.22363        908730241124.12024        908709894792.02307        908244679635.97949     
     6 dust_mass                        324.19385997885132        9.3577293158506472E-006   9.3577312787075566E-006   9.3577300683319967E-006   9.3577140152716434E-006   9.3576374666951889E-006   9.3573734748091863E-006   9.3566521625726684E-006   9.3415891149355426E-006
     7 u   7.0162859775779278        7.0583328857319874        7.1693033167409466        7.3965987936610089        7.7658238566395656        8.2177899931297080        8.5806634750712991        9.9754019280278978     
     8 v   0.0000000000000000        0.0000000000000000        0.0000000000000000        0.0000000000000000        0.0000000000000000        0.0000000000000000        0.0000000000000000        0.0000000000000000     
     9 teta   147.61966732093049        168.54757895560138        168.62487570105756        168.83484726114642        169.32134171225377        170.24056991939565        171.67221112253304        173.47542102796041        186.22344609723632     

  • trunk/LMDZ.MARS/libf/phymars/dyn1d/init_testphys1d_mod.F90

    r3996 r4043  
    313313! ------------------------
    314314psurf = 610. ! Default value for psurf
     315pa = 20.     ! Default value for transition pressure (for hybrid coord.)
     316preff = 610. ! Default value for reference surface pressure
    315317write(*,*) 'Surface pressure (Pa)?'
    316318if (.not. therestart1D) then
    … …  
    318320else
    319321    open(3,file = trim(start1Dname),status = "old",action = "read")
    320     read(3,*) header, psurf
    321 endif
    322 write(*,*) " psurf = ",psurf
     322    read(3,*) header, psurf, pa, preff
     323endif
     324write(*,*) " psurf = ", psurf
     325write(*,*) " pa    = ", pa
     326write(*,*) " preff = ", preff
    323327
    324328! Compute pressures and altitudes of atmospheric levels
    … …  
    330334call getin("hybrid",hybrid)
    331335write(*,*) " hybrid = ", hybrid
    332 
    333 ! Reference pressures
    334 pa = 20.     ! Transition pressure (for hybrid coord.)
    335 preff = 610. ! Reference surface pressure
    336336
    337337call disvert_noterre

  • trunk/LMDZ.MARS/libf/phymars/dyn1d/testphys1d.F90

    r3912 r4043  
    1111use dimradmars_mod,      only: tauvis, totcloudfrac, albedo
    1212use dust_param_mod,      only: tauscaling
    13 use comvert_mod,         only: ap, bp, aps, bps
     13use comvert_mod,         only: ap, bp, aps, bps, pa, preff
    1414use physiq_mod,          only: physiq
    1515use turb_mod,            only: q2
    … …  
    261261
    262262! Writing the "restart1D.txt" file for the next run
    263 if (startfiles_1D) call writerestart1D('restart1D.txt',psurf,tsurf(1,:),nlayer,size(tsurf,2),temp,u,v,nq,noms,qsurf(1,:,:),q)
     263if (startfiles_1D) call writerestart1D('restart1D.txt',psurf,pa,preff,tsurf(1,:),nlayer,size(tsurf,2),temp,u,v,nq,noms,qsurf(1,:,:),q)
    264264
    265265write(*,*) "testphys1d: everything is cool!"

  • trunk/LMDZ.MARS/libf/phymars/dyn1d/writerestart1D_mod.F90

    r3074 r4043  
    55contains
    66
    7 SUBROUTINE writerestart1D(filename,psurf,tsurf,nlayer,nslope,temp,u,v,nq,qnames,qsurf,q)
     7SUBROUTINE writerestart1D(filename,psurf,pa,preff,tsurf,nlayer,nslope,temp,u,v,nq,qnames,qsurf,q)
    88
    99implicit none
    1010
    1111! Arguments
    12 character(len = *),                intent(in) :: filename
    13 integer,                           intent(in) :: nlayer, nq, nslope
    14 real,                              intent(in) :: psurf
    15 real, dimension(nslope),           intent(in) :: tsurf
    16 real, dimension(nlayer),           intent(in) :: temp, u, v
    17 real, dimension(nlayer,nq),        intent(in) :: q
    18 real, dimension(nq,nslope),        intent(in) :: qsurf
    19 character(len = *), dimension(nq), intent(in) :: qnames
     12character(*),                intent(in) :: filename
     13integer,                     intent(in) :: nlayer, nq, nslope
     14real,                        intent(in) :: psurf, pa, preff
     15real, dimension(nslope),     intent(in) :: tsurf
     16real, dimension(nlayer),     intent(in) :: temp, u, v
     17real, dimension(nlayer,nq),  intent(in) :: q
     18real, dimension(nq,nslope),  intent(in) :: qsurf
     19character(*), dimension(nq), intent(in) :: qnames
    2020
    2121! Local variables
    … …  
    2424! Write the data needed for a restart in "restart1D.txt"
    2525open(1,file = filename,status = "replace",action = "write")
    26 write(1,*) 'ps', psurf
     26write(1,*) 'ps', psurf, pa, preff
    2727do i = 1,nq
    2828    write(1,*) qnames(i), (qsurf(i,j), j = 1,nslope), (q(il,i), il = 1,nlayer)
Note: See TracChangeset for help on using the changeset viewer.