Changeset 5107 for trunk/NEMOGCM


Ignore:
Timestamp:
2015-02-26T18:18:47+01:00 (6 years ago)
Author:
smasson
Message:

update IO for VVL compatibility, see ticket #1474

Location:
trunk/NEMOGCM
Files:
5 deleted
12 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMOGCM/CONFIG/AMM12/EXP00/iodef.xml

    r4695 r5107  
    4040 
    4141   <file id="file2" name_suffix="_grid_U" description="ocean U grid variables" > 
    42      <field field_ref="suoce"        name="uos"     long_name="sea_surface_x_velocity"    /> 
     42     <field field_ref="ssu"          name="uos"     long_name="sea_surface_x_velocity"    /> 
    4343   </file> 
    4444    
    4545   <file id="file3" name_suffix="_grid_V" description="ocean V grid variables" > 
    46      <field field_ref="svoce"        name="vos"     long_name="sea_surface_y_velocity"    /> 
     46     <field field_ref="ssv"          name="vos"     long_name="sea_surface_y_velocity"    /> 
    4747   </file> 
    4848    
     
    7070   <file id="file5" name_suffix="_grid_U" description="ocean U grid variables" > 
    7171     <field field_ref="uoce"         name="uo"      long_name="sea_water_x_velocity"      /> 
    72      <field field_ref="suoce"        name="uos"     long_name="sea_surface_x_velocity"    /> 
     72     <field field_ref="ssu"          name="uos"     long_name="sea_surface_x_velocity"    /> 
    7373     <field field_ref="utau"         name="tauuo"   long_name="surface_downward_x_stress" /> 
    7474   </file> 
     
    7676   <file id="file6" name_suffix="_grid_V" description="ocean V grid variables" > 
    7777     <field field_ref="voce"         name="vo"      long_name="sea_water_y_velocity"      /> 
    78      <field field_ref="svoce"        name="vos"     long_name="sea_surface_y_velocity"    /> 
     78     <field field_ref="ssv"          name="vos"     long_name="sea_surface_y_velocity"    /> 
    7979     <field field_ref="vtau"         name="tauvo"   long_name="surface_downward_y_stress" /> 
    8080   </file> 
  • trunk/NEMOGCM/CONFIG/ORCA2_LIM/EXP00/iodef.xml

    r4696 r5107  
    3535  
    3636   <file id="file1" name_suffix="_grid_T" description="ocean T grid variables" > 
     37     <field field_ref="e3t"  /> 
     38     <field field_ref="toce"         name="thetao"   long_name="sea_water_potential_temperature"                    operation="instant" freq_op="5d" > @toce_e3t / @e3t </field> 
     39     <field field_ref="soce"         name="so"       long_name="sea_water_salinity"                                 operation="instant" freq_op="5d" > @soce_e3t / @e3t </field> 
    3740     <field field_ref="sst"          name="tos"      long_name="sea_surface_temperature"                       /> 
    3841     <field field_ref="sss"          name="sos"      long_name="sea_surface_salinity"                          /> 
    3942     <field field_ref="ssh"          name="zos"      long_name="sea_surface_height_above_geoid"                /> 
    40      <field field_ref="toce"         name="thetao"   long_name="sea_water_potential_temperature"               /> 
    41      <field field_ref="soce"         name="so"       long_name="sea_water_salinity"                            /> 
    42      <field field_ref="sst2"         name="tossq"    long_name="square_of_sea_surface_temperature"             /> 
    43      <field field_ref="ssh2"         name="zossq"    long_name="square_of_sea_surface_height_above_geoid"      /> 
     43     <field field_ref="sst"          name="tosstd"   long_name="sea_surface_temperature_standard_deviation"         operation="instant" freq_op="5d" > sqrt( @sst2 - @sst * @sst ) </field> 
     44     <field field_ref="ssh"          name="zosstd"   long_name="sea_surface_height_above_geoid_standard_deviation"  operation="instant" freq_op="5d" > sqrt( @ssh2 - @ssh * @ssh ) </field> 
     45     <field field_ref="sst"          name="sstdcy"   long_name="amplitude of sst diurnal cycle" operation="average" freq_op="1d" > @sstmax - @sstmin </field> 
    4446     <field field_ref="mldkz5"       /> 
    4547     <field field_ref="mldr10_1"     /> 
     48     <field field_ref="mldr10_1"     name="mldr10_1dcy"  long_name="amplitude of mldr10_1 diurnal cycle" operation="average" freq_op="1d" > @mldr10_1max - @mldr10_1min </field> 
    4649     <field field_ref="empmr"        name="wfo"      long_name="water_flux_into_sea_water"                     /> 
    4750     <field field_ref="qsr"          name="rsntds"   long_name="surface_net_downward_shortwave_flux"           /> 
    4851     <field field_ref="qt"           name="tohfls"   long_name="surface_net_downward_total_heat_flux"          /> 
    49      <field field_ref="saltflx"      name="sosflxdo"  /> 
    50      <field field_ref="taum"         name="taum" /> 
    51      <field field_ref="wspd"         name="sowindsp"  /> 
     52     <field field_ref="saltflx"      name="sosflxdo" /> 
     53     <field field_ref="taum"         name="taum"     /> 
     54     <field field_ref="wspd"         name="sowindsp" /> 
    5255          <field field_ref="precip"       name="soprecip" /> 
     56     <field field_ref="sbt"                          /> 
    5357   </file> 
    5458 
    5559   <file id="file3" name_suffix="_grid_U" description="ocean U grid variables" > 
    56      <field field_ref="suoce"        name="uos"     long_name="sea_surface_x_velocity"    /> 
    57      <field field_ref="uoce"         name="uo"      long_name="sea_water_x_velocity"      /> 
     60     <field field_ref="e3u"  /> 
     61     <field field_ref="ssu"          name="uos"     long_name="sea_surface_x_velocity"    /> 
     62     <field field_ref="uoce"         name="uo"      long_name="sea_water_x_velocity" operation="instant" freq_op="5d" > @uoce_e3u / @e3u </field> 
    5863     <field field_ref="utau"         name="tauuo"   long_name="surface_downward_x_stress" /> 
    5964          <!-- variables available with MLE 
     
    6368    
    6469   <file id="file4" name_suffix="_grid_V" description="ocean V grid variables" > 
    65      <field field_ref="svoce"        name="vos"     long_name="sea_surface_y_velocity"    /> 
    66      <field field_ref="voce"         name="vo"      long_name="sea_water_y_velocity"      /> 
     70     <field field_ref="e3v"  /> 
     71     <field field_ref="ssv"          name="vos"     long_name="sea_surface_y_velocity"    /> 
     72     <field field_ref="voce"         name="vo"      long_name="sea_water_y_velocity" operation="instant" freq_op="5d" > @voce_e3v / @e3v </field> 
    6773     <field field_ref="vtau"         name="tauvo"   long_name="surface_downward_y_stress" /> 
    6874          <!-- variables available with MLE 
     
    7278    
    7379   <file id="file5" name_suffix="_grid_W" description="ocean W grid variables" > 
    74      <field field_ref="woce"         name="wo"      long_name="ocean vertical velocity"         /> 
     80     <field field_ref="e3w"  /> 
     81     <field field_ref="woce"         name="wo"      long_name="ocean vertical velocity" operation="instant" freq_op="5d" > @woce_e3w / @e3w </field> 
    7582     <field field_ref="avt"          name="difvho"  long_name="ocean_vertical_heat_diffusivity" /> 
    7683   </file> 
     
    156163        We must have buffer_size > jpi*jpj*jpk*8 (with jpi and jpj the subdomain size) 
    157164--> 
    158      <variable id="buffer_size"               type="integer">5000000</variable> 
     165     <variable id="buffer_size"               type="integer">50000000</variable> 
    159166     <variable id="buffer_server_factor_size" type="integer">2</variable> 
    160167     <variable id="info_level"                type="integer">0</variable> 
  • trunk/NEMOGCM/CONFIG/ORCA2_LIM/EXP00/iodef_ar5.xml

    r4153 r5107  
    4949     <field field_ref="ssh"          name="zos"      long_name="sea_surface_height_above_geoid"                                 /> 
    5050     <field field_ref="ssh2"         name="zossq"    long_name="square_of_sea_surface_height_above_geoid"             level="2" /> 
    51      <!-- masscello : sea_water_mass_per_unit_area = cellthc*rau0                                          no time changes --> 
    52      <field field_ref="cellthc"      name="thkcello" long_name="cell_thickness"                                                 /> <!-- no time changes --> 
     51     <!-- masscello : sea_water_mass_per_unit_area = e3t*rau0                                          no time changes --> 
     52     <field field_ref="e3t"          name="thkcello" long_name="cell_thickness"                                                 /> <!-- no time changes --> 
    5353     <field field_ref="toce"         name="thetao"   long_name="sea_water_potential_temperature"                                /> 
    5454     <field field_ref="sst"          name="tos"      long_name="sea_surface_temperature"                              level="1" /> 
  • trunk/NEMOGCM/CONFIG/ORCA2_LIM/EXP00/iodef_crs.xml

    r4152 r5107  
    4040 
    4141   <file id="file2" name_suffix="_grid_U" description="ocean U grid variables" > 
    42      <field field_ref="suoce"        name="uos"     long_name="sea_surface_x_velocity"    /> 
     42     <field field_ref="ssu"          name="uos"     long_name="sea_surface_x_velocity"    /> 
    4343   </file> 
    4444    
    4545   <file id="file3" name_suffix="_grid_V" description="ocean V grid variables" > 
    46      <field field_ref="svoce"        name="vos"     long_name="sea_surface_y_velocity"    /> 
     46     <field field_ref="ssv"          name="vos"     long_name="sea_surface_y_velocity"    /> 
    4747   </file> 
    4848    
     
    7171   <file id="file5" name_suffix="_grid_U" description="ocean U grid variables" > 
    7272     <field field_ref="uoce"         name="uo"      long_name="sea_water_x_velocity"      /> 
    73      <field field_ref="suoce"        name="uos"     long_name="sea_surface_x_velocity"    /> 
     73     <field field_ref="ssu"          name="uos"     long_name="sea_surface_x_velocity"    /> 
    7474     <field field_ref="utau"         name="tauuo"   long_name="surface_downward_x_stress" /> 
    7575   </file> 
     
    7777   <file id="file6" name_suffix="_grid_V" description="ocean V grid variables" > 
    7878     <field field_ref="voce"         name="vo"      long_name="sea_water_y_velocity"      /> 
    79      <field field_ref="svoce"        name="vos"     long_name="sea_surface_y_velocity"    /> 
     79     <field field_ref="ssv"          name="vos"     long_name="sea_surface_y_velocity"    /> 
    8080     <field field_ref="vtau"         name="tauvo"   long_name="surface_downward_y_stress" /> 
    8181   </file> 
  • trunk/NEMOGCM/CONFIG/ORCA2_LIM3/EXP00/iodef.xml

    r4990 r5107  
    3838 
    3939   <file id="file1" name_suffix="_grid_T" description="ocean T grid variables" > 
     40     <field field_ref="e3t"  /> 
     41     <field field_ref="toce"         name="thetao"   long_name="sea_water_potential_temperature"                    operation="instant" freq_op="5d" > @toce_e3t / @e3t </field> 
     42     <field field_ref="soce"         name="so"       long_name="sea_water_salinity"                                 operation="instant" freq_op="5d" > @soce_e3t / @e3t </field> 
    4043     <field field_ref="sst"          name="tos"      long_name="sea_surface_temperature"                       /> 
    4144     <field field_ref="sss"          name="sos"      long_name="sea_surface_salinity"                          /> 
    4245     <field field_ref="ssh"          name="zos"      long_name="sea_surface_height_above_geoid"                /> 
    43      <field field_ref="toce"         name="thetao"   long_name="sea_water_potential_temperature"               /> 
    44      <field field_ref="soce"         name="so"       long_name="sea_water_salinity"                            /> 
    45      <field field_ref="sst2"         name="tossq"    long_name="square_of_sea_surface_temperature"             /> 
    46      <field field_ref="ssh2"         name="zossq"    long_name="square_of_sea_surface_height_above_geoid"      /> 
     46     <field field_ref="sst"          name="tosstd"   long_name="sea_surface_temperature_standard_deviation"         operation="instant" freq_op="5d" > sqrt( @sst2 - @sst * @sst ) </field> 
     47     <field field_ref="ssh"          name="zosstd"   long_name="sea_surface_height_above_geoid_standard_deviation"  operation="instant" freq_op="5d" > sqrt( @ssh2 - @ssh * @ssh ) </field> 
     48     <field field_ref="sst"          name="sstdcy"   long_name="amplitude of sst diurnal cycle" operation="average" freq_op="1d" > @sstmax - @sstmin </field> 
    4749     <field field_ref="mldkz5"       /> 
    4850     <field field_ref="mldr10_1"     /> 
    49       <field field_ref="heatc"        name="heatc"    long_name="Heat content vertically integrated"            /> 
    50       <field field_ref="saltc"        name="saltc"    long_name="Salt content vertically integrated"            /> 
     51     <field field_ref="mldr10_1"     name="mldr10_1dcy"  long_name="amplitude of mldr10_1 diurnal cycle" operation="average" freq_op="1d" > @mldr10_1max - @mldr10_1min </field> 
     52     <field field_ref="sbt"                          /> 
     53     <field field_ref="heatc"        name="heatc"    long_name="Heat content vertically integrated"            /> 
     54     <field field_ref="saltc"        name="saltc"    long_name="Salt content vertically integrated"            /> 
    5155   </file> 
    5256 
     
    7276 
    7377   <file id="file3" name_suffix="_grid_U" description="ocean U grid variables" > 
     78     <field field_ref="e3u"  /> 
    7479     <field field_ref="ssu"          name="uos"     long_name="sea_surface_x_velocity"    /> 
    75      <field field_ref="uoce"         name="uo"      long_name="sea_water_x_velocity"      /> 
     80     <field field_ref="uoce"         name="uo"      long_name="sea_water_x_velocity" operation="instant" freq_op="5d" > @uoce_e3u / @e3u </field> 
    7681     <field field_ref="utau"         name="tauuo"   long_name="surface_downward_x_stress" /> 
    7782          <!-- available with key_diaar5 --> 
     
    8287    
    8388   <file id="file4" name_suffix="_grid_V" description="ocean V grid variables" > 
     89     <field field_ref="e3v"  /> 
    8490     <field field_ref="ssv"          name="vos"     long_name="sea_surface_y_velocity"    /> 
    85      <field field_ref="voce"         name="vo"      long_name="sea_water_y_velocity"      /> 
     91     <field field_ref="voce"         name="vo"      long_name="sea_water_y_velocity" operation="instant" freq_op="5d" > @voce_e3v / @e3v </field> 
    8692     <field field_ref="vtau"         name="tauvo"   long_name="surface_downward_y_stress" /> 
    8793          <!-- available with key_diaar5 --> 
     
    9298    
    9399   <file id="file5" name_suffix="_grid_W" description="ocean W grid variables" > 
    94      <field field_ref="woce"         name="wo"      long_name="ocean vertical velocity"         /> 
     100     <field field_ref="e3w"  /> 
     101     <field field_ref="woce"         name="wo"      long_name="ocean vertical velocity" operation="instant" freq_op="5d" > @woce_e3w / @e3w </field> 
    95102     <field field_ref="avt"          name="difvho"  long_name="ocean_vertical_heat_diffusivity" /> 
    96103     <field field_ref="w_masstr"     name="vovematr"  /> 
     
    323330        We must have buffer_size > jpi*jpj*jpk*8 (with jpi and jpj the subdomain size) 
    324331--> 
    325      <variable id="buffer_size"               type="integer">5000000</variable> 
     332     <variable id="buffer_size"               type="integer">50000000</variable> 
    326333     <variable id="buffer_server_factor_size" type="integer">2</variable> 
    327334     <variable id="info_level"                type="integer">0</variable> 
  • trunk/NEMOGCM/CONFIG/ORCA2_LIM_CFC_C14b/EXP00/iodef.xml

    r4153 r5107  
    5353   <file id="file2" name_suffix="_grid_U" description="ocean U grid variables" > 
    5454     <field field_ref="uoce"         name="uo"      long_name="sea_water_x_velocity"      /> 
    55      <field field_ref="suoce"        name="uos"     long_name="sea_surface_x_velocity"    /> 
     55     <field field_ref="ssu"          name="uos"     long_name="sea_surface_x_velocity"    /> 
    5656     <field field_ref="utau"         name="tauuo"   long_name="surface_downward_x_stress" /> 
    5757   </file> 
     
    5959   <file id="file3" name_suffix="_grid_V" description="ocean V grid variables" > 
    6060     <field field_ref="voce"         name="vo"      long_name="sea_water_y_velocity"      /> 
    61      <field field_ref="svoce"        name="vos"     long_name="sea_surface_y_velocity"    /> 
     61     <field field_ref="ssv"          name="vos"     long_name="sea_surface_y_velocity"    /> 
    6262     <field field_ref="vtau"         name="tauvo"   long_name="surface_downward_y_stress" /> 
    6363   </file> 
  • trunk/NEMOGCM/CONFIG/ORCA2_LIM_PISCES/EXP00/iodef.xml

    r4996 r5107  
    4040 
    4141   <file id="file2" name_suffix="_grid_U" description="ocean U grid variables" > 
    42      <field field_ref="suoce"        name="uos"     long_name="sea_surface_x_velocity"    /> 
     42     <field field_ref="ssu"          name="uos"     long_name="sea_surface_x_velocity"    /> 
    4343   </file> 
    4444    
    4545   <file id="file3" name_suffix="_grid_V" description="ocean V grid variables" > 
    46      <field field_ref="svoce"        name="vos"     long_name="sea_surface_y_velocity"    /> 
     46     <field field_ref="ssv"          name="vos"     long_name="sea_surface_y_velocity"    /> 
    4747   </file> 
    4848 
     
    8686   <file id="file6" name_suffix="_grid_U" description="ocean U grid variables" > 
    8787     <field field_ref="uoce"         name="uo"      long_name="sea_water_x_velocity"      /> 
    88      <field field_ref="suoce"        name="uos"     long_name="sea_surface_x_velocity"    /> 
     88     <field field_ref="ssu"          name="uos"     long_name="sea_surface_x_velocity"    /> 
    8989     <field field_ref="utau"         name="tauuo"   long_name="surface_downward_x_stress" /> 
    9090          <!-- variables available with MLE 
     
    9595   <file id="file7" name_suffix="_grid_V" description="ocean V grid variables" > 
    9696     <field field_ref="voce"         name="vo"      long_name="sea_water_y_velocity"      /> 
    97      <field field_ref="svoce"        name="vos"     long_name="sea_surface_y_velocity"    /> 
     97     <field field_ref="ssv"          name="vos"     long_name="sea_surface_y_velocity"    /> 
    9898     <field field_ref="vtau"         name="tauvo"   long_name="surface_downward_y_stress" /> 
    9999          <!-- variables available with MLE 
  • trunk/NEMOGCM/CONFIG/ORCA2_SAS_LIM/EXP00/iodef.xml

    r4696 r5107  
    4343 
    4444   <file id="file2" name_suffix="_grid_U" description="ocean U grid variables" > 
    45      <field field_ref="suoce"        name="suoce"     long_name="sea_surface_x_velocity"    /> 
     45     <field field_ref="ssu"          name="suoce"     long_name="sea_surface_x_velocity"    /> 
    4646   </file> 
    4747    
    4848   <file id="file3" name_suffix="_grid_V" description="ocean V grid variables" > 
    49      <field field_ref="svoce"        name="svoce"     long_name="sea_surface_y_velocity"    /> 
     49     <field field_ref="ssv"          name="svoce"     long_name="sea_surface_y_velocity"    /> 
    5050   </file> 
    5151 
     
    6969    
    7070   <file id="file5" name_suffix="_grid_U" description="ocean U grid variables" > 
    71      <field field_ref="suoce"        name="uos"     long_name="sea_surface_x_velocity"    /> 
     71     <field field_ref="ssu"          name="uos"     long_name="sea_surface_x_velocity"    /> 
    7272     <field field_ref="utau"         name="tauuo"   long_name="surface_downward_x_stress" /> 
    7373   </file> 
    7474    
    7575   <file id="file6" name_suffix="_grid_V" description="ocean V grid variables" > 
    76      <field field_ref="svoce"        name="vos"     long_name="sea_surface_y_velocity"    /> 
     76     <field field_ref="ssv"          name="vos"     long_name="sea_surface_y_velocity"    /> 
    7777     <field field_ref="vtau"         name="tauvo"   long_name="surface_downward_y_stress" /> 
    7878   </file> 
  • trunk/NEMOGCM/CONFIG/SHARED/field_def.xml

    r4996 r5107  
    1616       
    1717     <field_group id="grid_T" grid_ref="grid_T_2D" > 
    18          <field id="toce"         long_name="temperature"                               unit="degC" grid_ref="grid_T_3D"/> 
    19          <field id="soce"         long_name="salinity"                                  unit="psu"  grid_ref="grid_T_3D"/> 
     18         <field id="e3t"          long_name="T-cell thickness"                          unit="m"      grid_ref="grid_T_3D"/> 
     19 
     20         <field id="toce"         long_name="temperature"                               unit="degC"   grid_ref="grid_T_3D"/> 
     21         <field id="toce_e3t"     long_name="temperature * e3t"                         unit="degC*m" grid_ref="grid_T_3D" > toce * e3t </field > 
     22         <field id="soce"         long_name="salinity"                                  unit="psu"    grid_ref="grid_T_3D"/> 
     23         <field id="soce_e3t"     long_name="salinity * e3t"                            unit="psu*m"  grid_ref="grid_T_3D" > soce * e3t </field > 
     24 
    2025         <field id="sst"          long_name="sea surface temperature"                   unit="degC"                     /> 
    21          <field id="sst2"         long_name="square of sea surface temperature"         unit="degC2"                    /> 
     26         <field id="sst2"         long_name="square of sea surface temperature"         unit="degC2"                     >  sst * sst </field > 
     27         <field id="sstmax"       long_name="max of sea surface temperature"       field_ref="sst"  operation="maximum" /> 
     28         <field id="sstmin"       long_name="min of sea surface temperature"       field_ref="sst"  operation="minimum" /> 
    2229         <field id="sstgrad"      long_name="module of sst gradient"                    unit="degC/m"                   /> 
    2330         <field id="sstgrad2"     long_name="square of module of sst gradient"          unit="degC2/m2"                 /> 
     31         <field id="sbt"          long_name="sea bottom temperature"                    unit="degC"                     /> 
     32     
    2433         <field id="sss"          long_name="sea surface salinity"                      unit="psu"                      /> 
    25          <field id="sss2"         long_name="square of sea surface salinity"            unit="psu2"                     /> 
     34         <field id="sss2"         long_name="square of sea surface salinity"            unit="psu2"                      >  sss * sss </field > 
     35         <field id="sssmax"       long_name="max of sea surface salinity"          field_ref="sss"  operation="maximum" /> 
     36         <field id="sssmin"       long_name="min of sea surface salinity"          field_ref="sss"  operation="minimum" /> 
     37         <field id="sbs"          long_name="sea bottom salinity"                       unit="psu"                      /> 
     38 
    2639         <field id="ssh"          long_name="sea surface height"                        unit="m"                        /> 
    27          <field id="ssh2"         long_name="square of sea surface height"              unit="m2"                       /> 
     40         <field id="ssh2"         long_name="square of sea surface height"              unit="m2"                        >  ssh * ssh </field > 
     41         <field id="sshmax"       long_name="max of sea surface height"            field_ref="ssh"  operation="maximum" /> 
     42 
    2843         <field id="mldkz5"       long_name="mixing layer depth (Turbocline)"           unit="m"                        /> 
    2944         <field id="mldr10_1"     long_name="Mixed Layer Depth 0.01 ref.10m"            unit="m"                        /> 
    30          <field id="heatc"        long_name="Heat content vertically integrated"        unit="J/m2"                      /> 
     45         <field id="mldr10_1max"  long_name="max of Mixed Layer Depth 0.01 ref.10m" field_ref="mldr10_1"  operation="maximum" /> 
     46         <field id="mldr10_1min"  long_name="min of Mixed Layer Depth 0.01 ref.10m" field_ref="mldr10_1"  operation="minimum" /> 
     47         <field id="heatc"        long_name="Heat content vertically integrated"        unit="J/m2"                     /> 
    3148         <field id="saltc"        long_name="Salt content vertically integrated"        unit="PSU*kg/m2"                /> 
    3249         <!-- EOS --> 
     
    5774         <field id="botpres"      long_name="Pressure at sea floor"                     unit="dbar"                     /> 
    5875         <!-- variables available with key_diaar5/key_vvl --> 
    59          <field id="cellthc"      long_name="Cell thickness"                            unit="m"    grid_ref="grid_T_3D"/> 
    6076         <!-- variables available with key_vvl --> 
    6177         <field id="tpt_dep"      long_name="T-point depth"                             unit="m"    grid_ref="grid_T_3D"/> 
     
    145161 
    146162         <field id="empmr"        long_name="Net Upward Water Flux"                                        unit="kg/m2/s"  /> 
    147          <field id="saltflx"      long_name="Downward salt flux"                                           unit="PSU/m2/s"  /> 
     163         <field id="saltflx"      long_name="Downward salt flux"                                           unit="PSU/m2/s" /> 
    148164         <field id="fmmflx"       long_name="Water flux due to freezing/melting"                           unit="kg/m2/s"  /> 
    149165         <field id="snowpre"      long_name="Snow precipitation"                                           unit="kg/m2/s"  /> 
    150166         <field id="runoffs"      long_name="River Runoffs"                                                unit="Kg/m2/s"  /> 
    151          <field id="precip"       long_name="Total precipitation"                                           unit="kg/m2/s"  /> 
     167         <field id="precip"       long_name="Total precipitation"                                          unit="kg/m2/s"  /> 
    152168 
    153169 
     
    169185         <field id="isfgammat"    long_name="transfert coefficient for isf (temperature) "                 unit="m/s"      /> 
    170186         <field id="isfgammas"    long_name="transfert coefficient for isf (salinity)    "                 unit="m/s"      /> 
    171         <field id="stbl"         long_name="salinity in the Losh tbl                    "                 unit="PSU"      /> 
    172         <field id="ttbl"         long_name="temperature in the Losh tbl                 "                 unit="C"        /> 
     187    <field id="stbl"         long_name="salinity in the Losh tbl                    "                 unit="PSU"      /> 
     188    <field id="ttbl"         long_name="temperature in the Losh tbl                 "                 unit="C"        /> 
    173189 
    174190         <!-- *_oce variables available with ln_blk_clio or ln_blk_core --> 
     
    202218         <field id="snowmel_cea"   long_name="Snow Melt Rate (cell average)"                               unit="kg/m2/s"  /> 
    203219         <field id="sntoice_cea"   long_name="Snow-Ice Formation Rate (cell average)"                      unit="kg/m2/s"  /> 
    204          <field id="ticemel_cea"      long_name="Rate of Melt at Upper Surface of Sea Ice (cell average)"     unit="kg/m2/s"  /> 
     220         <field id="ticemel_cea"   long_name="Rate of Melt at Upper Surface of Sea Ice (cell average)"     unit="kg/m2/s"  /> 
    205221 
    206222    <!-- ice fields --> 
     
    259275 
    260276         <field id="micet"        long_name="Mean ice temperature"                                         unit="degC"     /> 
    261          <field id="icehc"        long_name="ice total heat content"                                       unit="10^9J"   />  
     277         <field id="icehc"        long_name="ice total heat content"                                       unit="10^9J"    />  
    262278         <field id="isnowhc"      long_name="snow total heat content"                                      unit="10^9J"    /> 
    263279         <field id="icest"        long_name="ice surface temperature"                                      unit="degC"     /> 
     
    268284         <field id="idive"        long_name="divergence"                                                   unit="10-8s-1"  /> 
    269285         <field id="ishear"       long_name="shear"                                                        unit="10-8s-1"  /> 
    270          <field id="icevolu"      long_name="ice volume"                                                   unit="m"      /> 
    271          <field id="snowvol"      long_name="snow volume"                                                  unit="m"      /> 
     286         <field id="icevolu"      long_name="ice volume"                                                   unit="m"        /> 
     287         <field id="snowvol"      long_name="snow volume"                                                  unit="m"        /> 
    272288 
    273289         <field id="icetrp"       long_name="ice volume transport"                                         unit="m/day"   /> 
    274290         <field id="snwtrp"       long_name="snw volume transport"                                         unit="m/day"   /> 
    275          <field id="deitrp"       long_name="advected ice enhalpy"                                         unit="W/m2"   /> 
    276          <field id="destrp"       long_name="advected snw enhalpy"                                         unit="W/m2"   /> 
     291         <field id="deitrp"       long_name="advected ice enhalpy"                                         unit="W/m2"    /> 
     292         <field id="destrp"       long_name="advected snw enhalpy"                                         unit="W/m2"    /> 
    277293 
    278294         <field id="sfxbri"       long_name="brine salt flux"                                              unit="psu*kg/m2/day" /> 
     
    314330         <field id="hfxres"   long_name="heat fluxes from ice-ocean mass exchange during resultant"           unit="W/m2"  /> 
    315331         <field id="hfxsub"   long_name="heat fluxes from ice-atm. mass exchange during sublimation"          unit="W/m2"  /> 
    316          <field id="hfxspr"   long_name="heat fluxes from ice-atm. mass exchange during snow precip"          unit="W/m2" /> 
     332         <field id="hfxspr"   long_name="heat fluxes from ice-atm. mass exchange during snow precip"          unit="W/m2"  /> 
    317333 
    318334    <!-- diags --> 
    319335         <field id="hfxdhc"    long_name="Heat content variation in snow and ice"   unit="W/m2" /> 
    320          <field id="hfxtur"    long_name="turbulent heat flux at the ice base"      unit="W/m2"  /> 
     336         <field id="hfxtur"    long_name="turbulent heat flux at the ice base"      unit="W/m2" /> 
    321337 
    322338      </field_group> 
     
    325341       
    326342      <field_group id="grid_U"   grid_ref="grid_U_2D"> 
    327          <field id="utau"         long_name="Wind Stress along i-axis"                    unit="N/m2" /> 
    328          <field id="uoce"         long_name="ocean current along i-axis"                  unit="m/s"  grid_ref="grid_U_3D" /> 
    329          <field id="ssu"          long_name="ocean surface current along i-axis"          unit="m/s"     /> 
    330          <field id="uocetr_eff"   long_name="Effective ocean transport along i-axis"      unit="m3/s" grid_ref="grid_U_3D" /> 
     343         <field id="e3u"          long_name="U-cell thickness"                               unit="m"        grid_ref="grid_U_3D" /> 
     344         <field id="utau"         long_name="Wind Stress along i-axis"                       unit="N/m2"                          /> 
     345         <field id="uoce"         long_name="ocean current along i-axis"                     unit="m/s"      grid_ref="grid_U_3D" /> 
     346         <field id="uoce_e3u"     long_name="ocean current along i-axis * e3u"               unit="m2/s"     grid_ref="grid_U_3D"  > uoce * e3u </field> 
     347         <field id="ssu"          long_name="ocean surface current along i-axis"             unit="m/s"                           /> 
     348         <field id="sbu"          long_name="ocean bottom current along i-axis"              unit="m/s"                           /> 
     349         <field id="uocetr_eff"   long_name="Effective ocean transport along i-axis"         unit="m3/s"     grid_ref="grid_U_3D" /> 
    331350         <field id="uocet"        long_name="ocean transport along i-axis times temperature" unit="degC.m/s" grid_ref="grid_U_3D" /> 
    332351         <field id="uoces"        long_name="ocean transport along i-axis times salinity"    unit="psu.m/s"  grid_ref="grid_U_3D" /> 
    333352         <!-- variables available with MLE --> 
    334          <field id="psiu_mle"     long_name="MLE streamfunction along i-axis"             unit="m3/s"  grid_ref="grid_U_3D" /> 
     353         <field id="psiu_mle"     long_name="MLE streamfunction along i-axis"                unit="m3/s"     grid_ref="grid_U_3D" /> 
    335354         <!-- uoce_eiv: available with key_traldf_eiv and key_diaeiv --> 
    336          <field id="uoce_eiv"     long_name="EIV ocean current along i-axis"              unit="m/s"  grid_ref="grid_U_3D" /> 
     355         <field id="uoce_eiv"     long_name="EIV ocean current along i-axis"                 unit="m/s"      grid_ref="grid_U_3D" /> 
    337356         <!-- uoce_eiv: available with key_trabbl --> 
    338          <field id="uoce_bbl"     long_name="BBL ocean current along i-axis"              unit="m/s"  grid_ref="grid_U_3D" /> 
    339         <field id="ahu_bbl"      long_name="BBL diffusive flux along i-axis"             unit="m3/s" /> 
     357         <field id="uoce_bbl"     long_name="BBL ocean current along i-axis"                 unit="m/s"      grid_ref="grid_U_3D" /> 
     358    <field id="ahu_bbl"      long_name="BBL diffusive flux along i-axis"                unit="m3/s"                          /> 
    340359         <!-- variable for ice shelves --> 
    341         <field id="utbl"         long_name="zonal current in the Losh tbl"               unit="m/s"  axis_ref="none" /> 
     360    <field id="utbl"         long_name="zonal current in the Losh tbl"                  unit="m/s"                          /> 
    342361         <!-- variables available with key_diaar5 --> 
    343          <field id="u_masstr"     long_name="ocean eulerian mass transport along i-axis"  unit="kg/s" grid_ref="grid_U_3D" /> 
    344          <field id="u_heattr"     long_name="ocean eulerian heat transport along i-axis"  unit="W"    /> 
    345          <field id="u_salttr"     long_name="ocean eulerian salt transport along i-axis"  unit="PSU*kg/s"    /> 
    346          <field id="ueiv_heattr"  long_name="ocean bolus heat transport along i-axis"     unit="W"    /> 
    347          <field id="udiff_heattr" long_name="ocean diffusion heat transport along i-axis" unit="W"    /> 
     362         <field id="u_masstr"     long_name="ocean eulerian mass transport along i-axis"     unit="kg/s"    grid_ref="grid_U_3D" /> 
     363         <field id="u_heattr"     long_name="ocean eulerian heat transport along i-axis"     unit="W"                             /> 
     364         <field id="u_salttr"     long_name="ocean eulerian salt transport along i-axis"     unit="PSU*kg/s"                      /> 
     365         <field id="ueiv_heattr"  long_name="ocean bolus heat transport along i-axis"        unit="W"                             /> 
     366         <field id="udiff_heattr" long_name="ocean diffusion heat transport along i-axis"    unit="W"                             /> 
    348367     </field_group> 
    349368       
     
    351370       
    352371      <field_group id="grid_V"   grid_ref="grid_V_2D"> 
    353          <field id="vtau"         long_name="Wind Stress along j-axis"                    unit="N/m2" /> 
    354          <field id="voce"         long_name="ocean current along j-axis"                  unit="m/s"  grid_ref="grid_V_3D" /> 
    355          <field id="ssv"          long_name="ocean surface current along j-axis"          unit="m/s"  /> 
    356          <field id="vocetr_eff"   long_name="Effective ocean transport along j-axis"      unit="m3/s" grid_ref="grid_V_3D" /> 
     372         <field id="e3v"          long_name="V-cell thickness"                               unit="m"        grid_ref="grid_V_3D" /> 
     373         <field id="vtau"         long_name="Wind Stress along j-axis"                       unit="N/m2"                          /> 
     374         <field id="voce"         long_name="ocean current along j-axis"                     unit="m/s"      grid_ref="grid_V_3D" /> 
     375         <field id="voce_e3v"     long_name="ocean current along j-axis * e3v"               unit="m2/s"     grid_ref="grid_V_3D"  > voce * e3v </field> 
     376         <field id="ssv"          long_name="ocean surface current along j-axis"             unit="m/s"                           /> 
     377         <field id="sbv"          long_name="ocean bottom current along j-axis"              unit="m/s"                           /> 
     378         <field id="vocetr_eff"   long_name="Effective ocean transport along j-axis"         unit="m3/s"     grid_ref="grid_V_3D" /> 
    357379         <field id="vocet"        long_name="ocean transport along j-axis times temperature" unit="degC.m/s" grid_ref="grid_V_3D" /> 
    358380         <field id="voces"        long_name="ocean transport along j-axis times salinity"    unit="psu.m/s"  grid_ref="grid_V_3D" /> 
    359381         <!-- variables available with MLE --> 
    360          <field id="psiv_mle"     long_name="MLE streamfunction along j-axis"             unit="m3/s"  grid_ref="grid_V_3D" /> 
     382         <field id="psiv_mle"     long_name="MLE streamfunction along j-axis"                unit="m3/s"     grid_ref="grid_V_3D" /> 
    361383         <!-- voce_eiv: available with key_traldf_eiv and key_diaeiv --> 
    362          <field id="voce_eiv"     long_name="EIV ocean current along j-axis"              unit="m/s"  grid_ref="grid_V_3D" /> 
     384         <field id="voce_eiv"     long_name="EIV ocean current along j-axis"                 unit="m/s"      grid_ref="grid_V_3D" /> 
    363385         <!-- voce_eiv: available with key_trabbl --> 
    364          <field id="voce_bbl"     long_name="BBL ocean current along j-axis"              unit="m/s"  grid_ref="grid_V_3D" /> 
    365         <field id="ahv_bbl"      long_name="BBL diffusive flux along j-axis"             unit="m3/s"                   /> 
     386         <field id="voce_bbl"     long_name="BBL ocean current along j-axis"                 unit="m/s"      grid_ref="grid_V_3D" /> 
     387    <field id="ahv_bbl"      long_name="BBL diffusive flux along j-axis"                unit="m3/s"                          /> 
    366388         <!-- variable for ice shelves --> 
    367         <field id="vtbl"         long_name="meridional current in the Losh tbl"        unit="m/s"  axis_ref="none" /> 
     389    <field id="vtbl"         long_name="meridional current in the Losh tbl"             unit="m/s"                          /> 
    368390         <!-- variables available with key_diaar5 --> 
    369          <field id="v_masstr"     long_name="ocean eulerian mass transport along j-axis"  unit="kg/s" grid_ref="grid_V_3D" /> 
    370          <field id="v_heattr"     long_name="ocean eulerian heat transport along j-axis"  unit="W"     /> 
    371          <field id="v_salttr"     long_name="ocean eulerian salt transport along i-axis"  unit="PSU*kg/s"    /> 
    372          <field id="veiv_heattr"  long_name="ocean bolus heat transport along j-axis"     unit="W"     /> 
    373          <field id="vdiff_heattr" long_name="ocean diffusion heat transport along j-axis" unit="W"     /> 
     391         <field id="v_masstr"     long_name="ocean eulerian mass transport along j-axis"     unit="kg/s"    grid_ref="grid_V_3D" /> 
     392         <field id="v_heattr"     long_name="ocean eulerian heat transport along j-axis"     unit="W"                             /> 
     393         <field id="v_salttr"     long_name="ocean eulerian salt transport along i-axis"     unit="PSU*kg/s"                      /> 
     394         <field id="veiv_heattr"  long_name="ocean bolus heat transport along j-axis"        unit="W"                             /> 
     395         <field id="vdiff_heattr" long_name="ocean diffusion heat transport along j-axis"    unit="W"                             /> 
    374396      </field_group> 
    375397       
     
    377399       
    378400      <field_group id="grid_W" grid_ref="grid_W_3D"> 
     401        <field id="e3w"          long_name="W-cell thickness"                            unit="m"   /> 
    379402        <field id="woce"         long_name="ocean vertical velocity"                  unit="m/s"    /> 
     403        <field id="woce_e3w"     long_name="ocean vertical velocity * e3v"            unit="m2/s"    > woce * e3w </field> 
    380404        <field id="wocetr_eff"   long_name="effective ocean vertical transport"       unit="m3/s"   /> 
    381405        <!-- woce_eiv: available with key_traldf_eiv and key_diaeiv --> 
     
    430454      <field_group id="SBC_scalar"  domain_ref="1point" > 
    431455         <!-- available with ln_limdiaout --> 
    432     <field id="ibgvoltot"    long_name="global mean ice volume"                   unit="km3"   /> 
    433     <field id="sbgvoltot"    long_name="global mean snow volume"                  unit="km3"   /> 
    434     <field id="ibgarea"      long_name="global mean ice area"                     unit="km2"   /> 
    435     <field id="ibgsaline"    long_name="global mean ice salinity"                 unit="psu"   /> 
    436     <field id="ibgtemper"    long_name="global mean ice temperature"              unit="degC"   /> 
     456    <field id="ibgvoltot"    long_name="global mean ice volume"                   unit="km3"      /> 
     457    <field id="sbgvoltot"    long_name="global mean snow volume"                  unit="km3"      /> 
     458    <field id="ibgarea"      long_name="global mean ice area"                     unit="km2"      /> 
     459    <field id="ibgsaline"    long_name="global mean ice salinity"                 unit="psu"      /> 
     460    <field id="ibgtemper"    long_name="global mean ice temperature"              unit="degC"     /> 
    437461    <field id="ibgheatco"    long_name="global mean ice heat content"             unit="10^20J"   /> 
    438462    <field id="sbgheatco"    long_name="global mean snow heat content"            unit="10^20J"   /> 
    439     <field id="ibgsaltco"    long_name="global mean ice salt content"             unit="psu*km3"   /> 
     463    <field id="ibgsaltco"    long_name="global mean ice salt content"             unit="psu*km3"  /> 
    440464 
    441465    <field id="ibgvfx"       long_name="global mean volume flux (emp)"               unit="m/day"   /> 
     
    466490 
    467491        <field id="ibghfxthd"   long_name="heat fluxes from ice-ocean exchange during thermo"              unit="W"  /> 
    468         <field id="ibghfxsum"   long_name="heat fluxes causing surface ice melt"              unit="W"  /> 
    469         <field id="ibghfxbom"   long_name="heat fluxes causing bottom ice melt"              unit="W"  /> 
    470         <field id="ibghfxbog"   long_name="heat fluxes causing bottom ice growth"              unit="W"  /> 
    471         <field id="ibghfxdif"   long_name="heat fluxes causing ice temperature change"              unit="W"  /> 
    472         <field id="ibghfxopw"   long_name="heat fluxes causing open water ice formation"              unit="W"  /> 
     492        <field id="ibghfxsum"   long_name="heat fluxes causing surface ice melt"                           unit="W"  /> 
     493        <field id="ibghfxbom"   long_name="heat fluxes causing bottom ice melt"                            unit="W"  /> 
     494        <field id="ibghfxbog"   long_name="heat fluxes causing bottom ice growth"                          unit="W"  /> 
     495        <field id="ibghfxdif"   long_name="heat fluxes causing ice temperature change"                     unit="W"  /> 
     496        <field id="ibghfxopw"   long_name="heat fluxes causing open water ice formation"                   unit="W"  /> 
    473497        <field id="ibghfxdyn"   long_name="heat fluxes from ice-ocean exchange during dynamic"             unit="W"  /> 
    474498        <field id="ibghfxres"   long_name="heat fluxes from ice-ocean exchange during resultant"           unit="W"  /> 
     
    479503 
    480504    <field id="ibgfrcvol"    long_name="global mean forcing volume (emp)"         unit="km3"      /> 
    481     <field id="ibgfrcsfx"    long_name="global mean forcing salt   (sfx)"         unit="psu*km3"   /> 
     505    <field id="ibgfrcsfx"    long_name="global mean forcing salt   (sfx)"         unit="psu*km3"  /> 
    482506    <field id="ibgvolgrm"    long_name="global mean ice growth+melt volume"       unit="km3"      /> 
    483507      </field_group> 
  • trunk/NEMOGCM/NEMO/OPA_SRC/DIA/diaar5.F90

    r4990 r5107  
    8282      CALL wrk_alloc( jpi , jpj , jpk        , zrhd      , zrhop    ) 
    8383      CALL wrk_alloc( jpi , jpj , jpk , jpts , ztsn                 ) 
    84  
    85       CALL iom_put( 'cellthc', fse3t(:,:,:) ) 
    8684 
    8785      zarea_ssh(:,:) = area(:,:) * sshn(:,:) 
  • trunk/NEMOGCM/NEMO/OPA_SRC/DIA/diawri.F90

    r4990 r5107  
    142142      ENDIF 
    143143 
    144       IF( lk_vvl ) THEN 
    145          z3d(:,:,:) = tsn(:,:,:,jp_tem) * fse3t_n(:,:,:) 
    146          CALL iom_put( "toce" , z3d                        )   ! heat content 
     144      IF( .NOT.lk_vvl ) THEN 
     145         CALL iom_put( "e3t" , fse3t_n(:,:,:) ) 
     146         CALL iom_put( "e3u" , fse3u_n(:,:,:) ) 
     147         CALL iom_put( "e3v" , fse3v_n(:,:,:) ) 
     148         CALL iom_put( "e3w" , fse3w_n(:,:,:) ) 
     149      ENDIF 
     150       
     151      CALL iom_put( "toce", tsn(:,:,:,jp_tem) )    ! 3D temperature 
     152      CALL iom_put(  "sst", tsn(:,:,1,jp_tem) )    ! surface temperature 
     153      IF ( iom_use("sbt") ) THEN 
    147154         DO jj = 1, jpj 
    148155            DO ji = 1, jpi 
    149                z2d(ji,jj) = tsn(ji,jj,mikt(ji,jj),jp_tem) * fse3t_n(ji,jj,mikt(ji,jj)) 
    150             END DO 
    151          END DO   
    152          CALL iom_put( "sst"  , z2d(:,:)                 )   ! sea surface heat content       
     156               z2d(ji,jj) = tsn(ji,jj,MAX(mbathy(ji,jj),1),jp_tem) 
     157            END DO 
     158         END DO 
     159         CALL iom_put( "sbt", z2d )                ! bottom temperature 
     160      ENDIF 
     161       
     162      CALL iom_put( "soce", tsn(:,:,:,jp_sal) )    ! 3D salinity 
     163      CALL iom_put(  "sss", tsn(:,:,1,jp_sal) )    ! surface salinity 
     164      IF ( iom_use("sbs") ) THEN 
    153165         DO jj = 1, jpj 
    154166            DO ji = 1, jpi 
    155                z2d(ji,jj) = tsn(ji,jj,mikt(ji,jj),jp_tem)**2 * fse3t_n(ji,jj,mikt(ji,jj)) 
    156             END DO 
    157          END DO   
    158          CALL iom_put( "sst2" , z2d(:,:)      )   ! sea surface content of squared temperature 
    159          z3d(:,:,:) = tsn(:,:,:,jp_sal) * fse3t_n(:,:,:)             
    160          CALL iom_put( "soce" , z3d                        )   ! salinity content 
     167               z2d(ji,jj) = tsn(ji,jj,MAX(mbathy(ji,jj),1),jp_sal) 
     168            END DO 
     169         END DO 
     170         CALL iom_put( "sbs", z2d )                ! bottom salinity 
     171      ENDIF 
     172          
     173      CALL iom_put( "uoce", un(:,:,:)         )    ! 3D i-current 
     174      CALL iom_put(  "ssu", un(:,:,1)         )    ! surface i-current 
     175      IF ( iom_use("sbu") ) THEN 
    161176         DO jj = 1, jpj 
    162177            DO ji = 1, jpi 
    163                z2d(ji,jj) = tsn(ji,jj,mikt(ji,jj),jp_sal) * fse3t_n(ji,jj,mikt(ji,jj)) 
    164             END DO 
    165          END DO   
    166          CALL iom_put( "sss"  , z2d(:,:)                 )   ! sea surface salinity content 
     178               z2d(ji,jj) = un(ji,jj,MAX(mbathy(ji,jj),1)) 
     179            END DO 
     180         END DO 
     181         CALL iom_put( "sbu", z2d )                ! bottom i-current 
     182      ENDIF 
     183       
     184      CALL iom_put( "voce", vn(:,:,:)         )    ! 3D j-current 
     185      CALL iom_put(  "ssv", vn(:,:,1)         )    ! surface j-current 
     186      IF ( iom_use("sbv") ) THEN 
    167187         DO jj = 1, jpj 
    168188            DO ji = 1, jpi 
    169                z2d(ji,jj) = tsn(ji,jj,mikt(ji,jj),jp_sal)**2 * fse3t_n(ji,jj,mikt(ji,jj)) 
    170             END DO 
    171          END DO   
    172          CALL iom_put( "sss2" , z2d(:,:)                 )   ! sea surface content of squared salinity 
    173       ELSE 
    174          CALL iom_put( "toce" , tsn(:,:,:,jp_tem)        )   ! temperature 
    175          IF ( iom_use("sst") ) THEN 
    176             DO jj = 1, jpj 
    177                DO ji = 1, jpi 
    178                   z2d(ji,jj) = tsn(ji,jj,mikt(ji,jj),jp_tem) 
    179                END DO 
    180             END DO 
    181             CALL iom_put( "sst"  , z2d(:,:)            ) ! sea surface temperature 
    182          ENDIF 
    183          IF ( iom_use("sst2") )   CALL iom_put( "sst2" , z2d(:,:) * z2d(:,:) ) ! square of sea surface temperature 
    184          CALL iom_put( "soce" , tsn(:,:,:,jp_sal)          )   ! salinity 
    185          IF ( iom_use("sss") ) THEN 
    186             DO jj = 1, jpj 
    187                DO ji = 1, jpi 
    188                   z2d(ji,jj) = tsn(ji,jj,mikt(ji,jj),jp_sal) 
    189                END DO 
    190             END DO 
    191             CALL iom_put( "sss"  , z2d(:,:)            ) ! sea surface salinity 
    192          ENDIF 
    193          CALL iom_put( "sss2" , z2d(:,:) * z2d(:,:) ) ! square of sea surface salinity 
    194       END IF 
    195       IF( lk_vvl .AND. (.NOT. ln_dynadv_vec) ) THEN 
    196          CALL iom_put( "uoce" , umask(:,:,:) * un(:,:,:) * fse3u_n(:,:,:) )    ! i-transport 
    197          CALL iom_put( "voce" , vmask(:,:,:) * vn(:,:,:) * fse3v_n(:,:,:) )    ! j-transport 
    198       ELSE 
    199          CALL iom_put( "uoce" , umask(:,:,:) * un(:,:,:)                  )    ! i-current 
    200          CALL iom_put( "voce" , vmask(:,:,:) * vn(:,:,:)                  )    ! j-current 
    201          IF ( iom_use("ssu") ) THEN 
    202             DO jj = 1, jpj 
    203                DO ji = 1, jpi 
    204                   z2d(ji,jj) = un(ji,jj,miku(ji,jj)) 
    205                END DO 
    206             END DO 
    207             CALL iom_put( "ssu"   , z2d                                    )    ! i-current 
    208          ENDIF 
    209          IF ( iom_use("ssv") ) THEN 
    210             DO jj = 1, jpj 
    211                DO ji = 1, jpi 
    212                   z2d(ji,jj) = vn(ji,jj,mikv(ji,jj)) 
    213                END DO 
    214             END DO 
    215             CALL iom_put( "ssv"   , z2d                                    )    ! j-current 
    216          ENDIF 
    217       ENDIF 
    218       CALL iom_put(    "avt"  , avt                        )    ! T vert. eddy diff. coef. 
    219       CALL iom_put(    "avm"  , avmu                       )    ! T vert. eddy visc. coef. 
    220       IF( lk_zdfddm ) THEN 
    221          CALL iom_put( "avs" , fsavs(:,:,:)                          )    ! S vert. eddy diff. coef. 
    222       ENDIF 
    223  
    224       IF ( iom_use("sstgrad2") .OR. iom_use("sstgrad2") ) THEN 
     189               z2d(ji,jj) = vn(ji,jj,MAX(mbathy(ji,jj),1)) 
     190            END DO 
     191         END DO 
     192         CALL iom_put( "sbv", z2d )                ! bottom j-current 
     193      ENDIF 
     194 
     195      CALL iom_put( "avt" , avt                        )    ! T vert. eddy diff. coef. 
     196      CALL iom_put( "avm" , avmu                       )    ! T vert. eddy visc. coef. 
     197      CALL iom_put( "avs" , fsavs(:,:,:)               )    ! S vert. eddy diff. coef. (useful only with key_zdfddm) 
     198 
     199      IF ( iom_use("sstgrad") .OR. iom_use("sstgrad2") ) THEN 
    225200         DO jj = 2, jpjm1                                    ! sst gradient 
    226201            DO ji = fs_2, fs_jpim1   ! vector opt. 
     
    234209         CALL lbc_lnk( z2d, 'T', 1. ) 
    235210         CALL iom_put( "sstgrad2",  z2d               )    ! square of module of sst gradient 
    236          !CDIR NOVERRCHK< 
    237211         z2d(:,:) = SQRT( z2d(:,:) ) 
    238212         CALL iom_put( "sstgrad" ,  z2d               )    ! module of sst gradient 
     
    243217         z2d(:,:)  = 0._wp  
    244218         DO jk = 1, jpkm1 
    245             DO jj = 2, jpjm1 
    246                DO ji = fs_2, fs_jpim1   ! vector opt. 
     219            DO jj = 1, jpj 
     220               DO ji = 1, jpi 
    247221                  z2d(ji,jj) = z2d(ji,jj) + fse3t(ji,jj,jk) * tsn(ji,jj,jk,jp_tem) * tmask(ji,jj,jk) 
    248222               END DO 
    249223            END DO 
    250224         END DO 
    251          CALL lbc_lnk( z2d, 'T', 1. ) 
    252225         CALL iom_put( "heatc", (rau0 * rcp) * z2d )    ! vertically integrated heat content (J/m2) 
    253226      ENDIF 
     
    256229         z2d(:,:)  = 0._wp  
    257230         DO jk = 1, jpkm1 
    258             DO jj = 2, jpjm1 
    259                DO ji = fs_2, fs_jpim1   ! vector opt. 
     231            DO jj = 1, jpj 
     232               DO ji = 1, jpi 
    260233                  z2d(ji,jj) = z2d(ji,jj) + fse3t(ji,jj,jk) * tsn(ji,jj,jk,jp_sal) * tmask(ji,jj,jk) 
    261234               END DO 
    262235            END DO 
    263236         END DO 
    264          CALL lbc_lnk( z2d, 'T', 1. ) 
    265237         CALL iom_put( "saltc", rau0 * z2d )   ! vertically integrated salt content (PSU*kg/m2) 
    266238      ENDIF 
  • trunk/NEMOGCM/NEMO/OPA_SRC/DOM/domvvl.F90

    r4998 r5107  
    588588      INTEGER, INTENT( in )               :: kt       ! time step 
    589589      !! * Local declarations 
    590       REAL(wp), POINTER, DIMENSION(:,:,:) :: z_e3t_def 
    591590      INTEGER                             :: ji,jj,jk       ! dummy loop indices 
    592591      !!---------------------------------------------------------------------- 
    593592 
    594593      IF( nn_timing == 1 )  CALL timing_start('dom_vvl_sf_swp') 
    595       ! 
    596       CALL wrk_alloc( jpi, jpj, jpk, z_e3t_def                ) 
    597594      ! 
    598595      IF( kt == nit000 )   THEN 
     
    679676      ! Write outputs 
    680677      ! ============= 
    681       z_e3t_def(:,:,:) = ( ( fse3t_n(:,:,:) - e3t_0(:,:,:) ) / e3t_0(:,:,:) * 100 * tmask(:,:,:) ) ** 2 
    682       CALL iom_put( "cellthc" , fse3t_n  (:,:,:) ) 
     678      CALL iom_put(     "e3t" , fse3t_n  (:,:,:) ) 
     679      CALL iom_put(     "e3u" , fse3u_n  (:,:,:) ) 
     680      CALL iom_put(     "e3v" , fse3v_n  (:,:,:) ) 
     681      CALL iom_put(     "e3w" , fse3w_n  (:,:,:) ) 
    683682      CALL iom_put( "tpt_dep" , fsde3w_n (:,:,:) ) 
    684       CALL iom_put( "e3tdef"  , z_e3t_def(:,:,:) ) 
     683      IF( iom_use("e3tdef") )   & 
     684         CALL iom_put( "e3tdef"  , ( ( fse3t_n(:,:,:) - e3t_0(:,:,:) ) / e3t_0(:,:,:) * 100 * tmask(:,:,:) ) ** 2 ) 
    685685 
    686686      ! write restart file 
    687687      ! ================== 
    688688      IF( lrst_oce ) CALL dom_vvl_rst( kt, 'WRITE' ) 
    689       ! 
    690       CALL wrk_dealloc( jpi, jpj, jpk, z_e3t_def ) 
    691689      ! 
    692690      IF( nn_timing == 1 )  CALL timing_stop('dom_vvl_sf_swp') 
Note: See TracChangeset for help on using the changeset viewer.