;####### !!! WARNING !!! ### !!! WARNING !!! ### !!! WARNING !!! ######### ;## system procedure in use line 139 --> system("rm test_file.nc") ## ;######################################################################### data = asciiread ("TOA5_STUP_10s_5.dat",-1,"string") variable_5 = (str_get_field(data(1::),5,",")) print(variable_5(0)); -> is variable name of culom 5 ;variable_5(0)=variable_5(3:end) this TIMESTAMP=(str_get_field(data(4::),1,",")) RECOR=stringtofloat(str_get_field(data(4::),2,",")) STUP_1_Rain_1h=stringtofloat(str_get_field(data(4::),3,",")) STUP_1_Tlak=stringtofloat(str_get_field(data(4::),4,",")) STUP_2_T=stringtofloat(str_get_field(data(4::),5,",")) STUP_2_RH=stringtofloat(str_get_field(data(4::),6,",")) STUP_4_T=stringtofloat(str_get_field(data(4::),7,",")) STUP_4_RH=stringtofloat(str_get_field(data(4::),8,",")) STUP_4_PAR=stringtofloat(str_get_field(data(4::),9,",")) STUP_4_PARW=stringtofloat(str_get_field(data(4::),10,",")) STUP_4_Vh=stringtofloat(str_get_field(data(4::),11,",")) STUP_10_U_Avg=stringtofloat(str_get_field(data(4::),12,",")) STUP_10_V_Avg=stringtofloat(str_get_field(data(4::),13,",")) STUP_10_W_Avg=stringtofloat(str_get_field(data(4::),14,",")) STUP_10_Vh_Avg=stringtofloat(str_get_field(data(4::),15,",")) STUP_10_V3_Avg=stringtofloat(str_get_field(data(4::),16,",")) STUP_10_Dir_Avg=stringtofloat(str_get_field(data(4::),17,",")) TRS1_T=stringtofloat(str_get_field(data(4::),18,",")) TRS1_RH=stringtofloat(str_get_field(data(4::),19,",")) TRS1_UV=stringtofloat(str_get_field(data(4::),20,",")) TRS1_U_Avg=stringtofloat(str_get_field(data(4::),21,",")) TRS1_V_Avg=stringtofloat(str_get_field(data(4::),22,",")) TRS1_W_Avg=stringtofloat(str_get_field(data(4::),23,",")) TRS1_Vh_Avg=stringtofloat(str_get_field(data(4::),24,",")) TRS1_V3_Avg=stringtofloat(str_get_field(data(4::),25,",")) TRS1_Dir_Avg=stringtofloat(str_get_field(data(4::),26,",")) TRS2_T=stringtofloat(str_get_field(data(4::),27,",")) TRS2_RH=stringtofloat(str_get_field(data(4::),28,",")) TRS2_UV=stringtofloat(str_get_field(data(4::),29,",")) TRS2_U_Avg=stringtofloat(str_get_field(data(4::),30,",")) TRS2_V_Avg=stringtofloat(str_get_field(data(4::),31,",")) TRS2_W_Avg=stringtofloat(str_get_field(data(4::),32,",")) TRS2_Vh_Avg=stringtofloat(str_get_field(data(4::),33,",")) TRS2_V3_Avg=stringtofloat(str_get_field(data(4::),34,",")) TRS2_Dir_Avg=stringtofloat(str_get_field(data(4::),35,",")) dims=dimsizes(TIMESTAMP) dtime_map = (/1,4,1,2,1,2,1,2,1,2,1,2/) dt = str_split_by_length(TIMESTAMP, dtime_map) print(dims) time=new((/dims/),double) time(:) = cd_inv_calendar( tointeger(dt(:,1)), tointeger(dt(:,3)), tointeger(dt(:,5)), tointeger(dt(:,7)), tointeger(dt(:,9)), tointeger(dt(:,11)), "seconds since 2017-06-01 00:00:00", 0 ) time!0 = "time" time@long_name = "time" time@units = "seconds since 2017-06-01 00:00:00" print(time) STUP_1_Rain_1h !0 ="time" STUP_1_Tlak!0 ="time" STUP_2_T!0 ="time" STUP_2_RH!0 ="time" STUP_4_T!0 ="time" STUP_4_RH!0 ="time" STUP_4_PAR!0 ="time" STUP_4_PARW!0 ="time" STUP_4_Vh!0 ="time" STUP_10_U_Avg!0 ="time" STUP_10_V_Avg!0 ="time" STUP_10_W_Avg!0 ="time" STUP_10_Vh_Avg!0 ="time" STUP_10_V3_Avg!0 ="time" STUP_10_Dir_Avg!0 ="time" TRS1_T!0 ="time" TRS1_RH!0 ="time" TRS1_UV!0 ="time" TRS1_U_Avg!0 ="time" TRS1_V_Avg!0 ="time" TRS1_W_Avg!0 ="time" TRS1_Vh_Avg!0 ="time" TRS1_V3_Avg!0 ="time" TRS1_Dir_Avg!0 ="time" TRS2_T!0 ="time" TRS2_RH!0 ="time" TRS2_UV!0 ="time" TRS2_U_Avg!0 ="time" TRS2_V_Avg!0 ="time" TRS2_W_Avg!0 ="time" TRS2_Vh_Avg!0 ="time" TRS2_V3_Avg!0 ="time" TRS2_Dir_Avg!0 ="time" STUP_1_Rain_1h&time = time STUP_1_Tlak&time = time STUP_2_T&time = time STUP_2_RH&time = time STUP_4_T&time = time STUP_4_RH&time = time STUP_4_PAR&time = time STUP_4_PARW&time = time STUP_4_Vh&time = time STUP_10_U_Avg&time = time STUP_10_V_Avg&time = time STUP_10_W_Avg&time = time STUP_10_Vh_Avg&time = time STUP_10_V3_Avg&time = time STUP_10_Dir_Avg&time = time TRS1_T&time = time TRS1_RH&time = time TRS1_UV&time = time TRS1_U_Avg&time = time TRS1_V_Avg&time = time TRS1_W_Avg&time = time TRS1_Vh_Avg&time = time TRS1_V3_Avg&time = time TRS1_Dir_Avg&time = time TRS2_T&time = time TRS2_RH&time = time TRS2_UV&time = time TRS2_U_Avg&time = time TRS2_V_Avg&time = time TRS2_W_Avg&time = time TRS2_Vh_Avg&time = time TRS2_V3_Avg&time = time TRS2_Dir_Avg&time = time system("rm test_file.nc") a = addfile("test_file.nc","c") a->STUP_1_Rain_1h=STUP_1_Rain_1h a->STUP_1_Tlak=STUP_1_Tlak a->STUP_2_T=STUP_2_T a->STUP_2_RH=STUP_2_RH a->STUP_4_T=STUP_4_T a->STUP_4_RH=STUP_4_RH a->STUP_4_PAR=STUP_4_PAR a->STUP_4_PARW=STUP_4_PARW a->STUP_4_Vh=STUP_4_Vh a->STUP_10_U_Avg=STUP_10_U_Avg a->STUP_10_V_Avg=STUP_10_V_Avg a->STUP_10_W_Avg=STUP_10_W_Avg a->STUP_10_Vh_Avg=STUP_10_Vh_Avg a->STUP_10_V3_Avg=STUP_10_V3_Avg a->STUP_10_Dir_Avg=STUP_10_Dir_Avg a->TRS1_T=TRS1_T a->TRS1_RH=TRS1_RH a->TRS1_UV=TRS1_UV a->TRS1_U_Avg=TRS1_U_Avg a->TRS1_V_Avg=TRS1_V_Avg a->TRS1_W_Avg=TRS1_W_Avg a->TRS1_Vh_Avg=TRS1_Vh_Avg a->TRS1_V3_Avg=TRS1_V3_Avg a->TRS1_Dir_Avg=TRS1_Dir_Avg a->TRS2_T=TRS2_T a->TRS2_RH=TRS2_RH a->TRS2_UV=TRS2_UV a->TRS2_U_Avg=TRS2_U_Avg a->TRS2_V_Avg=TRS2_V_Avg a->TRS2_W_Avg=TRS2_W_Avg a->TRS2_Vh_Avg=TRS2_Vh_Avg a->TRS2_V3_Avg=TRS2_V3_Avg a->TRS2_Dir_Avg=TRS2_Dir_Avg exit ;""""""" END OF CODE """""""""""""""""""""""""""" ;################################################ lon=fspan(38,49,1110); lat=fspan(6,20.5,1450); lat!0 = "lat" lat@long_name = "latitude" lat@units = "degrees_north" lon!0 = "lon" lon@long_name = "longitude" lon@units = "degrees_east" data!0 ="lat" data!1 ="lon" data&lat = lat data&lon = lon data1 = data(lon|:,lat|:)