Monitoring

Load the ferret extension



In [1]:

    
%load_ext ferretmagic

Build a list of files to open

Get from http://vesg.ipsl.upmc.fr/thredds/catalog/work/p86caub/IPSLCM5A-MR/PROD/historicalGHGNOLU/catalog.html



In [4]:

    
listfiles = (
"http://vesg.ipsl.upmc.fr/thredds/dodsC/work/p86caub/IPSLCM6/PROD/historical/CM61-LR-hist-03.1950-REDO-1850/MONITORING/files/ATM_precip_global_ave.nc",
"http://vesg.ipsl.upmc.fr/thredds/dodsC/work/p86caub/IPSLCM6/PROD/historical/CM61-LR-hist-03.1950-REDO-1900/MONITORING/files/ATM_precip_global_ave.nc",
"http://vesg.ipsl.upmc.fr/thredds/dodsC/work/p86caub/IPSLCM6/PROD/historical/CM61-LR-hist-03.1950-REDO-1950/MONITORING/files/ATM_precip_global_ave.nc",
)
print len(listfiles)



In [5]:

    
for f in listfiles:
    print f
    %ferret_run 'use "%s"' % f









    



http://vesg.ipsl.upmc.fr/thredds/dodsC/work/p86caub/IPSLCM6/PROD/historical/CM61-LR-hist-03.1950-REDO-1850/MONITORING/files/ATM_precip_global_ave.nc
http://vesg.ipsl.upmc.fr/thredds/dodsC/work/p86caub/IPSLCM6/PROD/historical/CM61-LR-hist-03.1950-REDO-1900/MONITORING/files/ATM_precip_global_ave.nc
http://vesg.ipsl.upmc.fr/thredds/dodsC/work/p86caub/IPSLCM6/PROD/historical/CM61-LR-hist-03.1950-REDO-1950/MONITORING/files/ATM_precip_global_ave.nc



In [6]:

    
%%ferret
show data









    




     currently SET data sets:
    1> http://vesg.ipsl.upmc.fr/thredds/dodsC/work/p86caub/IPSLCM6/PROD/historical/CM61-LR-hist-03.1950-REDO-1850/MONITORING/files/ATM_precip_global_ave.nc
 name     title                             I         J         K         L
 PRECIP_GLOBAL
                                           ...       ...       ...       1:600
 
    2> http://vesg.ipsl.upmc.fr/thredds/dodsC/work/p86caub/IPSLCM6/PROD/historical/CM61-LR-hist-03.1950-REDO-1900/MONITORING/files/ATM_precip_global_ave.nc
 name     title                             I         J         K         L
 PRECIP_GLOBAL
                                           ...       ...       ...       1:600
 
    3> http://vesg.ipsl.upmc.fr/thredds/dodsC/work/p86caub/IPSLCM6/PROD/historical/CM61-LR-hist-03.1950-REDO-1950/MONITORING/files/ATM_precip_global_ave.nc  (default)
 name     title                             I         J         K         L
 PRECIP_GLOBAL
                                           ...       ...       ...       1:600



In [7]:

    
%%ferret -s 800,400

cancel mode logo
set text/font=arial

ppl color 2, 0, 50, 100, 50
ppl color 3, 100, 50, 0, 50
ppl color 4, 80, 0, 80, 50
ppl color 6, 70, 70, 70

let/title="GHGNOLUMR1" var1 = PRECIP_GLOBAL[l=@sbx:120,d=1]
let/title="GHGNOLUMR2" var2 = PRECIP_GLOBAL[l=@sbx:120,d=2]
let/title="GHGNOLUMR3" var3 = PRECIP_GLOBAL[l=@sbx:120,d=3]

plot/color/thick=3/grat=(dash,color=6)/key=title var1, var2, var3



In [ ]: