Ignore:
Timestamp:
May 11, 2023, 5:40:03 PM (21 months ago)
Author:
romain.vande
Message:

Mars PEM :

Adapt PEM to the subslope PCM configuration, it is now fully compatible.

Create a PEM folder in deftank that contains:

run_pem1: a bash file that runs chained simulation of PEM as well as running a parameterizable number of PCM simulation in between.

It also takes care of reshaping XIOS output as well as renaming outputs etc… in the spirit of run_month1.

It is written for Irene machine and the header needs to be adapted for other machines.

run_PEM.def: A text file that shows the possible parameters to choose before a PEM simulation.

It should be included at the end of run.def just like callphys.def

ob_ex_lsp.asc: An ascii file containing the obliquity, eccentricity, ls_peri data from Laskar in Martian year.
README: A txt file explaining the content of the folder

Adapt field_def_physics_mars.xml to consider the case with 7 subslopes in the PCM.
Change context_lmdz_physics.xml to be able to output the file needed by the PEM.

Correct a few other minor bugs.

RV

Location:
trunk/LMDZ.MARS/deftank
Files:
5 added
2 edited

Legend:

Unmodified
Added
Removed

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

    r2934 r2963  
    129129
    130130        <!-- diurnal averages outputs; output_freq is every day -->
    131         <file id="diurnalave"
    132               name="Xdiurnalave"
     131        <file id="diurnalave_"
     132              name="Xdiurnalave_"
    133133              output_freq="1d"
    134134              time_units="days"
    … …  
    157157            </field_group>
    158158        </file>
     159
     160        <file id="diurnalave"
     161              name="Xdiurnalave"
     162              output_freq="1d"
     163              type="one_file"
     164              time_units="days"
     165              enabled=".false.">
     166                   
     167            <!-- VARS 0D -->
     168            <field_group operation="average"
     169                         freq_op="1ts">
     170                <field field_ref="Ls" />
     171            </field_group>
     172
     173            <field_group operation="average"
     174                         freq_op="1ts">
     175                <field field_ref="area" operation="once" />
     176                <field field_ref="ps" />
     177                <field field_ref="tsurf" />
     178                <field field_ref="tsurf_slope01" />
     179                <field field_ref="tsurf_slope02" />
     180                <field field_ref="tsurf_slope03" />
     181                <field field_ref="tsurf_slope04" />
     182                <field field_ref="tsurf_slope05" />
     183                <field field_ref="tsurf_slope06" />
     184                <field field_ref="tsurf_slope07" />
     185                <field field_ref="Waterdensity_surface01" />
     186                <field field_ref="Waterdensity_surface02" />
     187                <field field_ref="Waterdensity_surface03" />
     188                <field field_ref="Waterdensity_surface04" />
     189                <field field_ref="Waterdensity_surface05" />
     190                <field field_ref="Waterdensity_surface06" />
     191                <field field_ref="Waterdensity_surface07" />
     192                <field field_ref="h2o_layer1" />
     193                <field field_ref="co2_layer1" />
     194            </field_group>
     195
     196            <field_group operation="minimum"
     197                         freq_op="1ts">
     198                <field field_ref="co2ice" />
     199                <field field_ref="co2ice_slope01" />
     200                <field field_ref="co2ice_slope02" />
     201                <field field_ref="co2ice_slope03" />
     202                <field field_ref="co2ice_slope04" />
     203                <field field_ref="co2ice_slope05" />
     204                <field field_ref="co2ice_slope06" />
     205                <field field_ref="co2ice_slope07" />
     206                <field field_ref="h2o_ice_s" />
     207                <field field_ref="h2o_ice_s_slope01" />
     208                <field field_ref="h2o_ice_s_slope02" />
     209                <field field_ref="h2o_ice_s_slope03" />
     210                <field field_ref="h2o_ice_s_slope04" />
     211                <field field_ref="h2o_ice_s_slope05" />
     212                <field field_ref="h2o_ice_s_slope06" />
     213                <field field_ref="h2o_ice_s_slope07" />
     214                <field field_ref="watercap_slope01" />
     215                <field field_ref="watercap_slope02" />
     216                <field field_ref="watercap_slope03" />
     217                <field field_ref="watercap_slope04" />
     218                <field field_ref="watercap_slope05" />
     219                <field field_ref="watercap_slope06" />
     220                <field field_ref="watercap_slope07" />
     221            </field_group>
     222
     223            <!-- VARS soil -->
     224            <field_group operation="average"
     225                         freq_op="1ts">
     226                <field field_ref="soiltemp" />
     227                <field field_ref="soiltemp_slope01" />
     228                <field field_ref="soiltemp_slope02" />
     229                <field field_ref="soiltemp_slope03" />
     230                <field field_ref="soiltemp_slope04" />
     231                <field field_ref="soiltemp_slope05" />
     232                <field field_ref="soiltemp_slope06" />
     233                <field field_ref="soiltemp_slope07" />
     234                <field field_ref="Waterdensity_soil_slope01" />
     235                <field field_ref="Waterdensity_soil_slope02" />
     236                <field field_ref="Waterdensity_soil_slope03" />
     237                <field field_ref="Waterdensity_soil_slope04" />
     238                <field field_ref="Waterdensity_soil_slope05" />
     239                <field field_ref="Waterdensity_soil_slope06" />
     240                <field field_ref="Waterdensity_soil_slope07" />
     241            </field_group>
     242        </file>
     243
     244
    159245    </file_definition>
    160246

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

    r2934 r2963  
    6868                   long_name="Surface emissivity of slope 01" 
    6969                   unit="" />
     70            <field id="emis_slope02"
     71                   long_name="Surface emissivity of slope 02" 
     72                   unit="" />
     73            <field id="emis_slope03"
     74                   long_name="Surface emissivity of slope 03" 
     75                   unit="" />
     76            <field id="emis_slope04"
     77                   long_name="Surface emissivity of slope 04" 
     78                   unit="" />
     79            <field id="emis_slope05"
     80                   long_name="Surface emissivity of slope 05" 
     81                   unit="" />
     82            <field id="emis_slope06"
     83                   long_name="Surface emissivity of slope 06" 
     84                   unit="" />
     85            <field id="emis_slope07"
     86                   long_name="Surface emissivity of slope 07" 
     87                   unit="" />
    7088            <field id="albedo"
    7189                   long_name="Albedo of the surface" 
    … …  
    7492                   long_name="Albedo of the surface for slope 01" 
    7593                   unit="" />
    76 
     94            <field id="albedo_slope02"
     95                   long_name="Albedo of the surface for slope 02" 
     96                   unit="" />
     97            <field id="albedo_slope03"
     98                   long_name="Albedo of the surface for slope 03" 
     99                   unit="" />
     100            <field id="albedo_slope04"
     101                   long_name="Albedo of the surface for slope 04" 
     102                   unit="" />
     103            <field id="albedo_slope05"
     104                   long_name="Albedo of the surface for slope 05" 
     105                   unit="" />
     106            <field id="albedo_slope06"
     107                   long_name="Albedo of the surface for slope 06" 
     108                   unit="" />
     109            <field id="albedo_slope07"
     110                   long_name="Albedo of the surface for slope 07" 
     111                   unit="" />
    77112            <field id="local_time"
    78113                   long_name="Local time"
    … …  
    88123                   long_name="Surface Temperature of slope 01"
    89124                   unit="K" />
     125            <field id="tsurf_slope02"
     126                   long_name="Surface Temperature of slope 02"
     127                   unit="K" />
     128            <field id="tsurf_slope03"
     129                   long_name="Surface Temperature of slope 03"
     130                   unit="K" />
     131            <field id="tsurf_slope04"
     132                   long_name="Surface Temperature of slope 04"
     133                   unit="K" />
     134            <field id="tsurf_slope05"
     135                   long_name="Surface Temperature of slope 05"
     136                   unit="K" />
     137            <field id="tsurf_slope06"
     138                   long_name="Surface Temperature of slope 06"
     139                   unit="K" />
     140            <field id="tsurf_slope07"
     141                   long_name="Surface Temperature of slope 07"
     142                   unit="K" />
    90143            <field id="temp_layer1"
    91144                   long_name="Temperature in layer 1"
    … …  
    102155                   long_name="Longwave radiation at the surface on slope 01"
    103156                   unit="W.m-2" />
     157            <field id="fluxsurf_lw_slope02"
     158                   long_name="Longwave radiation at the surface on slope 02"
     159                   unit="W.m-2" />
     160            <field id="fluxsurf_lw_slope03"
     161                   long_name="Longwave radiation at the surface on slope 03"
     162                   unit="W.m-2" />
     163            <field id="fluxsurf_lw_slope04"
     164                   long_name="Longwave radiation at the surface on slope 04"
     165                   unit="W.m-2" />
     166            <field id="fluxsurf_lw_slope05"
     167                   long_name="Longwave radiation at the surface on slope 05"
     168                   unit="W.m-2" />
     169            <field id="fluxsurf_lw_slope06"
     170                   long_name="Longwave radiation at the surface on slope 06"
     171                   unit="W.m-2" />
     172            <field id="fluxsurf_lw_slope07"
     173                   long_name="Longwave radiation at the surface on slope 07"
     174                   unit="W.m-2" />
    104175            <field id="fluxtop_lw"
    105176                   long_name="Longwave radiation at the top of the atmosphere"
    … …  
    111182            <field id="fluxsurf_dn_sw_slope01"
    112183                   long_name="Incoming shortwave radiation at the surface on slope 01" 
     184                   unit="W.m-2" />
     185            <field id="fluxsurf_dn_sw_slope02"
     186                   long_name="Incoming shortwave radiation at the surface on slope 02" 
     187                   unit="W.m-2" />
     188            <field id="fluxsurf_dn_sw_slope03"
     189                   long_name="Incoming shortwave radiation at the surface on slope 03" 
     190                   unit="W.m-2" />
     191            <field id="fluxsurf_dn_sw_slope04"
     192                   long_name="Incoming shortwave radiation at the surface on slope 04" 
     193                   unit="W.m-2" />
     194            <field id="fluxsurf_dn_sw_slope05"
     195                   long_name="Incoming shortwave radiation at the surface on slope 05" 
     196                   unit="W.m-2" />
     197            <field id="fluxsurf_dn_sw_slope06"
     198                   long_name="Incoming shortwave radiation at the surface on slope 06" 
     199                   unit="W.m-2" />
     200            <field id="fluxsurf_dn_sw_slope07"
     201                   long_name="Incoming shortwave radiation at the surface on slope 07" 
    113202                   unit="W.m-2" />
    114203            <field id="fluxtop_up_sw"
    … …  
    145234                   long_name="tendency due to vertical diffusion of background dust on surface"
    146235                   unit="kg.m-2.s-1" />
     236            <field id="dqsdifdustq_slope02"
     237                   long_name="tendency due to vertical diffusion of background dust on surface"
     238                   unit="kg.m-2.s-1" />
     239            <field id="dqsdifdustq_slope03"
     240                   long_name="tendency due to vertical diffusion of background dust on surface"
     241                   unit="kg.m-2.s-1" />
     242            <field id="dqsdifdustq_slope04"
     243                   long_name="tendency due to vertical diffusion of background dust on surface"
     244                   unit="kg.m-2.s-1" />
     245            <field id="dqsdifdustq_slope05"
     246                   long_name="tendency due to vertical diffusion of background dust on surface"
     247                   unit="kg.m-2.s-1" />
     248            <field id="dqsdifdustq_slope06"
     249                   long_name="tendency due to vertical diffusion of background dust on surface"
     250                   unit="kg.m-2.s-1" />
     251            <field id="dqsdifdustq_slope07"
     252                   long_name="tendency due to vertical diffusion of background dust on surface"
     253                   unit="kg.m-2.s-1" />
    147254            <field id="dqsdifrdsq"
    148255                   long_name="tendency due to vertical diffusion of stormdust on surface"
    149256                   unit="kg.m-2.s-1" />
    150257            <field id="dqsdifrdsq_slope01"
     258                   long_name="tendency due to vertical diffusion of stormdust on surface"
     259                   unit="kg.m-2.s-1" />
     260            <field id="dqsdifrdsq_slope02"
     261                   long_name="tendency due to vertical diffusion of stormdust on surface"
     262                   unit="kg.m-2.s-1" />
     263            <field id="dqsdifrdsq_slope03"
     264                   long_name="tendency due to vertical diffusion of stormdust on surface"
     265                   unit="kg.m-2.s-1" />
     266            <field id="dqsdifrdsq_slope04"
     267                   long_name="tendency due to vertical diffusion of stormdust on surface"
     268                   unit="kg.m-2.s-1" />
     269            <field id="dqsdifrdsq_slope05"
     270                   long_name="tendency due to vertical diffusion of stormdust on surface"
     271                   unit="kg.m-2.s-1" />
     272            <field id="dqsdifrdsq_slope06"
     273                   long_name="tendency due to vertical diffusion of stormdust on surface"
     274                   unit="kg.m-2.s-1" />
     275            <field id="dqsdifrdsq_slope07"
    151276                   long_name="tendency due to vertical diffusion of stormdust on surface"
    152277                   unit="kg.m-2.s-1" />
    … …  
    185310            <field id="watercap_slope01"
    186311                   long_name="Perennial water ice thickness of slope 01"
    187                    unit="kg.m-2" />           
     312                   unit="kg.m-2" /> 
     313            <field id="watercap_slope02"
     314                   long_name="Perennial water ice thickness of slope 02"
     315                   unit="kg/m2" />
     316            <field id="watercap_slope03"
     317                   long_name="Perennial water ice thickness of slope 03"
     318                   unit="kg/m2" />
     319            <field id="watercap_slope04"
     320                   long_name="Perennial water ice thickness of slope 04"
     321                   unit="kg/m2" />
     322            <field id="watercap_slope05"
     323                   long_name="Perennial water ice thickness of slope 05"
     324                   unit="kg/m2" />
     325            <field id="watercap_slope06"
     326                   long_name="Perennial water ice thickness of slope 06"
     327                   unit="kg/m2" />
     328            <field id="watercap_slope07"
     329                   long_name="Perennial water ice thickness of slope 07"
     330                   unit="kg/m2" />
     331   
    188332            <field id="surf_h2o_lh"
    189333                   long_name="Ground ice latent heat flux"
    … …  
    213357                   long_name="Mass of water ice on the surface of slope 01"
    214358                   unit="kg.m-2" />
     359            <field id="h2o_ice_s_slope02"
     360                   long_name="Mass of water ice on the surface of slope 02"
     361                   unit="kg/m2" />
     362            <field id="h2o_ice_s_slope03"
     363                   long_name="Mass of water ice on the surface of slope 03"
     364                   unit="kg/m2" />
     365            <field id="h2o_ice_s_slope04"
     366                   long_name="Mass of water ice on the surface of slope 04"
     367                   unit="kg/m2" />
     368            <field id="h2o_ice_s_slope05"
     369                   long_name="Mass of water ice on the surface of slope 05"
     370                   unit="kg/m2" />
     371            <field id="h2o_ice_s_slope06"
     372                   long_name="Mass of water ice on the surface of slope 06"
     373                   unit="kg/m2" />
     374            <field id="h2o_ice_s_slope07"
     375                   long_name="Mass of water ice on the surface of slope 07"
     376                   unit="kg/m2" />
     377            <field id="Waterdensity_surface01"
     378                   long_name="Waterdensity_surface of slope 01"
     379                   unit="XX" />
     380            <field id="Waterdensity_surface02"
     381                   long_name="Waterdensity_surface of slope 02"
     382                   unit="XX" />
     383            <field id="Waterdensity_surface03"
     384                   long_name="Waterdensity_surface of slope 03"
     385                   unit="XX" />
     386            <field id="Waterdensity_surface04"
     387                   long_name="Waterdensity_surface of slope 04"
     388                   unit="XX" />
     389            <field id="Waterdensity_surface05"
     390                   long_name="Waterdensity_surface of slope 05"
     391                   unit="XX" />
     392            <field id="Waterdensity_surface06"
     393                   long_name="Waterdensity_surface of slope 06"
     394                   unit="XX" />
     395            <field id="Waterdensity_surface07"
     396                   long_name="Waterdensity_surface of slope 07"
     397                   unit="XX" />
     398            <field id="h2o_layer1"
     399                   long_name="h2o in the first layer" 
     400                   unit="kg/kg" />
     401            <field id="co2_layer1"
     402                   long_name="co2 in the first layer" 
     403                   unit="kg/kg" />
     404
     405
    215406           
    216407            <!-- CO2 cycle -->
    … …  
    221412                   long_name="CO2 ice thickness on slope 01"
    222413                   unit="kg.m-2" />
     414            <field id="co2ice_slope02"
     415                   long_name="CO2 ice thickness on slope 02"
     416                   unit="kg/m2" />
     417            <field id="co2ice_slope03"
     418                   long_name="CO2 ice thickness on slope 03"
     419                   unit="kg/m2" />
     420            <field id="co2ice_slope04"
     421                   long_name="CO2 ice thickness on slope 04"
     422                   unit="kg/m2" />
     423            <field id="co2ice_slope05"
     424                   long_name="CO2 ice thickness on slope 05"
     425                   unit="kg/m2" />
     426            <field id="co2ice_slope06"
     427                   long_name="CO2 ice thickness on slope 06"
     428                   unit="kg/m2" />
     429            <field id="co2ice_slope07"
     430                   long_name="CO2 ice thickness on slope 07"
     431                   unit="kg/m2" />
     432
    223433           
    224434            <field id="co2condens_zfallice"
    … …  
    463673                   long_name="Soil temperature for slope 01"
    464674                   unit="K" />
     675            <field id="soiltemp_slope02"
     676                   long_name="Soil temperature for slope 02"
     677                   unit="K" />
     678            <field id="soiltemp_slope03"
     679                   long_name="Soil temperature for slope 03"
     680                   unit="K" />
     681            <field id="soiltemp_slope04"
     682                   long_name="Soil temperature for slope 04"
     683                   unit="K" />
     684            <field id="soiltemp_slope05"
     685                   long_name="Soil temperature for slope 05"
     686                   unit="K" />
     687            <field id="soiltemp_slope06"
     688                   long_name="Soil temperature for slope 06"
     689                   unit="K" />
     690            <field id="soiltemp_slope07"
     691                   long_name="Soil temperature for slope 07"
     692                   unit="K" />
    465693            <field id="inertiedat"
    466694                   long_name="Soil thermal inertia"
    467695                   unit="J/kg/K" />
     696            <field id="inertiesoil_slope01"
     697                   long_name="Soil thermal inertia for slope 01"
     698                   unit="J/kg/K" />
     699            <field id="inertiesoil_slope02"
     700                   long_name="Soil thermal inertia for slope 02"
     701                   unit="J/kg/K" />
     702            <field id="inertiesoil_slope03"
     703                   long_name="Soil thermal inertia for slope 03"
     704                   unit="J/kg/K" />
     705            <field id="inertiesoil_slope04"
     706                   long_name="Soil thermal inertia for slope 04"
     707                   unit="J/kg/K" />
     708            <field id="inertiesoil_slope05"
     709                   long_name="Soil thermal inertia for slope 05"
     710                   unit="J/kg/K" />
     711            <field id="inertiesoil_slope06"
     712                   long_name="Soil thermal inertia for slope 06"
     713                   unit="J/kg/K" />
     714            <field id="inertiesoil_slope07"
     715                   long_name="Soil thermal inertia for slope 07"
     716                   unit="J/kg/K" />
     717            <field id="Waterdensity_soil_slope01 for slope 01"
     718                   long_name="rhowater_soil"
     719                   unit="J/kg/K" />
     720            <field id="Waterdensity_soil_slope02 for slope 02"
     721                   long_name="rhowater_soil"
     722                   unit="J/kg/K" />
     723            <field id="Waterdensity_soil_slope03 for slope 03"
     724                   long_name="rhowater_soil"
     725                   unit="J/kg/K" />
     726            <field id="Waterdensity_soil_slope04 for slope 04"
     727                   long_name="rhowater_soil"
     728                   unit="J/kg/K" />
     729            <field id="Waterdensity_soil_slope05 for slope 05"
     730                   long_name="rhowater_soil"
     731                   unit="J/kg/K" />
     732            <field id="Waterdensity_soil_slope06 for slope 06"
     733                   long_name="rhowater_soil"
     734                   unit="J/kg/K" />
     735            <field id="Waterdensity_soil_slope07 for slope 07"
     736                   long_name="rhowater_soil"
     737                   unit="J/kg/K" />
     738
    468739        </field_group>
    469740
Note: See TracChangeset for help on using the changeset viewer.