New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 5260 for branches/2014/dev_r4650_UKMO10_Tidally_Meaned_Diagnostics/NEMOGCM/CONFIG/ORCA2_LIM – NEMO

Ignore:
Timestamp:
2015-05-12T12:37:15+02:00 (9 years ago)
Author:
deazer
Message:

Merged branch with Trunk at revision 5253.
Checked with SETTE, passes modified iodef.xml for AMM12 experiment

Location:
branches/2014/dev_r4650_UKMO10_Tidally_Meaned_Diagnostics/NEMOGCM/CONFIG/ORCA2_LIM
Files:
1 deleted
16 edited
3 copied

Legend:

Unmodified
Added
Removed

  • branches/2014/dev_r4650_UKMO10_Tidally_Meaned_Diagnostics/NEMOGCM/CONFIG/ORCA2_LIM/EXP00/1_namelist_cfg

    r4329 r5260  
    5555/ 
    5656!----------------------------------------------------------------------- 
     57&namsplit      !   time splitting parameters                            ("key_dynspg_ts") 
     58!----------------------------------------------------------------------- 
     59/ 
     60!----------------------------------------------------------------------- 
     61&namcrs        !   Grid coarsening for dynamics output and/or 
     62               !   passive tracer coarsened online simulations 
     63!----------------------------------------------------------------------- 
     64/ 
     65!----------------------------------------------------------------------- 
    5766&namtsd    !   data : Temperature  & Salinity 
    5867!----------------------------------------------------------------------- 
    … …  
    7281/ 
    7382!----------------------------------------------------------------------- 
    74 &namsbc_core   !   namsbc_core  CORE bulk formulea 
    75 !----------------------------------------------------------------------- 
     83&namsbc_core   !   namsbc_core  CORE bulk formulae 
     84!----------------------------------------------------------------------- 
     85!              !  file name                    ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights                               ! rotation ! land/sea mask ! 
     86!              !                               !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename                              ! pairing  ! filename      ! 
     87   sn_wndi     = 'u_10.15JUNE2009_fill'        ,         6         , 'U_10_MOD',   .false.    , .true. , 'yearly'  , 'weights_core2_agrif_bicubic.nc'      , 'Uwnd'   , '' 
     88   sn_wndj     = 'v_10.15JUNE2009_fill'        ,         6         , 'V_10_MOD',   .false.    , .true. , 'yearly'  , 'weights_core2_agrif_bicubic.nc'      , 'Vwnd'   , '' 
     89   sn_qsr      = 'ncar_rad.15JUNE2009_fill'    ,        24         , 'SWDN_MOD',   .false.    , .true. , 'yearly'  , 'weights_core2_agrif_bilinear.nc'     , ''       , '' 
     90   sn_qlw      = 'ncar_rad.15JUNE2009_fill'    ,        24         , 'LWDN_MOD',   .false.    , .true. , 'yearly'  , 'weights_core2_agrif_bilinear.nc'     , ''       , '' 
     91   sn_tair     = 't_10.15JUNE2009_fill'        ,         6         , 'T_10_MOD',   .false.    , .true. , 'yearly'  , 'weights_core2_agrif_bilinear.nc'     , ''       , '' 
     92   sn_humi     = 'q_10.15JUNE2009_fill'        ,         6         , 'Q_10_MOD',   .false.    , .true. , 'yearly'  , 'weights_core2_agrif_bilinear.nc'     , ''       , '' 
     93   sn_prec     = 'ncar_precip.15JUNE2009_fill' ,        -1         , 'PRC_MOD1',   .false.    , .true. , 'yearly'  , 'weights_core2_agrif_bilinear.nc'     , ''       , '' 
     94   sn_snow     = 'ncar_precip.15JUNE2009_fill' ,        -1         , 'SNOW'    ,   .false.    , .true. , 'yearly'  , 'weights_core2_agrif_bilinear.nc'     , ''       , '' 
     95   sn_tdif     = 'taudif_core'                 ,        24         , 'taudif'  ,   .false.    , .true. , 'yearly'  , 'weights_core2_agrif_bilinear.nc'     , ''       , '' 
     96   cn_dir      = './'      !  root directory for the location of the bulk files 
     97   ln_2m       = .false.   !  air temperature and humidity referenced at 2m (T) instead 10m (F) 
     98   ln_taudif   = .false.   !  HF tau contribution: use "mean of stress module - module of the mean stress" data 
     99   ln_bulk2z   = .false.   !  Air temperature/humidity and wind vectors are referenced at heights rn_zqt and rn_zu 
     100   rn_zqt      = 3.        !  Air temperature and humidity reference height (m) (ln_bulk2z) 
     101   rn_zu       = 4.        !  Wind vector reference height (m)                  (ln_bulk2z) 
     102   rn_pfac     = 1.        !  multiplicative factor for precipitation (total & snow) 
     103   rn_efac     = 1.        !  multiplicative factor for evaporation (0. or 1.) 
     104   rn_vfac     = 0.        !  multiplicative factor for ocean/ice velocity  
     105                           !  in the calculation of the wind stress (0.=absolute winds or 1.=relative winds) 
    76106/ 
    77107!----------------------------------------------------------------------- 
    78108&namtra_qsr    !   penetrative solar radiation 
    79109!----------------------------------------------------------------------- 
     110!              !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! 
     111!              !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! 
     112   sn_chl      ='chlorophyll',        -1         , 'CHLA'    ,   .true.     , .true. , 'yearly'  , 'weights_bilin.nc'       , '' 
    80113/ 
    81114!----------------------------------------------------------------------- 
    … …  
    106139&nameos        !   ocean physical parameters 
    107140!----------------------------------------------------------------------- 
     141/ 
     142!----------------------------------------------------------------------- 
     143&namtra_dmp    !   tracer: T & S newtonian damping 
     144!----------------------------------------------------------------------- 
     145   ln_tradmp   =  .false.   !  add a damping termn (T) or not (F) 
    108146/ 
    109147!----------------------------------------------------------------------- 
    … …  
    183221!----------------------------------------------------------------------- 
    184222/ 
     223!----------------------------------------------------------------------- 
     224&namobs       !  observation usage                                      ('key_diaobs') 
     225!----------------------------------------------------------------------- 
     226/ 
     227!----------------------------------------------------------------------- 
     228&nam_asminc   !   assimilation increments                               ('key_asminc') 
     229!----------------------------------------------------------------------- 
     230/ 

  • branches/2014/dev_r4650_UKMO10_Tidally_Meaned_Diagnostics/NEMOGCM/CONFIG/ORCA2_LIM/EXP00/iodef.xml

    r4272 r5260  
    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> 
    … …  
    8289     <field field_ref="iceprod_cea"  name="sip"     long_name="sea_ice_thickness"        /> 
    8390     <field field_ref="ist_ipa"      /> 
    84      <field field_ref="ioceflxb"     /> 
    8591     <field field_ref="uice_ipa"     /> 
    8692     <field field_ref="vice_ipa"     /> 
    … …  
    157163        We must have buffer_size > jpi*jpj*jpk*8 (with jpi and jpj the subdomain size) 
    158164--> 
    159      <variable id="buffer_size"               type="integer">5000000</variable> 
     165     <variable id="buffer_size"               type="integer">50000000</variable> 
    160166     <variable id="buffer_server_factor_size" type="integer">2</variable> 
    161167     <variable id="info_level"                type="integer">0</variable> 

  • branches/2014/dev_r4650_UKMO10_Tidally_Meaned_Diagnostics/NEMOGCM/CONFIG/ORCA2_LIM/EXP00/iodef_ar5.xml

    r4153 r5260  
    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" /> 

  • branches/2014/dev_r4650_UKMO10_Tidally_Meaned_Diagnostics/NEMOGCM/CONFIG/ORCA2_LIM/EXP00/iodef_crs.xml

    r4152 r5260  
    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> 

  • branches/2014/dev_r4650_UKMO10_Tidally_Meaned_Diagnostics/NEMOGCM/CONFIG/ORCA2_LIM/EXP00/namelist_cfg

    r4370 r5260  
    188188!----------------------------------------------------------------------- 
    189189/ 
     190!----------------------------------------------------------------------- 
     191&namobs       !  observation usage                                      ('key_diaobs') 
     192!----------------------------------------------------------------------- 
     193/ 
     194!----------------------------------------------------------------------- 
     195&nam_asminc   !   assimilation increments                               ('key_asminc') 
     196!----------------------------------------------------------------------- 
     197/ 

  • branches/2014/dev_r4650_UKMO10_Tidally_Meaned_Diagnostics/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/README

    r4394 r5260  
    22# prerequired 
    33#---------------------------------------------------------------------- 
    4 To use these idl tools, you need to download some climatogies and mask files; 
     4First, you need either : 
     5- IDL (version 6.4 or above), see : http://www.exelisvis.com/ProductsServices/IDL.aspx 
     6  In this case, you also need to download SAXO which is a free package of IDL scripts: 
     7  define $PATH_SAXO, the path where you will download SAXO and get it through svn with the following command. 
     8  > PATH_SAXO=... 
     9  > svn checkout http://forge.ipsl.jussieu.fr/saxo/svn/trunk/SRC $PATH_SAXO/SAXO_DIR/SRC 
     10 
     11- or the IDL Virtual Machine which is free to use and does not require a license to run , see : 
     12  http://www.exelisvis.com/Support/HelpArticlesDetail/TabId/219/ArtMID/900/ArticleID/12395/The-IDL-Virtual-Machine.aspx 
     13  the virtual machine requites std_main.sav that is distributed with this README. 
     14 
     15Next, to use these idl tools, you need to download some climatogies and mask files; 
    516that you can find here: http://dodsp.idris.fr/reee512/NEMO_OUT/ORCA2_LIM/ 
    617 
    … …  
    1324# define your std_plot_vardef.sh or std_ts_vardef.sh file 
    1425#---------------------------------------------------------------------- 
    15 Use the examples provided in : 
     26These files are needed to define your PATHs, the experiments and variables names you used in your experiment. 
     27 - std_plot_vardef.sh is used to do all plots based on temporal mean (maps, vertical profiles...). 
     28 - std_ts_vardef.sh is used to do all time-series type of plot. 
     29 
     30To build you own std_plot_vardef.sh or std_ts_vardef.sh file; use the examples provided such as: 
    1631 - std_ts_vardef.sh_example1 or std_ts_vardef.sh_example2 
    1732 - std_plot_vardef.sh_example1 or std_plot_vardef.sh_example2 
    18 to build your own std_plot_vardef.sh or std_ts_vardef.sh file. 
    1933 
    20 This file is needed to define you PATH, the experiments and variables names 
     34Note that if you use the IDL Virtual Machine, the variable SAXO_DIR defined in std_plot_vardef.sh or std_ts_vardef.sh is not used. Any definition will be ok. 
    2135 
    22 you can copy std_ts_vardef.sh_example1 con std_ts_vardef.sh 
    23 and std_plot_vardef.sh_example1 std_plot_vardef.sh 
     36#---------------------------------------------------------------------- 
     37# HOW TO USE  
     38#---------------------------------------------------------------------- 
     39./std_main.sh -plot -pdf 
     40or 
     41./std_main.sh -ts -pdf 
    2442 
    2543 
    2644 
    2745 
    28  
    29  
    30 #---------------------------------------------------------------------- 
    31 # short note on: How to build IDL virtual Machine: 
     46########################################################################################################## 
     47# short note for developers of this package on: 
     48#                   How to build the tarball required for IDL virtual Machine: 
     49########################################################################################################## 
     50# 
     51# we need to recreate std_main.sav as soon as we change IDL programmes files as 
     52# std_main.sav contains all ".pro" files aready compiled to be used with the virtual machine 
    3253# 
    3354. ./std_plot_vardef.sh # or . ./std_ts_vardef.sh  

  • branches/2014/dev_r4650_UKMO10_Tidally_Meaned_Diagnostics/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_IceFrac.pro

    r4559 r5260  
    2121   
    2222  IF keyword_set(arc) THEN BEGIN  
    23     domdef, 20, 380, 50, 90 
     23    domdef, 20, 380, 40, 90 
    2424    map = [90, 0, 0] 
    2525  ENDIF 
    2626  IF keyword_set(ant) THEN BEGIN  
    27     domdef, 20, 380, -90, -50 
     27    domdef, 20, 380, -90, -40 
    2828    map = [-90, 0, 0] 
    2929  ENDIF 

  • branches/2014/dev_r4650_UKMO10_Tidally_Meaned_Diagnostics/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_all.pro

    r4559 r5260  
    172172  Ithi_april_1 = {arr:1./float(cnt) * total(reform(Ithi_1.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Ithi_1.unit} 
    173173  jan = where(mm EQ 1, cnt) 
     174  Ithi_jan_1 = {arr:1./float(cnt) * total(reform(Ithi_1.arr[*, *, temporary(jan)],nxt,nyt,cnt), 3), unit:Ithi_1.unit}  
    174175  sept = where(mm EQ 9, cnt) 
    175   Ithi_jan_1 = {arr:1./float(cnt) * total(reform(Ithi_1.arr[*, *, temporary(jan)],nxt,nyt,cnt), 3), unit:Ithi_1.unit}  
    176176  Ithi_sept_1 = {arr:1./float(cnt) * total(reform(Ithi_1.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Ithi_1.unit}  
    177177  undefine, Ithi_1 
    … …  
    288288     Isal_sept_2 = {arr:-1} 
    289289   ENDELSE 
    290 ; 
     290 
    291291  IvelU_1 = read_ncdf(getenv('VAR1_IvelU'), allrecords = allrec, filename = std_file1_I )  
    292292  IvelV_1 = read_ncdf(getenv('VAR1_IvelV'), allrecords = allrec, filename = std_file1_I )  
    … …  
    335335    Ivelo_sept_2 = {arr:-1} 
    336336  ENDELSE 
    337 ; 
     337 
    338338  Iage_1 = read_ncdf(getenv('VAR1_Iage'), allrecords = allrec, filename = std_file1_I )  
    339339  caldat, time, mm 
    … …  
    396396  cnt = cnt+1   &   blabla = 'Erp salinity damping term' 
    397397  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_erp, ERP1, ERP2, _extra = ex 
    398 ; 
     398;; 
    399399  cnt = cnt+1   &   blabla = 'Evaporation - Precipitation - Runoff term' 
    400400  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_emp, EMP1, EMP2, _extra = ex 
    … …  
    439439  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_mld, MLD1, MLD2, MLD, _extra = ex 
    440440; 
     441  cnt = cnt+1   &   blabla = 'Mixed layer depth ortho plan ARCTIC' 
     442  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_MLD_ortho, MLD1, MLD2, /ARC, _extra = ex 
     443;sf 
     444  cnt = cnt+1   &   blabla = 'Mixed layer depth ortho plan ANTARTIC' 
     445  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_MLD_ortho, MLD1, MLD2, /ANT, _extra = ex 
     446;sf 
    441447  cnt = cnt+1   &   blabla = 'Zonal mean Mixed layer depth' 
    442448  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_ZonMld, MLD1, MLD2, MLD, _extra = ex 
    … …  
    564570  cnt = cnt+1   &   blabla = 'Antartic Ice Velocity: SEPT' 
    565571  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceVel, IvelU_sept_1, IvelU_sept_2, IvelV_sept_1, IvelV_sept_2, Ivelo_sept_1, Ivelo_sept_2, /ANT,  /SEPT, _extra = ex 
    566 ; 
     572 
    567573  cnt = cnt+1   &   blabla = 'Arctic Ice Age: MARCH' 
    568574  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceAge, Iage_march_1, Iage_march_2, Ifra_march_1, Ifra_march_2, /ARC, /MARCH, _extra = ex 

  • branches/2014/dev_r4650_UKMO10_Tidally_Meaned_Diagnostics/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_vardef.sh

    r4562 r5260  
    2121#idl_command=idl71 
    2222# 
    23 PS_DIR=$( pwd )/OR2LPv36_ps_plot 
    24 PDF_DIR=$( pwd )/OR2LPv36_pdf_plot 
     23PS_DIR=$( pwd )/p4H25a50-testht_ps_plot 
     24PDF_DIR=$( pwd )/p4H25a50-testht_pdf_plot 
    2525HTML_DIR=$( pwd )/html_plot 
    2626SAXO_DIR=/Users/sflod/SAXO_DIR 
    2727# 
    28 DIR_DATA=/Users/sflod/idl_PLOTS/DATA_STORE/RUN_CLIMATO/500y_pisces     # path of data in NetCDF format 
     28DIR_DATA=/Users/sflod/idl_PLOTS/DATA_STORE/RUN_CLIMATO/lim3_ada     # path of data in NetCDF format 
    2929DIR_CLIMATO=/Users/sflod/idl_PLOTS/CLIMATOLOGIES    # path of climatological data 
    3030DIR_MASK=/Users/sflod/idl_PLOTS/MASK  # path of mask files (ex: subbasins) 
    … …  
    5454# 
    5555#===================== EXP1 ===================== 
    56 # 23610101_23701231 
    57 FILE1_T=OR2LPv36_23610101_23701231_1Y_grid_T.nc   # exp1 grid T input file 
    58 FILE1_U=OR2LPv36_23610101_23701231_1Y_grid_U.nc   # exp1 grid U input file 
    59 FILE1_V=OR2LPv36_23610101_23701231_1Y_grid_V.nc   # exp1 grid V input file 
    60 FILE1_I=OR2LPv36_23610101_23701231_1M_icemod.nc   # exp1 ice    input file 
     56# 
     57FILE1_T=p4H25a50_20410101_20501231_1Y_grid_T.nc   # exp1 grid T input file 
     58FILE1_U=p4H25a50_20410101_20501231_1Y_grid_U.nc   # exp1 grid U input file 
     59FILE1_V=p4H25a50_20410101_20501231_1Y_grid_V.nc   # exp1 grid V input file 
     60FILE1_I=p4H25a50_20410101_20501231_1M_icemod.nc   # exp1 ice    input file 
    6161VAR1_T=thetao 
    6262VAR1_S=so         
    63 VAR1_QNET=tohfls 
    64 VAR1_ERP=sowafldp       
     63VAR1_QNET=qt 
     64VAR1_ERP=wfcorr 
    6565VAR1_EMP=wfo 
    6666VAR1_MLD=mldr10_1 
    6767VAR1_U=uocetr_eff 
    6868VAR1_V=vocetr_eff 
    69 VAR1_Ithick=iicethic 
    70 VAR1_Ifrac=soicecov 
    71 VAR1_Isnow=isnowthi 
    72 #VAR1_Isal=sisali 
    73 #VAR1_Iage=siages 
    74 VAR1_IvelU=iicevelu 
    75 VAR1_IvelV=iicevelv 
    76 #VAR1_Ivelo=sivelo 
     69VAR1_Ithick=sithic 
     70VAR1_Ifrac=siconc 
     71VAR1_Isnow=snvolu 
     72VAR1_Isal=sisali 
     73VAR1_Iage=siages 
     74VAR1_IvelU=sivelu 
     75VAR1_IvelV=sivelv 
     76VAR1_Ivelo=sivelo 
    7777# 
    7878#===================== EXP2 ===================== 
    7979# 
    80 FILE2_T=OR2LPv36_23610101_23701231_1Y_grid_T.nc   # exp1 grid T input file 
    81 FILE2_U=OR2LPv36_23610101_23701231_1Y_grid_U.nc   # exp1 grid U input file 
    82 FILE2_V=OR2LPv36_23610101_23701231_1Y_grid_V.nc   # exp1 grid V input file 
    83 FILE2_I=OR2LPv36_23610101_23701231_1M_icemod.nc   # exp1 ice    input file 
     80FILE2_T=testht_20410101_20501231_1Y_grid_T.nc   # exp1 grid T input file 
     81FILE2_U=testht_20410101_20501231_1Y_grid_U.nc   # exp1 grid U input file 
     82FILE2_V=testht_20410101_20501231_1Y_grid_V.nc   # exp1 grid V input file 
     83FILE2_I=testht_20410101_20501231_1M_icemod.nc   # exp1 ice    input file 
    8484VAR2_T=thetao 
    8585VAR2_S=so 
    86 VAR2_QNET=tohfls 
    87 VAR2_ERP=sowafldp 
     86VAR2_QNET=qt 
     87VAR2_ERP=wfcorr 
    8888VAR2_EMP=wfo 
    8989VAR2_MLD=mldr10_1 
    9090VAR2_U=uocetr_eff 
    9191VAR2_V=vocetr_eff 
    92 VAR2_Ithick=iicethic 
    93 VAR2_Ifrac=soicecov 
    94 VAR2_Isnow=isnowthi 
    95 #VAR2_Isal=sisali 
    96 #VAR2_Iage=siages 
    97 VAR2_IvelU=iicevelu 
    98 VAR2_IvelV=iicevelv 
    99 #VAR2_Ivelo=sivelo 
     92VAR2_Ithick=sithic 
     93VAR2_Ifrac=siconc 
     94VAR2_Isnow=snvolu 
     95VAR2_Isal=sisali 
     96VAR2_Iage=siages 
     97VAR2_IvelU=sivelu 
     98VAR2_IvelV=sivelv 
     99VAR2_Ivelo=sivelo 
    100100# 
    101101######################### Export Variables ############################### 

  • branches/2014/dev_r4650_UKMO10_Tidally_Meaned_Diagnostics/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_vardef.sh.old_names

    r4559 r5260  
    7474VAR1_IvelU=iicevelu 
    7575VAR1_IvelV=iicevelv 
    76 VAR1_Ivelo=iicevelo 
     76#VAR1_Ivelo=iicevelo 
    7777# 
    7878#===================== EXP2 ===================== 
    … …  
    9999VAR2_IvelU=iicevelu 
    100100VAR2_IvelV=iicevelv 
    101 VAR2_Ivelo=iicevelo 
     101#VAR2_Ivelo=iicevelo 
    102102# 
    103103######################### Export Variables ############################### 

  • branches/2014/dev_r4650_UKMO10_Tidally_Meaned_Diagnostics/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_ICE.pro

    r4497 r5260  
    8888  ICE_S_15 = {arr:ICE_S_15 * 1.e-12, unit : '10^12 m^2'} 
    8989  ; 
    90   title = 'Northern Hemisphere'+'!C'+prefix+' (BLACK) '+d1_d2+'!C'+'OBSERVATION (light blue) '+'!C'+' Global Annual Mean Ice Area (CONTINUOUS) '+'!C'+ 'and Extend minus 15% (DASHED)' 
     90  ;;title = 'Northern Hemisphere'+'!C'+prefix+' (BLACK) '+d1_d2+'!C'+'OBSERVATION (light blue) '+'!C'+' Global Annual Mean Ice Area (CONTINUOUS) '+'!C'+ 'and Extend minus 15% (DASHED)' 
     91  title = 'Northern Hemisphere'+'!C'+prefix+' (BLACK) '+d1_d2+'!C'+'OBSERVATION (light blue) '+'!C'+' Global Annual Mean Ice Area (DASHED) '+'!C'+ 'and Extend minus 15% (CONTINUOUS)' 
    9192  jpt=12 
    9293  time=julday(1,15,1900)+30*lindgen(12) 
    93   pltt, ICE_N, 't', MIN = 4., MAX = 16., /REMPLI, /PORTRAIT, XGRIDSTYLE = 1, DATE_FORMAT = '%M' $ 
     94  pltt, ICE_N, 't', MIN = 4., MAX = 16., /REMPLI, /PORTRAIT, LINESTYLE=2, XGRIDSTYLE = 1, DATE_FORMAT = '%M' $ 
    9495       , COLOR = 000 , small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
    95   pltt, ICE_N_15, 't', /REMPLI, /PORTRAIT, LINESTYLE=2 $       ;;; dashed lines is LINESTYLE=2  $ 
     96  pltt, ICE_N_15, 't', /REMPLI, /PORTRAIT $       ;;; dashed lines is LINESTYLE=2  $ 
    9697        , /ov1d, COLOR = 000, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex  
    97   pltt, vICE_area_NH, 't',  /REMPLI, /PORTRAIT  $  
     98  pltt, vICE_area_NH, 't',  /REMPLI, /PORTRAIT, LINESTYLE=2  $  
    9899        , /ov1d, COLOR = 100, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
    99   pltt, vICE_ext_NH, 't', /REMPLI, /PORTRAIT, LINESTYLE=2 $   ;;; dashed lines is LINESTYLE=2  $ 
     100  pltt, vICE_ext_NH, 't', /REMPLI, /PORTRAIT $   ;;; dashed lines is LINESTYLE=2  $ 
    100101        , /ov1d, COLOR = 100, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
    101102; 
    102   title = 'Southern Hemisphere'+'!C'+prefix+' (BLACK) '+d1_d2+'!C'+'OBSERVATION (light blue) '+'!C'+' Global Annual Mean Ice Area (CONTINUOUS) '+'!C'+ 'and Extend minus 15% (DASHED)' 
    103   pltt, ICE_S, 't', MIN = 0., MAX = 20., /REMPLI, /NOERASE , XGRIDSTYLE = 1 , DATE_FORMAT = '%M' $ 
     103  title = 'Southern Hemisphere'+'!C'+prefix+' (BLACK) '+d1_d2+'!C'+'OBSERVATION (light blue) '+'!C'+' Global Annual Mean Ice Area (DASHED) '+'!C'+ 'and Extend minus 15% (CONTINUOUS)' 
     104  pltt, ICE_S, 't', MIN = 0., MAX = 20., /REMPLI, LINESTYLE=2, /NOERASE , XGRIDSTYLE = 1 , DATE_FORMAT = '%M' $ 
    104105        ,COLOR = 000, small = [1, 2, 2], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
    105   pltt, ICE_S_15, 't', /REMPLI, /PORTRAIT, LINESTYLE=2 $  
     106  pltt, ICE_S_15, 't', /REMPLI, /PORTRAIT $  
    106107        , /ov1d, COLOR = 000, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex  
    107   pltt,  vICE_area_SH, 't', /REMPLI, /PORTRAIT  $  
     108  pltt,  vICE_area_SH, 't', /REMPLI, /PORTRAIT, LINESTYLE=2  $  
    108109         , /ov1d, COLOR = 100, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
    109   pltt,  vICE_ext_SH, 't', /REMPLI, /PORTRAIT, LINESTYLE=2 $  
     110  pltt,  vICE_ext_SH, 't', /REMPLI, /PORTRAIT $  
    110111        , /ov1d, COLOR = 100, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
    111112; 
    … …  
    159160 
    160161 
    161   title = 'Northern Hemisphere'+'!C'+prefix+' (BLACK) - '+prefix2+' (RED) '+d1_d2_2+'!C'+'OBSERVATION (light blue) '+'!C'+' Global Annual Mean Ice Area (CONTINUOUS) '+'!C'+ 'and Extend minus 15% (DASHED)' 
     162  ;;title = 'Northern Hemisphere'+'!C'+prefix+' (BLACK) - '+prefix2+' (RED) '+d1_d2_2+'!C'+'OBSERVATION (light blue) '+'!C'+' Global Annual Mean Ice Area (CONTINUOUS) '+'!C'+ 'and Extend minus 15% (DASHED)' 
     163  title = 'Northern Hemisphere'+'!C'+prefix+' (BLACK) - '+prefix2+' (RED) '+d1_d2_2+'!C'+'OBSERVATION (light blue) '+'!C'+' Global Annual Mean Ice Area (DASHED) '+'!C'+ 'and Extend minus 15% (CONTINUOUS)' 
    162164  jpt=12 
    163165  time=julday(1,15,1900)+30*lindgen(12) 
    164   pltt, ICE_N, 't', MIN = 4, MAX = 16,  /REMPLI, /PORTRAIT, XGRIDSTYLE = 1, window = 2, DATE_FORMAT = '%M' $ 
     166  pltt, ICE_N, 't', MIN = 4, MAX = 16,  /REMPLI, /PORTRAIT, LINESTYLE=2, XGRIDSTYLE = 1, window = 2, DATE_FORMAT = '%M' $ 
    165167        , COLOR = 000, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex ; BLACK 
    166   pltt, ICE_N2, 't', /REMPLI, /PORTRAIT $ 
     168  pltt, ICE_N2, 't', /REMPLI, /PORTRAIT , LINESTYLE=2 $ 
    167169        , /ov1d, COLOR = 250, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex  ; RED 
    168   pltt, ICE_N_15, 't', /REMPLI, /PORTRAIT, LINESTYLE=2 $ ; linee tratteggiate LINESTYLE=2  $ 
     170  pltt, ICE_N_15, 't', /REMPLI, /PORTRAIT $ ; linee tratteggiate LINESTYLE=2  $ 
    169171        , /ov1d, COLOR = 000, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex  
    170   pltt, ICE_N2_15, 't', /REMPLI, /PORTRAIT, LINESTYLE=2 $ ; linee tratteggiate LINESTYLE=2  $ 
     172  pltt, ICE_N2_15, 't', /REMPLI, /PORTRAIT $ ; linee tratteggiate LINESTYLE=2  $ 
    171173        , /ov1d, COLOR = 250, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
    172   pltt, vICE_area_NH, 't', /REMPLI, /PORTRAIT  $  
     174  pltt, vICE_area_NH, 't', /REMPLI, /PORTRAIT, LINESTYLE=2  $  
    173175        , /ov1d, COLOR = 100, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex  ; light blue 
    174   pltt, vICE_ext_NH, 't', /REMPLI, /PORTRAIT, LINESTYLE=2 $  
     176  pltt, vICE_ext_NH, 't', /REMPLI, /PORTRAIT $  
    175177        , /ov1d, COLOR = 100, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex  ; blu scuro 
    176178; 
    177   title ='Southern Hemisphere'+'!C'+prefix+' (BLACK) - '+prefix2+' (RED) '+d1_d2_2+'!C'+'OBSERVATION (light blue) '+'!C'+'Global Annual Mean Ice Area (CONTINUS)'+'!C'+ 'and Extend minus 15% (DASHED)' 
     179  title ='Southern Hemisphere'+'!C'+prefix+' (BLACK) - '+prefix2+' (RED) '+d1_d2_2+'!C'+'OBSERVATION (light blue) '+'!C'+'Global Annual Mean Ice Area (DASHED)'+'!C'+ 'and Extend minus 15% (CONTINUOUS)' 
    178180;  title ='Southern Hemisphere'+'!C' 
    179   pltt, ICE_S, 't', MIN = 0., MAX = 20.,  /REMPLI, /NOERASE, XGRIDSTYLE = 1, DATE_FORMAT = '%M' $ 
     181  pltt, ICE_S, 't', MIN = 0., MAX = 20., /REMPLI, LINESTYLE=2, /NOERASE, XGRIDSTYLE = 1, DATE_FORMAT = '%M' $ 
    180182         , COLOR = 000, small = [1, 2, 2], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
    181   pltt, ICE_S2, 't', /REMPLI, /NOERASE $ 
     183  pltt, ICE_S2, 't', /REMPLI, /NOERASE, LINESTYLE=2 $ 
    182184        , /ov1d, COLOR = 250, small = [1, 2, 2], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
    183   pltt, ICE_S_15 , 't', /REMPLI, /PORTRAIT, LINESTYLE=2 $ ; linee tratteggiate LINESTYLE=2  $ 
     185  pltt, ICE_S_15 , 't', /REMPLI, /PORTRAIT $ ; linee tratteggiate LINESTYLE=2  $ 
    184186        , /ov1d, COLOR = 000, small = [1, 2, 2], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
    185   pltt, ICE_S2_15, 't', /REMPLI, /PORTRAIT, LINESTYLE=2 $ ; linee tratteggiate LINESTYLE=2  $ 
     187  pltt, ICE_S2_15, 't', /REMPLI, /PORTRAIT $ ; linee tratteggiate LINESTYLE=2  $ 
    186188        , /ov1d, COLOR = 250, small = [1, 2, 2], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
    187   pltt,  vICE_area_SH, 't', /REMPLI, /PORTRAIT $  
     189  pltt,  vICE_area_SH, 't', /REMPLI, /PORTRAIT, LINESTYLE=2 $  
    188190          , /ov1d, COLOR = 100, small = [1, 2, 2], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
    189   pltt,  vICE_ext_SH, 't', /REMPLI, /PORTRAIT, LINESTYLE=2 $  
     191  pltt,  vICE_ext_SH, 't', /REMPLI, /PORTRAIT $  
    190192        , /ov1d, COLOR = 100, small = [1, 2, 2], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
    191193; 

  • branches/2014/dev_r4650_UKMO10_Tidally_Meaned_Diagnostics/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_ICE_FRAM.pro

    r4514 r5260  
    55arr2d = arr2d[t1:t2] 
    66arr2d = reform(arr2d,12,nyear) ; put in 2D array 
    7 arr2d = total(arr2d,2)/nyear ; total over 2th dimension (i.e.years) 
     7;arr2d = total(arr2d,2)/nyear ; total over 2th dimension (i.e.years) 
     8arr2d = arr2d[*, nyear-1] ;  select last year 
    89 
    910return, arr2d 
    … …  
    3435; 
    3536  d1_d2 = '('+strtrim(date1, 1)+' - '+strtrim(date2, 1)+')' 
     37  d2 = '('+strtrim(date2, 1)+')' 
    3638; 
    3739  iodir = std_iodir_data 
    … …  
    7476  ;; AREA 
    7577  ICE_area_export = reform(ICE_area_export, 12, nyr) 
    76   ICE_area_export = total(ICE_area_export,2)/nyr 
     78  ;ICE_area_export = total(ICE_area_export,2)/nyr  ; old version monthly mean average  over all years 
     79  ICE_area_export = ICE_area_export[*,nyr-1]  ; dim= 12, index 0-11; last year choosen 
    7780  ; ICE_area_export = {arr:ICE_area_export * 1.e-12 * 86400 * 365 , unit : '10^6 Km^2/year'}  ; annual mean 
    7881  ICE_area_export = {arr:ICE_area_export * 1.e-12 * 86400 * 30 , unit : '10^6 Km^2/month'}    ; monthly mean 
    … …  
    8083  ; 
    8184  ICE_vol_export = reform(ICE_vol_export, 12, nyr) 
    82   ICE_vol_export = total(ICE_vol_export,2)/nyr 
     85  ;ICE_vol_export = total(ICE_vol_export,2)/nyr 
     86  ICE_vol_export = ICE_vol_export[*,nyr-1] 
    8387  ; ICE_vol_export = {arr:ICE_vol_export * 1.e-9 * 86400 * 365 , unit : '10^3 Km^3/year'}      ; annual mean 
    8488  ICE_vol_export = {arr:ICE_vol_export * 1.e-9 * 86400 * 30 , unit : '10^3 Km^3/month'}      ; monthly mean 
    8589 
    8690  ; 
    87   title = 'Fram Strait Areal Export'+'!C'+prefix+' '+d1_d2 
     91  ;title = 'Fram Strait Areal Export: LAST YEAR'+'!C'+prefix+' '+d1_d2 
     92  title = 'Fram Strait Areal Export: LAST YEAR'+'!C'+prefix+' '+d2 
    8893  jpt=12 
    8994  time=julday(1,15,1900)+30*lindgen(12) 
    … …  
    99104  xyouts, julday(7,15,1900), 0.41, 'Data from Kwok et al.(2004), 1979-2002 ', ALIGN = 0, CHARTHICK = 2, CHARSIZE=0.8, COLOR=2 
    100105  ; 
    101   title = 'Fram Strait Volume Export'+'!C'+prefix+' '+d1_d2 
     106  ;title = 'Fram Strait Volume Export LAST YEAR'+'!C'+prefix+' '+d1_d2 
     107  title = 'Fram Strait Volume Export LAST YEAR'+'!C'+prefix+' '+d2 
    102108  pltt, ICE_vol_export, 't',  /REMPLI, MIN = 60., MAX = 500. , /NOERASE, XGRIDSTYLE = 1 $ 
    103109       , small = [1, 2, 2], YTITLE = '10^3 Km^3/month', TITLE = title, DATE_FORMAT = '%M', _extra = ex 
    … …  
    118124 
    119125    d1_d2_2 = '('+strtrim(date1_2, 1)+' - '+strtrim(date2_2, 1)+')' 
     126    d2_2 = '('+strtrim(date2_2, 1)+')' 
    120127    tsave = time 
    121128    domdef, indx1, indx2, indy1, indy2, /xindex, /yindex,/memeindices 
    … …  
    140147    ;; AREA 
    141148    ICE_area_export_2 = reform(ICE_area_export_2, 12, nyr) 
    142     ICE_area_export_2 = total(ICE_area_export_2,2)/nyr 
     149    ;ICE_area_export_2 = total(ICE_area_export_2,2)/nyr 
     150    ICE_area_export_2 = ICE_area_export_2[*,nyr-1] 
    143151    ; ICE_area_export_2 = {arr:ICE_area_export_2 * 1.e-12 * 86400 * 365 , unit : '10^6 Km^2/year'}   ; annual mean 
    144152    ICE_area_export_2 = {arr:ICE_area_export_2 * 1.e-12 * 86400 * 30 , unit : '10^6 Km^2/month'}     ; monthly mean 
    145153    ; 
    146154    ICE_vol_export_2 = reform(ICE_vol_export_2, 12, nyr) 
    147     ICE_vol_export_2 = total(ICE_vol_export_2,2)/nyr 
     155    ;ICE_vol_export_2 = total(ICE_vol_export_2,2)/nyr 
     156    ICE_vol_export_2 = ICE_vol_export_2[*,nyr-1] 
    148157    ; ICE_vol_export_2 = {arr:ICE_vol_export_2 * 1.e-12 * 86400 * 365 , unit : '10^3 Km^3/year'}     ; annual mean 
    149158    ICE_vol_export_2 = {arr:ICE_vol_export_2 * 1.e-9 * 86400 * 30 , unit : '10^3 Km^3/month'}       ; monthly mean 
    … …  
    152161    if KEYWORD_SET(postscript) then openps, filename+'_2.ps', portrait = 1 
    153162 
    154     title = 'Fram Strait Areal Export'+'!C'+prefix+' (BLACK) - '+prefix2+' (RED) '+d1_d2_2 
     163    ;title = 'Fram Strait Areal Export'+'!C'+prefix+' (BLACK) - '+prefix2+' (RED) '+d1_d2_2 
     164    title = 'Fram Strait Areal Export LAST YEAR'+'!C'+prefix+' (BLACK) - '+prefix2+' (RED) '+d2_2 
    155165    jpt=12 
    156166    time=julday(1,15,1900)+30*lindgen(12) 
    … …  
    170180    xyouts, julday(7,15,1900), 0.41, 'Data from Kwok et al.(2004), 1979-2002 ', ALIGN = 0, CHARTHICK = 2, CHARSIZE=0.8, COLOR=2 
    171181    ; 
    172     title = 'Fram Strait Volume Export'+'!C'+prefix+' (BLACK) - '+prefix2+' (RED) '+d1_d2_2 
     182    ;title = 'Fram Strait Volume Export'+'!C'+prefix+' (BLACK) - '+prefix2+' (RED) '+d1_d2_2 
     183    title = 'Fram Strait Volume Export LAST YEAR'+'!C'+prefix+' (BLACK) - '+prefix2+' (RED) '+d2_2 
    173184    pltt, ICE_vol_export, 't', /REMPLI,  MIN = 60., MAX = 500., /NOERASE, XGRIDSTYLE = 1 $ 
    174185        , small = [1, 2, 2], YTITLE = '10^3 Km^3/month', TITLE = title, DATE_FORMAT = '%M', _extra = ex 

  • branches/2014/dev_r4650_UKMO10_Tidally_Meaned_Diagnostics/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_ICE_Vol.pro

    r4504 r5260  
    7070; 
    7171  title ='Southern Hemisphere' +'!C'+prefix+' '+d1_d2+' - '+'!C'+'Global Annual Mean Ice Volume (BLACK)' 
    72   pltt, ICE_vol_S, 't', MIN = 0., MAX = 11000., /REMPLI, /NOERASE, XGRIDSTYLE = 1 , DATE_FORMAT = '%M' $ 
     72  pltt, ICE_vol_S, 't', MIN = 0., MAX = 12000., /REMPLI, /NOERASE, XGRIDSTYLE = 1 , DATE_FORMAT = '%M' $ 
    7373        , small = [1, 2, 2], YTITLE = '10^9 Km^3 ', TITLE = title, _extra = ex 
    7474; 
    … …  
    119119    ; 
    120120    title = 'Southern Hemisphere'+'!C'+prefix+' (BLACK) - '+prefix2+' (RED) '+d1_d2_2+'!C'+'Global Annual Mean Ice Volume' 
    121     pltt, ICE_vol_S, 't', MIN = 0., MAX = 11000., /REMPLI, /NOERASE, XGRIDSTYLE = 1, DATE_FORMAT = '%M' $ 
     121    pltt, ICE_vol_S, 't', MIN = 0., MAX = 12000., /REMPLI, /NOERASE, XGRIDSTYLE = 1, DATE_FORMAT = '%M' $ 
    122122          , small = [1, 2, 2], YTITLE = '10^9 Km^3 ', TITLE = title, _extra = ex    
    123123    pltt, ICE_vol_S2, 't',  /REMPLI, /NOERASE $ 

  • branches/2014/dev_r4650_UKMO10_Tidally_Meaned_Diagnostics/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_all.pro

    r4559 r5260  
    6060  IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_SSH, masknp, POSTSCRIPT = postscript, _extra = ex 
    6161; 
    62   cnt = cnt+1   &   blabla = 'Global Mean Q net' 
    63   IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_Q, masknp, POSTSCRIPT = postscript, _extra = ex 
     62;;;;  cnt = cnt+1   &   blabla = 'Global Mean Q net' 
     63;;;;  IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_Q, masknp, POSTSCRIPT = postscript, _extra = ex 
    6464; 
    6565  cnt = cnt+1   &   blabla = 'Global Mean EMP' 
    … …  
    7777  cnt = cnt+1   &   blabla = 'Sea-Ice Volume' 
    7878  IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_ICE_Vol, masknp, POSTSCRIPT = postscript, _extra = ex 
    79 ; 
    80   cnt = cnt+1   &   blabla = 'ICE Velocity ' 
    81   IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_ICE_Vel, masknp, POSTSCRIPT = postscript, _extra = ex 
    82 ; 
     79;SF; and probably useless 
     80;SF  cnt = cnt+1   &   blabla = 'ICE Velocity ' 
     81;SF  IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_ICE_Vel, masknp, POSTSCRIPT = postscript, _extra = ex 
     82;SF; 
    8383  cnt = cnt+1   &   blabla = 'Sea-Ice Volume Export at Fram Strait' 
    8484  IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_ICE_FRAM, masknp, POSTSCRIPT = postscript, _extra = ex 

  • branches/2014/dev_r4650_UKMO10_Tidally_Meaned_Diagnostics/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_vardef.sh

    r4559 r5260  
    1818idl_command=/Applications/itt/idl64/bin/idl 
    1919# 
    20 PS_DIR=$( pwd )/ts_heatbug_ps 
    21 PDF_DIR=$( pwd )/ts_heatbug_pdf 
     20PS_DIR=$( pwd )/ts_p4H25a50-testht_ps 
     21PDF_DIR=$( pwd )/ts_p4H25a50-testht_pdf 
    2222HTML_DIR=$( pwd )/html 
    2323SAXO_DIR=/Users/sflod/SAXO_DIR 
    2424# 
    25 DIR_DATA=/Users/sflod/idl_PLOTS/DATA_STORE/RUN_CLIMATO/lim3_ada      # path of data in NetCDF format 
     25DIR_DATA=/Users/sflod/idl_PLOTS/DATA_STORE/RUN_CLIMATO/lim3_ada     # path of data in NetCDF format 
    2626DIR_CLIMATO=/Users/sflod/idl_PLOTS/CLIMATOLOGIES    # path of climatological data 
    2727DIR_MASK=/Users/sflod/idl_PLOTS/MASK  # path of mask files (ex: subbasins) 
    … …  
    5252READ_ONLY_FIRST_RECORD=0   # if 0 then read all records in files else read only the first reacord in each file 
    5353# 
    54 #DATE1=20010101   ;   DATE2=25001231 
    55 DATE1=20310101   ;   DATE2=20401231 
    56 #DATE1=00010101      ;   DATE2=00101231 
     54DATE1=20010101         ;   DATE2=20501231 
    5755# 
    58 VAR1_T=thetao        ;   V1T_PREF=heatbug       ;   V1T_SUFF=_1Y_grid_T.nc 
    59 VAR1_S=so            ;   V1S_PREF=heatbug      ;   V1S_SUFF=_1Y_grid_T.nc 
    60 VAR1_SSH=zos         ;   V1SSH_PREF=heatbug    ;   V1SSH_SUFF=_1Y_grid_T.nc 
    61 VAR1_Q=tohfls        ;   V1Q_PREF=heatbug      ;   V1Q_SUFF=_1Y_grid_T.nc 
    62 VAR1_EMP=wfo         ;   V1EMP_PREF=heatbug    ;   V1EMP_SUFF=_1Y_grid_T.nc 
    63 VAR1_U=uocetr_eff    ;   V1U_PREF=heatbug      ;   V1U_SUFF=_1Y_grid_U.nc 
    64 VAR1_V=vocetr_eff    ;   V1V_PREF=heatbug      ;   V1V_SUFF=_1Y_grid_V.nc 
    65 VAR1_ICE=siconc      ;   V1ICE_PREF=heatbug    ;   V1ICE_SUFF=_1M_icemod.nc 
    66 VAR1_Ithick=sithic   ;   V1It_PREF=heatbug     ;   V1It_SUFF=_1M_icemod.nc 
    67 VAR1_SNOW=sndept     ;   V1SNOW_PREF=heatbug   ;   V1SNOW_SUFF=_1M_icemod.nc 
    68 VAR1_IvelV=sivelv    ;   V1IvV_PREF=heatbug    ;   V1IvV_SUFF=_1M_icemod.nc 
    69 VAR1_Ivel=sivelo     ;   V1Iv_PREF=heatbug     ;   V1Iv_SUFF=_1M_icemod.nc 
     56VAR1_T=thetao          ;   V1T_PREF=p4H25a50      ;   V1T_SUFF=_1Y_grid_T.nc 
     57VAR1_S=so              ;   V1S_PREF=p4H25a50      ;   V1S_SUFF=_1Y_grid_T.nc 
     58VAR1_SSH=zos         ;   V1SSH_PREF=p4H25a50    ;   V1SSH_SUFF=_1Y_grid_T.nc 
     59##VAR1_Q=qt              ;   V1Q_PREF=p4H25a50      ;   V1Q_SUFF=_1Y_grid_T.nc 
     60VAR1_EMP=wfo         ;   V1EMP_PREF=p4H25a50    ;   V1EMP_SUFF=_1Y_grid_T.nc 
     61VAR1_U=uocetr_eff      ;   V1U_PREF=p4H25a50      ;   V1U_SUFF=_1Y_grid_U.nc 
     62VAR1_V=vocetr_eff      ;   V1V_PREF=p4H25a50      ;   V1V_SUFF=_1Y_grid_V.nc 
     63VAR1_ICE=siconc      ;   V1ICE_PREF=p4H25a50    ;   V1ICE_SUFF=_1M_icemod.nc 
     64VAR1_Ithick=sithic    ;   V1It_PREF=p4H25a50     ;   V1It_SUFF=_1M_icemod.nc 
     65VAR1_SNOW=snthic    ;   V1SNOW_PREF=p4H25a50   ;   V1SNOW_SUFF=_1M_icemod.nc 
     66VAR1_IvelV=sivelv    ;   V1IvV_PREF=p4H25a50    ;   V1IvV_SUFF=_1M_icemod.nc  
     67VAR1_Ivel=sivelo      ;   V1Iv_PREF=p4H25a50     ;   V1Iv_SUFF=_1M_icemod.nc 
    7068# 
    7169#===================== EXP2 ===================== 
    7270# 
    73 #DATE1_2=20010101   ;   DATE2_2=20101231 
    74 DATE1_2=20310101   ;   DATE2_2=20401231 
    75 #DATE1_2=20010101   ;   DATE2_2=25001231 
     71DATE1_2=20010101       ;   DATE2_2=20501231 
    7672# 
    77 VAR2_T=thetao      ;   V2T_PREF=heatbug      ;   V2T_SUFF=_1Y_grid_T.nc 
    78 VAR2_S=so          ;   V2S_PREF=heatbug      ;   V2S_SUFF=_1Y_grid_T.nc 
    79 VAR2_SSH=zos       ;   V2SSH_PREF=heatbug    ;   V2SSH_SUFF=_1Y_grid_T.nc 
    80 VAR2_Q=tohfls      ;   V2Q_PREF=heatbug      ;   V2Q_SUFF=_1Y_grid_T.nc 
    81 VAR2_EMP=wfo       ;   V2EMP_PREF=heatbug    ;   V2EMP_SUFF=_1Y_grid_T.nc 
    82 VAR2_U=uocetr_eff  ;   V2U_PREF=heatbug      ;   V2U_SUFF=_1Y_grid_U.nc 
    83 VAR2_V=vocetr_eff  ;   V2V_PREF=heatbug      ;   V2V_SUFF=_1Y_grid_V.nc 
    84 VAR2_ICE=siconc    ;   V2ICE_PREF=heatbug    ;   V2ICE_SUFF=_1M_icemod.nc 
    85 VAR2_Ithick=sithic ;   V2It_PREF=heatbug     ;   V2It_SUFF=_1M_icemod.nc 
    86 VAR2_SNOW=sndept   ;   V2SNOW_PREF=heatbug   ;   V2SNOW_SUFF=_1M_icemod.nc 
    87 VAR2_IvelV=sivelv  ;   V2IvV_PREF=heatbug    ;   V2IvV_SUFF=_1M_icemod.nc 
    88 VAR2_Ivel=sivelo   ;   V2Iv_PREF=heatbug     ;   V2Iv_SUFF=_1M_icemod.nc 
     73VAR2_T=thetao          ;   V2T_PREF=testht     ;   V2T_SUFF=_1Y_grid_T.nc 
     74VAR2_S=so           ;   V2S_PREF=testht     ;   V2S_SUFF=_1Y_grid_T.nc 
     75VAR2_SSH=zos         ;   V2SSH_PREF=testht   ;   V2SSH_SUFF=_1Y_grid_T.nc 
     76##VAR2_Q=qt            ;   V2Q_PREF=testht     ;   V2Q_SUFF=_1Y_grid_T.nc 
     77VAR2_EMP=wfo            ;   V2EMP_PREF=testht   ;   V2EMP_SUFF=_1Y_grid_T.nc 
     78VAR2_U=uocetr_eff      ;   V2U_PREF=testht     ;   V2U_SUFF=_1Y_grid_U.nc 
     79VAR2_V=vocetr_eff      ;   V2V_PREF=testht     ;   V2V_SUFF=_1Y_grid_V.nc 
     80VAR2_ICE=siconc         ;   V2ICE_PREF=testht   ;   V2ICE_SUFF=_1M_icemod.nc 
     81VAR2_Ithick=sithic       ;   V2It_PREF=testht    ;   V2It_SUFF=_1M_icemod.nc 
     82VAR2_SNOW=snthic       ;   V2SNOW_PREF=testht  ;   V2SNOW_SUFF=_1M_icemod.nc 
     83VAR2_IvelV=sivelv       ;   V2IvV_PREF=testht   ;   V2IvV_SUFF=_1M_icemod.nc  
     84VAR2_Ivel=sivelo        ;   V2IvV_PREF=testht   ;   V2IvV_SUFF=_1M_icemod.nc  
    8985# 
    9086######################### Export Variables ############################### 

  • branches/2014/dev_r4650_UKMO10_Tidally_Meaned_Diagnostics/NEMOGCM/CONFIG/ORCA2_LIM/cpp_ORCA2_LIM.fcm

    r4523 r5260  
    1  bld::tool::fppkeys key_trabbl key_lim2 key_dynspg_flt key_diaeiv key_ldfslp key_traldf_c2d key_traldf_eiv key_dynldf_c3d key_zdftke key_zdfddm key_zdftmx key_iomput key_mpp_mpi 
     1bld::tool::fppkeys key_trabbl key_lim2 key_dynspg_flt key_diaeiv key_ldfslp key_traldf_c2d key_traldf_eiv key_dynldf_c3d key_zdftke key_zdfddm key_zdftmx key_iomput key_mpp_mpi key_diaobs key_asminc 
Note: See TracChangeset for help on using the changeset viewer.