Changes between Version 5 and Version 6 of Documentation/UserGuide/Perturb

Timestamp:

2017-01-19T10:36:55+01:00 (7 years ago)

Author:

luyssaert

Comment:

--

Documentation/UserGuide/Perturb

@@ -50 +50 @@
 fullfile_name = path+file_name
 ncfile = Dataset(fullfile_name, mode='r+')
-lon = ncfile.variables['rlonu']
-lat = ncfile.variables['rlatu']
-lev = ncfile.variables['nivsigs']
+n_lon = len(ncfile.variables['rlonu'])
+n_lat = len(ncfile.variables['rlatu'])
+n_lev = len(ncfile.variables['nivsigs'])
 teta = ncfile.variables['teta']
-n_lon = len(lon)
-n_lat = len(lat)
-n_lev = len(lev)
-print teta[:,:,:,:]
 perturbation = np.random.normal(mu, sigma,[1,n_lev,n_lat,n_lon])
 ncfile.variables['teta'][:,:,:,:] = teta[:,:,:,:] + perturbation[:,:,:,:]