Ignore:
Timestamp:
2020-07-08T12:47:03+02:00 (3 months ago)
Author:
laurent
Message:

Slight improvements in STATION_ASF test-case

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/trunk/tests/STATION_ASF/EXPREF/plot_station_asf.py

    r12629 r13264  
    2626sys.setdefaultencoding('utf8') 
    2727 
    28 cy1     = '2016' ; # First year 
     28cy1     = '2018' ; # First year 
    2929cy2     = '2018' ; # Last year 
    30  
    31 jt0 = 0 
    32 jt0 = 17519 
    33  
    3430 
    3531dir_figs='.' 
     
    104100# Getting time array from the first file: 
    105101id_in = Dataset(cf_in[0]) 
    106 vt = id_in.variables['time_counter'][jt0:] 
     102vt = id_in.variables['time_counter'][:] 
    107103cunit_t = id_in.variables['time_counter'].units ; print(' "time_counter" is in "'+cunit_t+'"') 
    108104id_in.close() 
     
    138134            if ctest == 'skin':   id_in = Dataset(cf_in[ja]) 
    139135            if ctest == 'noskin': id_in = Dataset(cf_in_ns[ja]) 
    140             xF[:,ja] = id_in.variables[L_VNEM[jv]][jt0:,1,1] # only the center point of the 3x3 spatial domain! 
     136            xF[:,ja] = id_in.variables[L_VNEM[jv]][:,1,1] # only the center point of the 3x3 spatial domain! 
    141137            if ja == 0: cvar_lnm = id_in.variables[L_VNEM[jv]].long_name 
    142138            id_in.close() 
     
    211207    for ja in range(nb_algos-1): 
    212208        id_in = Dataset(cf_in[ja]) 
    213         xF[:,ja]   = id_in.variables[L_VNEM[jv]][jt0:,1,1] # only the center point of the 3x3 spatial domain! 
     209        xF[:,ja]   = id_in.variables[L_VNEM[jv]][:,1,1] # only the center point of the 3x3 spatial domain! 
    214210        if ja == 0: cvar_lnm = id_in.variables[L_VNEM[jv]].long_name 
    215211        id_in.close() 
    216212        # 
    217213        id_in = Dataset(cf_in_ns[ja]) 
    218         xFns[:,ja] = id_in.variables[L_VNEM[jv]][jt0:,1,1] # only the center point of the 3x3 spatial domain! 
     214        xFns[:,ja] = id_in.variables[L_VNEM[jv]][:,1,1] # only the center point of the 3x3 spatial domain! 
    219215        if ja == 0: cvar_lnm = id_in.variables[L_VNEM[jv]].long_name 
    220216        id_in.close() 
Note: See TracChangeset for help on using the changeset viewer.