@@ -169,7 +169,7 @@ def create_mgcm_atmos_average(short=False):
169169
170170 # Shorten file length if wanted
171171 if short :
172- len_time = 2
172+ len_time = 5
173173 else :
174174 len_time = 133
175175
@@ -287,16 +287,16 @@ def create_mgcm_atmos_average(short=False):
287287
288288 dst_num_micro_var = nc_file .createVariable ('dst_num_micro' , 'f4' , ('time' , 'pfull' , 'lat' , 'lon' ))
289289 dst_num_micro_var .long_name = 'dust_number'
290- dst_num_micro_var [:] = np .random .uniform (- 3.8e-15 , 6.3e+10 , size = (133 , 30 , 48 , 96 ))
290+ dst_num_micro_var [:] = np .random .uniform (- 3.8e-15 , 6.3e+10 , size = (len_time , 30 , 48 , 96 ))
291291
292292 ice_mass_micro_var = nc_file .createVariable ('ice_mass_micro' , 'f4' , ('time' , 'pfull' , 'lat' , 'lon' ))
293293 ice_mass_micro_var .long_name = 'ice_mass'
294- ice_mass_micro_var [:] = np .random .uniform (- 5.8e-34 , 3.1e-03 , size = (133 , 30 , 48 , 96 ))
294+ ice_mass_micro_var [:] = np .random .uniform (- 5.8e-34 , 3.1e-03 , size = (len_time , 30 , 48 , 96 ))
295295
296296 omega_var = nc_file .createVariable ('omega' , 'f4' , ('time' , 'pfull' , 'lat' , 'lon' ))
297297 omega_var .long_name = 'vertical wind'
298298 omega_var .units = 'Pa/s'
299- omega_var [:] = np .random .uniform (- 0.045597 , 0.0806756 , size = (133 , 30 , 48 , 96 ))
299+ omega_var [:] = np .random .uniform (- 0.045597 , 0.0806756 , size = (len_time , 30 , 48 , 96 ))
300300
301301 ps_var = nc_file .createVariable ('ps' , 'f4' , ('time' , 'lat' , 'lon' ))
302302 ps_var .long_name = 'surface pressure'
@@ -342,7 +342,7 @@ def create_mgcm_atmos_daily(short=False):
342342
343343 # Shorten file length if wanted
344344 if short :
345- len_time = 2
345+ len_time = 40
346346 else :
347347 len_time = 2672
348348
@@ -357,7 +357,10 @@ def create_mgcm_atmos_daily(short=False):
357357 time_var = nc_file .createVariable ('time' , 'f4' , ('time' ,))
358358 time_var .long_name = 'time'
359359 time_var .units = 'days'
360- time_var [:] = np .linspace (1336.2 , 2004.0 , len_time )
360+ if short :
361+ time_var [:] = np .linspace (1336.2 , 1336.2 + float (len_time )/ 4. , len_time )
362+ else :
363+ time_var [:] = np .linspace (1336.2 , 2004.0 , len_time )
361364
362365 lat_var = nc_file .createVariable ('lat' , 'f4' , ('lat' ,))
363366 lat_var .long_name = 'latitude'
@@ -405,8 +408,11 @@ def create_mgcm_atmos_daily(short=False):
405408 areo_var = nc_file .createVariable ('areo' , 'f4' , ('time' , 'scalar_axis' ))
406409 areo_var .long_name = 'areo'
407410 areo_var .units = 'degrees'
408- areo_vals = np .linspace (720.3 , 1079.8 , len_time )
409- areo_data = np .zeros ((len_time , 1 )) # Create a 2D array with shape (133, 1)
411+ if short :
412+ areo_vals = np .linspace (720.3 , 720.3 + 0.538 * float (len_time )/ 4. , len_time )
413+ else :
414+ areo_vals = np .linspace (720.3 , 1079.8 , len_time )
415+ areo_data = np .zeros ((len_time , 1 )) # Create a 2D array with shape (len_time, 1)
410416 for i in range (len_time ):
411417 areo_data [i , 0 ] = areo_vals [i ]
412418 areo_var [:] = areo_data
@@ -430,7 +436,7 @@ def create_mgcm_atmos_average_pstd(short=False):
430436
431437 # Shorten file length if wanted
432438 if short :
433- len_time = 2
439+ len_time = 5
434440 else :
435441 len_time = 133
436442
@@ -508,22 +514,22 @@ def create_mgcm_atmos_average_pstd(short=False):
508514 theta_var = nc_file .createVariable ('theta' , 'f4' , ('time' , 'pstd' , 'lat' , 'lon' ))
509515 theta_var .long_name = 'Potential temperature'
510516 theta_var .units = 'K'
511- theta_var [:] = np .random .uniform (104.113 , 3895.69 , size = (133 , 44 , 48 , 96 ))
517+ theta_var [:] = np .random .uniform (104.113 , 3895.69 , size = (len_time , 44 , 48 , 96 ))
512518
513519 rho_var = nc_file .createVariable ('rho' , 'f4' , ('time' , 'pstd' , 'lat' , 'lon' ))
514520 rho_var .long_name = 'Density'
515521 rho_var .units = 'kg/m^3'
516- rho_var [:] = np .random .uniform (7.05091e-07 , 0.0668856 , size = (133 , 44 , 48 , 96 ))
522+ rho_var [:] = np .random .uniform (7.05091e-07 , 0.0668856 , size = (len_time , 44 , 48 , 96 ))
517523
518524 omega_var = nc_file .createVariable ('omega' , 'f4' , ('time' , 'pstd' , 'lat' , 'lon' ))
519525 omega_var .long_name = 'vertical wind'
520526 omega_var .units = 'Pa/s'
521- omega_var [:] = np .random .uniform (- 0.045597 , 0.0806756 , size = (133 , 44 , 48 , 96 ))
527+ omega_var [:] = np .random .uniform (- 0.045597 , 0.0806756 , size = (len_time , 44 , 48 , 96 ))
522528
523529 w_var = nc_file .createVariable ('w' , 'f4' , ('time' , 'pstd' , 'lat' , 'lon' ))
524530 w_var .long_name = 'w'
525531 w_var .units = 'm/s'
526- w_var [:] = np .random .uniform (- 2.02603 , 1.58804 , size = (133 , 44 , 48 , 96 ))
532+ w_var [:] = np .random .uniform (- 2.02603 , 1.58804 , size = (len_time , 44 , 48 , 96 ))
527533
528534 ps_var = nc_file .createVariable ('ps' , 'f4' , ('time' , 'lat' , 'lon' ))
529535 ps_var .long_name = 'surface pressure'
@@ -569,7 +575,7 @@ def create_mgcm_atmos_diurn_pstd(short=False):
569575
570576 # Shorten file length if wanted
571577 if short :
572- len_time = 2
578+ len_time = 5
573579 else :
574580 len_time = 133
575581
@@ -662,7 +668,7 @@ def create_mgcm_atmos_diurn(short=False):
662668
663669 # Shorten file length if wanted
664670 if short :
665- len_time = 2
671+ len_time = 5
666672 else :
667673 len_time = 133
668674
@@ -773,7 +779,7 @@ def create_mgcm_atmos_average_pstd_c48(short=False):
773779
774780 # Shorten file length if wanted
775781 if short :
776- len_time = 2
782+ len_time = 5
777783 else :
778784 len_time = 133
779785
0 commit comments