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 13981 for NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/plot_station_asf_ICE.py – NEMO

Ignore:
Timestamp:
2020-12-02T11:50:03+01:00 (3 years ago)
Author:
laurent
Message:

Canceled changes in dynatf.F90, NEW "iom_put-able" "[uv]tau_oce" and "[uv]tau_ice" at T-points!

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/plot_station_asf_ICE.py

    r13893 r13981  
    3939############################################ 
    4040crealm = 'sea-ice' 
    41 L_VNEM = [   'Cd_ice',   'Ch_ice',   'qla_ice' ,   'qsb_ice'   ,  'qt_ice'    ,  'qlw_ice'  ,  'qsr_ice'  , 'taum_ai'   ] 
     41L_VNEM = [   'Cd_ice',   'Ch_ice',   'qla_ice' ,   'qsb_ice'   ,  'qt_ice'    ,  'qlw_ice'  ,  'qsr_ice'  , 'taum_ice'  ] 
    4242L_VARO = [     'Cd'  ,     'Ch'  ,   'Qlat'    ,    'Qsen'     ,   'Qnet'     ,   'Qlw'     ,    'Qsw'    ,  'Tau'      ] 
    4343L_VARL = [ r'$C_{D}$', r'$C_{H}$', r'$Q_{lat}$', r'$Q_{sens}$' , r'$Q_{net}$' , r'$Q_{lw}$' , r'$Q_{sw}$' , r'$|\tau|$' ] 
     
    144144    sys.exit(0) 
    145145 
    146 # >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    147 # >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
     146 
     147#----------------------------------------------------------------- 
    148148 
    149149 
     
    201201        id_toolarge, = nmp.where( xF[:,ja] > L_MAXT[jv] ) #  
    202202        xF[id_toolarge,ja] = L_MAXT[jv] 
    203         id_toosmall, = nmp.where( xF[:,ja] < L_MINT[jv] ) ; print("id_toosmall =", id_toosmall) 
     203        id_toosmall, = nmp.where( xF[:,ja] < L_MINT[jv] ) ; #print("id_toosmall =", id_toosmall) 
    204204        xF[id_toosmall,ja] = L_MINT[jv] 
    205205 
     
    264264 
    265265            for ja in range(nb_exp): 
    266                 fplot = nmp.ma.masked_where( xF[:,ja]==0., xF[:,ja] ) 
     266                fplot = nmp.ma.masked_where( xFa[:,ja]==0., xFa[:,ja] ) 
    267267                plt.plot(vtime, fplot, '-', color=l_color[ja], \ 
    268268                         linewidth=l_width[ja], label=list_exp[ja], alpha=0.6) #, zorder=10+ja) 
Note: See TracChangeset for help on using the changeset viewer.