Changeset 3611 for trunk


Ignore:
Timestamp:
Feb 7, 2025, 12:14:55 PM (5 months ago)
Author:
jbclement
Message:

PEM:
Few bug corrections:

  • global surface pressure was badlu updated (r3610);
  • XIOS slope variables were not consistent;
  • missing initialization for slope number (r3609);
  • taking into account that there may be no "startfi.nc" at the beginning of the PEM simulation.

JBC

Location:
trunk
Files:
6 edited

Legend:

Unmodified
Added
Removed

  • trunk/LMDZ.COMMON/libf/evolution/changelog.txt

    r3610 r3611  
    577577- The sublimation flux can now be modified by the new function 'recomp_tend_h2o' to account for the growth of a dust lag layer (see Eran Vos's note for the formula).
    578578- Addition of a function 'itp_tsoil' to compute the soil temperature at any point in the soil profile.
     579
     580== 07/02/2025 == JBC
     581Few bug corrections:
     582    - global surface pressure was badlu updated (r3610);
     583    - XIOS slope variables were not consistent;
     584    - missing initialization for slope number (r3609);
     585    - taking into account that there may be no "startfi.nc" at the beginning of the PEM simulation.

  • trunk/LMDZ.COMMON/libf/evolution/deftank/lib_launchPEM.sh

    r3579 r3611  
    144144    iPEM=1
    145145    iPCM=1
    146     cp startfi.nc starts/
     146    if [ -f "startfi.nc" ]; then
     147        cp startfi.nc starts/
     148    fi
    147149    if [ -f "start.nc" ]; then
    148150        cp start.nc starts/

  • trunk/LMDZ.COMMON/libf/evolution/pem.F90

    r3610 r3611  
    727727    do i = 1,ngrid
    728728        do islope = 1,nslope
    729             ps_avg_global_new = ps_avg_global_old - CO2cond_ps*g*cell_area(i)*d_co2ice(i,islope)*subslope_dist(i,islope)/cos(pi*def_slope_mean(islope)/180.)/total_surface
     729            ps_avg_global_new = ps_avg_global_new - CO2cond_ps*g*cell_area(i)*d_co2ice(i,islope)*subslope_dist(i,islope)/cos(pi*def_slope_mean(islope)/180.)/total_surface
    730730        enddo
    731731    enddo
    … …  
    733733    if (adsorption_pem) then
    734734        do i = 1,ngrid
    735             ps_avg_global_new = ps_avg_global_old - g*cell_area(i)*delta_co2_adsorbed(i)/total_surface
     735            ps_avg_global_new = ps_avg_global_new - g*cell_area(i)*delta_co2_adsorbed(i)/total_surface
    736736        enddo
    737737    endif
    … …  
    12511251call info_PEM(i_myear_leg,stopPEM,i_myear,n_myear)
    12521252
    1253 write(*,*) ">> The PEM has run for", i_myear_leg, "Martian years."
    1254 write(*,*) ">> The chained simulation has run for", i_myear, "Martian years =", i_myear*convert_years, "Earth years."
    1255 write(*,*) ">> The reached date is now", (year_bp_ini + i_myear)*convert_years, "Earth years."
    1256 write(*,*) ">> PEM: so far, so good!"
     1253write(*,*) '************************************************************************'
     1254write(*,*) "> The PEM has run for", i_myear_leg, "Martian years."
     1255write(*,*) "> The chained simulation has run for", i_myear, "Martian years =", i_myear*convert_years, "Earth years."
     1256write(*,*) "> The reached date is now", (year_bp_ini + i_myear)*convert_years, "Earth years."
     1257write(*,*) "> PEM: so far, so good!"
     1258write(*,*) '************************************************************************'
    12571259
    12581260if (layering_algo) then

  • trunk/LMDZ.COMMON/libf/evolution/read_data_PCM_mod.F90

    r3609 r3611  
    7676do islope = 1,nslope
    7777    if (nslope /= 1) then
     78        num='  '
    7879        write(num,'(i2.2)') islope
    7980        num = '_slope'//num
    … …  
    196197do islope = 1,nslope
    197198    if (nslope /= 1) then
     199        num='  '
    198200        write(num,'(i2.2)') islope
    199201        num = '_slope'//num

  • trunk/LMDZ.MARS/deftank/field_def_physics_mars.xml

    r3346 r3611  
    506506                   long_name="waterdensity_surface"
    507507                   unit="kg.m-3" />
    508             <field id="waterdensity_surface01"
     508            <field id="waterdensity_surface_slope01"
    509509                   long_name="waterdensity_surface of slope 01"
    510510                   unit="kg.m-3" />
    511             <field id="waterdensity_surface02"
     511            <field id="waterdensity_surface_slope02"
    512512                   long_name="waterdensity_surface of slope 02"
    513513                   unit="kg.m-3" />
    514             <field id="waterdensity_surface03"
     514            <field id="waterdensity_surface_slope03"
    515515                   long_name="waterdensity_surface of slope 03"
    516516                   unit="kg.m-3" />
    517             <field id="waterdensity_surface04"
     517            <field id="waterdensity_surface_slope04"
    518518                   long_name="waterdensity_surface of slope 04"
    519519                   unit="kg.m-3" />
    520             <field id="waterdensity_surface05"
     520            <field id="waterdensity_surface_slope05"
    521521                   long_name="waterdensity_surface of slope 05"
    522522                   unit="kg.m-3" />
    523             <field id="waterdensity_surface06"
     523            <field id="waterdensity_surface_slope06"
    524524                   long_name="waterdensity_surface of slope 06"
    525525                   unit="kg.m-3" />
    526             <field id="waterdensity_surface07"
     526            <field id="waterdensity_surface_slope07"
    527527                   long_name="waterdensity_surface of slope 07"
    528528                   unit="kg.m-3" />

  • trunk/LMDZ.MARS/deftank/file_def_physics_mars.xml

    r3305 r3611  
    148148                <field field_ref="tsurf_slope06" />
    149149                <field field_ref="tsurf_slope07" />
    150                 <field field_ref="waterdensity_surface01" />
    151                 <field field_ref="waterdensity_surface02" />
    152                 <field field_ref="waterdensity_surface03" />
    153                 <field field_ref="waterdensity_surface04" />
    154                 <field field_ref="waterdensity_surface05" />
    155                 <field field_ref="waterdensity_surface06" />
    156                 <field field_ref="waterdensity_surface07" />
     150                <field field_ref="waterdensity_surface_slope01" />
     151                <field field_ref="waterdensity_surface_slope02" />
     152                <field field_ref="waterdensity_surface_slope03" />
     153                <field field_ref="waterdensity_surface_slope04" />
     154                <field field_ref="waterdensity_surface_slope05" />
     155                <field field_ref="waterdensity_surface_slope06" />
     156                <field field_ref="waterdensity_surface_slope07" />
    157157                <field field_ref="h2o_layer1" />
    158158                <field field_ref="co2_layer1" />
Note: See TracChangeset for help on using the changeset viewer.