Info
Alle Inhalte des Nutzerportal sind nur auf Englisch verfügbar.

ex_netcdf_1.ncl

EGU 2014 Splinter Meeting SPM1.12 "NCL light hands-on"

Script:

;---------------------------------------------------------------
;-- DKRZ NCL Light Hands-On                   EGU 2014   SPM1.12
;-- 
;-- Example script:       ex_netcdf_1.ncl        (basic script)
;-- 
;-- Settings:             convert temperature in Kelvin to Celsius,
;--                       write variable to new netcdf file
;--
;-- 2014-03-22   Karin Meier-Fleischer ([Email protection active, please enable JavaScript.])
;---------------------------------------------------------------
begin
  outfile = "t_in_Celsius_1.nc"
  
  if (isfilepresent(outfile)) then
    system("rm -rf "+outfile)          ;-- make sure that file does not exist
  end if
  
  fin  = addfile("ECHAM5_OM_A1B_t_20010101.nc","r")  ;-- open data file
  fout = addfile(outfile,"c")          ;-- create new file (netcdf 3)

  filedimdef(fout,"time",-1,True)      ;-- make time and UNLIMITED dimension

  tK           =  fin->t               ;-- get variable
  tC           =  tK                   ;-- copy variable and its related dimensions
                                       ;   and attributes
  tC           =  tK - 273.15          ;-- convert from Kelvin to Celsius
  tC@units     = "Celsius"             ;-- define new units

  fout->tC     =  tC                   ;-- write variable to new netCDF file
  
end

 

Artikelaktionen