Ignore:
Timestamp:
Mar 26, 2025, 11:05:58 AM (5 months ago)
Author:
emillour
Message:

Pluto PCM:
Some fixes to enable runnnig with dynamico:

  • add "strictboundcorrk" flag in callcorrk_pluto to enable running even if outside of kmatrix temperatures (when strictboundcorrk=.true.)
  • add premature exiting of writediagsoil if not with lon-lat grid
  • while at it, turned surfprop.F90 into a module
  • in physiq, enforce the possibility to output subsurface-related field in most cases, not only when "fast=.true."
  • adapt reference xml files: subsurface quantities need to be defined on a dedicated grid, otherwise XIOS will generate misleading garbage values. Updated files are put in "deftank/dynamico" for now.

EM

File:
1 edited

Legend:

Unmodified
Added
Removed

  • trunk/LMDZ.PLUTO/deftank/dynamico/context_pcm_physics.xml

    r3544 r3698  
    1010
    1111    </domain_definition>
    12 
    13     <domain_definition>
    14         <domain id="dom_64_48" type="rectilinear" ni_glo="64" nj_glo="48"   >
    15           <generate_rectilinear_domain/>
    16           <interpolate_domain order="1"/>
    17         </domain>
    18 
    19         <domain id="dom_96_95" type="rectilinear"  ni_glo="96" nj_glo="95"  >
    20           <generate_rectilinear_domain/>
    21           <interpolate_domain order="1"/>
    22         </domain>
    23 
    24         <domain id="dom_144_142" type="rectilinear" ni_glo="144" nj_glo="142"   >
    25           <generate_rectilinear_domain/>
    26           <interpolate_domain order="1"/>
    27         </domain>
    28 
    29         <domain id="dom_512_360" type="rectilinear" ni_glo="512" nj_glo="360"   >
    30           <generate_rectilinear_domain/>
    31           <interpolate_domain order="1"/>
    32         </domain>
    33 
    34         <domain id="dom_720_360" type="rectilinear" ni_glo="720" nj_glo="360" >
    35           <generate_rectilinear_domain/>
    36           <interpolate_domain order="1"/>
    37         </domain>
    38 
    39         <domain id="dom_128_96" type="rectilinear" ni_glo="128" nj_glo="96" >
    40           <generate_rectilinear_domain/>
    41           <interpolate_domain order="1"/>
    42         </domain>
    43 
    44         <domain id="dom_256_192" type="rectilinear" ni_glo="256" nj_glo="192" >
    45           <generate_rectilinear_domain/>
    46           <interpolate_domain order="1"/>
    47         </domain>
    48 
    49         <domain id="dom_out" domain_ref="dom_128_96"/>
    50     </domain_definition>
    51     <!-- <domain_definition src="./domain_def.xml"/> -->
    5212
    5313    <!-- Definition of vertical axes
    … …  
    9959            <axis axis_ref="altitude" />
    10060        </grid>
     61        <grid id="grid_soil">
     62            <domain domain_ref="dom_glo" />
     63            <axis axis_ref="soil_layers" />
     64        </grid>
    10165        <grid id="grid_2D">
    10266            <domain domain_ref="dom_glo" />
    … …  
    12589            <domain domain_ref="dom_out" />
    12690        </grid>
     91        <grid id="grid_soil_out">
     92            <domain domain_ref="dom_out" />
     93            <axis axis_ref="soil_layers" />
     94        </grid>
    12795        <!--added for diagspecIR-->
    12896        <grid id="grid_IR_3D_out">
    … …  
    143111    </grid_definition>
    144112
     113    <!-- various predefined grids to output onto -->
     114    <domain_definition>
     115        <domain id="dom_64_48" type="rectilinear" ni_glo="64" nj_glo="48"   >
     116          <generate_rectilinear_domain/>
     117          <interpolate_domain order="1"/>
     118        </domain>
     119
     120        <domain id="dom_96_95" type="rectilinear"  ni_glo="96" nj_glo="95"  >
     121          <generate_rectilinear_domain/>
     122          <interpolate_domain order="1"/>
     123        </domain>
     124
     125        <domain id="dom_144_142" type="rectilinear" ni_glo="144" nj_glo="142"   >
     126          <generate_rectilinear_domain/>
     127          <interpolate_domain order="1"/>
     128        </domain>
     129
     130        <domain id="dom_512_360" type="rectilinear" ni_glo="512" nj_glo="360"   >
     131          <generate_rectilinear_domain/>
     132          <interpolate_domain order="1"/>
     133        </domain>
     134
     135        <domain id="dom_720_360" type="rectilinear" ni_glo="720" nj_glo="360" >
     136          <generate_rectilinear_domain/>
     137          <interpolate_domain order="1"/>
     138        </domain>
     139
     140        <domain id="dom_128_96" type="rectilinear" ni_glo="128" nj_glo="96" >
     141          <generate_rectilinear_domain/>
     142          <interpolate_domain order="1"/>
     143        </domain>
     144
     145        <domain id="dom_256_192" type="rectilinear" ni_glo="256" nj_glo="192" >
     146          <generate_rectilinear_domain/>
     147          <interpolate_domain order="1"/>
     148        </domain>
     149
     150        <!-- set dom_out to any of the reference grid above,
     151              or dom_glo to output on native grid -->
     152        <domain id="dom_out" domain_ref="dom_128_96"/>
     153
     154    </domain_definition>
     155
     156
    145157    <!-- Definition of model variables -->
    146158    <field_definition src="./field_def_physics.xml" />
Note: See TracChangeset for help on using the changeset viewer.