Changeset 3179


Ignore:
Timestamp:
Jan 18, 2024, 5:23:14 PM (10 months ago)
Author:
jbclement
Message:

Mars PCM:
Improvement of the error message for tracers initialization with a 1D start file + update of "start1D.txt" in the deftank + small cleanings.
JBC

Location:
trunk/LMDZ.MARS
Files:
4 edited

Legend:

Unmodified
Added
Removed

  • trunk/LMDZ.MARS/changelog.txt

    r3176 r3179  
    44394439While at it also from some callphys.def.* files and rename
    44404440callphys.def.watercycle as callphys.def.GCM5 to (hopefully) avoid confusions.
     4441
     4442== 18/01/2024 == JBC
     4443Improvement of the error message for tracers initialization with a 1D start file + update of "start1D.txt" in the deftank + small cleanings.

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

    r3092 r3179  
    1  ps   15844.910967918002     
    2  co2                              0.0000000000000000       0.99806873526322104       0.99806873526322093       0.99806873526322071       0.99806873526322049       0.99806873526321993       0.99806873526321904       0.99806873526321660       0.99806873526321194       0.99806873526320383       0.99806873526318929       0.99806873526316264       0.99806873526311735       0.99806873526304407       0.99806873526292728       0.99806873526274265       0.99806873526245643       0.99806873526200479       0.99806873526124928       0.99806873525982065       0.99806873525618467       0.99806873523022577       0.99806873520941064       0.99806873519925610       0.99806873518649675       0.99806873517068784       0.99806873515200578       0.99806873513094219       0.99806873510863758       0.99806873510863758     
    3  h2o_vap                          7.2725246193759237E-059   5.8820100507640916E-005   5.8820100688353599E-005   5.8820102755633127E-005   5.8820116010573458E-005   5.8820174882493682E-005   5.8820377808845918E-005   5.8820963205398175E-005   5.8822437714700769E-005   5.8825760304503722E-005   5.8832611449481862E-005   5.8845802606717229E-005   5.8869894659825426E-005   5.8912403186450224E-005   5.8986517484320695E-005   5.9118018232116458E-005   5.9339869879177149E-005   5.9118164371551599E-005   4.8106282131383099E-005   2.5397278973802698E-005   6.3061481836226679E-006   9.3802884750568872E-007   9.7516958519586652E-007   1.0089444459472081E-006   1.0459331125368584E-006   1.0833310711125330E-006   1.1209015686700208E-006   1.1577207352404221E-006   1.1925864567746269E-006   1.1925864567746269E-006
    4  h2o_ice                          7.2725246193759237E-059   1.2278899576264690E-057   1.6315181302798064E-057   2.5441790828268903E-057   4.2239110487880329E-055   2.8557485049211623E-053   2.9688260186503111E-049   1.7949152547667348E-043   9.3631379898865452E-040   2.6379228226787604E-038   0.0000000000000000        1.8758830304057299E-025   2.3609586616864940E-016   1.0472188694732230E-014   2.0788822907050202E-013   2.9182716021300108E-011   1.6472227685943954E-010   5.7106584131899193E-008   3.4526182417955551E-006   8.5543632485348008E-006   6.0999521119023684E-006   4.5510207393867132E-007   5.5488154372924425E-009   1.0221427037701964E-010   8.9779605439961759E-013   5.5538516514912759E-015   1.7859437496423380E-016   0.0000000000000000        0.0000000000000000        0.0000000000000000     
    5  ccn_number                       211059535298932.69        1849.4576440198794        1849.4582566053705        1849.4586744268302        1849.4589956111556        1849.4593008861555        1849.4596232822048        1849.4598351982863        1849.4597297222563        1849.4595214127285        1849.4594914571794        1849.3724116475139        1849.5916529580068        1850.0560517435595        1856.5748077355204        279.98827629407265        366.29923119992930        47081.990892276313        12841391.040254256        238487943.23776314        726442991.38934684        1562680885.2700438        24277639.690400742        514973.30873089907        5114.8672520673344        134.33051200098083        333.61003114880236        1755.9927731048235        75.257675414930645        75.257675414930645     
    6  ccn_mass                         984.06104526159436        1.8494576440198791E-017   1.8494582566053700E-017   1.8494586744268300E-017   1.8494589956111556E-017   1.8494593008861554E-017   1.8494596232822046E-017   1.8494598351982859E-017   1.8494597297222565E-017   1.8494595214127284E-017   1.8494594914571792E-017   1.8493724116475137E-017   1.8495916529580064E-017   1.8500560517435601E-017   1.4711524959196136E-015   1.9831509427745010E-013   5.5817990092541446E-013   6.9670583653778338E-011   1.9003293251800345E-008   3.5269180055120037E-007   1.0736789870862268E-006   2.3052678060375135E-006   3.5814177946918155E-008   7.5960801734492301E-010   7.4446918734484585E-012   4.4207908660675703E-014   1.4793308207097557E-015   1.6746694382287890E-017   1.3699270266221628E-018   1.3699270266221628E-018
    7  dust_number                      5978578170980206.0        114369652204.76767        114356549748.81723        114323177293.41280        114250666838.47505        114115899392.29398        113892552929.69894        113552275477.26718        113063989597.62381        112389451883.25011        111509512210.17575        110430121000.77000        109153391663.72237        107673187440.13013        105965302756.97717        103972323657.46086        101572442894.97495        98503858010.630646        94094766154.813080        85779619300.476898        65148281211.421646        6229648260.9592543        814981632.96982491        197504796.99926114        35613192.480257481        4526764.2866460187        407349.75764916284        25148.155974179361        1604.9267530796474        1604.9267530796474     
    8  dust_mass                        501.78869604461329        1.4555946664115001E-003   1.4544322931485786E-003   1.4516334906308965E-003   1.4455755160826605E-003   1.4343809409522598E-003   1.4159936615391615E-003   1.3883839989329550E-003   1.3496583475955031E-003   1.2980270209245006E-003   1.2339122767718103E-003   1.1600958433408263E-003   1.0793391808118474E-003   9.9404666452155318E-004   9.0587859819659137E-004   8.1544923989119642E-004   7.2195266840599168E-004   6.2232326175664016E-004   5.0875482685035374E-004   3.6208336229652212E-004   1.6009516955363409E-004   7.8111971715851268E-006   3.2512390916310862E-007   2.5437139874963564E-008   1.4501610853348091E-009   5.8744110938064371E-011   1.7209978498883730E-012   3.6676996075393140E-014   6.4901243491176011E-016   6.4901243491176011E-016
    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.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

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

    r3174 r3179  
    406406        read(3,*) header, (qsurf(1,iq,j), j = 1,size(qsurf,3)), (q(1,ilayer,iq), ilayer = 1,nlayer)
    407407        if (trim(tname(iq)) /= trim(header)) then
    408             write(*,*) 'Tracer names not compatible for initialization with "'//trim(start1Dname)//'"!'
     408            write(*,*) 'Tracer names between "traceur.def" and "'//trim(start1Dname)//'" do not match!'
     409            write(*,*) 'Please, write the tracer names in the same order for both files.'
    409410            error stop
    410411        endif
    … …  
    498499! Initialize winds for first time step
    499500if (.not. therestart1D) then
    500     u(:) = gru
    501     v(:) = grv
     501    u = gru
     502    v = grv
    502503else
    503504    read(3,*) header, (u(ilayer), ilayer = 1,nlayer)
    … …  
    542543call init_vertical_layers(nlayer,preff,scaleheight,ap,bp,aps,bps,presnivs,pseudoalt)
    543544
    544 plev(:) = ap(:) + psurf*bp(:)
    545 play(:) = aps(:) + psurf*bps(:)
    546 zlay(:) = -200.*r*log(play(:)/plev(1))/g
     545plev = ap + psurf*bp
     546play = aps + psurf*bps
     547zlay = -200.*r*log(play/plev(1))/g
    547548
    548549! Initialize temperature profile
    … …  
    552553! Altitude in km in profile: divide zlay by 1000
    553554tmp1(0) = 0.
    554 tmp1(1:) = zlay(:)/1000.
     555tmp1(1:) = zlay/1000.
    555556
    556557call profile(nlayer + 1,tmp1,tmp2)
    … …  
    609610! Initialize soil content
    610611! -----------------
    611 if (.not. therestartfi) qsoil(:,:,:,:) = 0.
     612if (.not. therestartfi) qsoil = 0.
    612613
    613614! Initialize depths

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

    r3142 r3179  
    187187    ! Compute geopotential
    188188    ! ~~~~~~~~~~~~~~~~~~~~
    189     s(:) = (aps(:)/psurf + bps(:))**rcp
    190     h(:) = cpp*temp(:)/(pks*s(:))
     189    s = (aps/psurf + bps)**rcp
     190    h = cpp*temp/(pks*s)
    191191
    192192    phi(1) = pks*h(1)*(1. - s(1))
    … …  
    203203        ! If atmospheric water is monitored
    204204            if (atm_wat_tau < 0.) then ! Prescribed atm_wat_profile: wet if >0, dry if =0
    205                 q(1,:,igcm_h2o_vap) = min(zqsat(:),atm_wat_profile*g/psurf)
     205                q(1,:,igcm_h2o_vap) = min(zqsat,atm_wat_profile*g/psurf)
    206206                q(1,:,igcm_h2o_ice) = 0. ! reset h2o ice
    207207            else ! Relaxation towards the value atm_wat_profile with relaxation time atm_wat_tau
    208208                q(1,:,igcm_h2o_vap) = atm_wat_profile*g/psurf + (q(1,:,igcm_h2o_vap) - atm_wat_profile*g/psurf)*dexp(-dttestphys/atm_wat_tau)
    209                 q(1,:,igcm_h2o_vap) = min(zqsat(:),q(1,:,igcm_h2o_vap))
     209                q(1,:,igcm_h2o_vap) = min(zqsat,q(1,:,igcm_h2o_vap))
    210210                q(1,:,igcm_h2o_ice) = 0. ! reset h2o ice
    211211            endif
    … …  
    240240    ! For some tests: No coriolis force at equator
    241241    !if (latitude(1) == 0.) then
    242     du(:) = du(:) + (gru - u(:))/1.e4
    243     dv(:) = dv(:) + (grv - v(:))/1.e4
     242    du = du + (gru - u)/1.e4
     243    dv = dv + (grv - v)/1.e4
    244244    !endif
    245245
    … …  
    262262    ! -----------------------------------
    263263    psurf = psurf + dttestphys*dpsurf(1) ! surface pressure change
    264     plev(:) = ap(:) + psurf*bp(:)
    265     play(:) = aps(:) + psurf*bps(:)
     264    plev = ap + psurf*bp
     265    play = aps + psurf*bps
    266266
    267267    ! Increment tracers
Note: See TracChangeset for help on using the changeset viewer.