CLM4.0 自己造强迫场单点run出错

jdy

以下是我造强迫场的文件~

begin
  filename = "/home/ash/clm4.0/clm4run-out/2013maduoqiangpochang/201308maduoqiangpo.csv"

;---Read in file as array of strings so we can parse each line.
  lines  = asciiread(filename,-1,"string")
  nlines = dimsizes(lines)   ; First line is a header

; The first line is a header, so skip this.
tbot=new(52512,float)
rh  =new(52512,float)
wind=new(52512,float)
pres=new(52512,float)
prate=new(52512,float)
fsds=new(52512,float)
flds=new(52512,float)

tbot1=new((/1488,1,1/),float)
rh1  =new((/1488,1,1/),float)
wind1=new((/1488,1,1/),float)
pres1=new((/1488,1,1/),float)
prate1=new((/1488,1,1/),float)
fsds1=new((/1488,1,1/),float)
flds1=new((/1488,1,1/),float)

tbot2=new((/1344,1,1/),float)
rh2  =new((/1344,1,1/),float)
wind2=new((/1344,1,1/),float)
pres2=new((/1344,1,1/),float)
prate2=new((/1344,1,1/),float)
fsds2=new((/1344,1,1/),float)
flds2=new((/1344,1,1/),float)

tbot3=new((/1488,1,1/),float)
rh3  =new((/1488,1,1/),float)
wind3=new((/1488,1,1/),float)
pres3=new((/1488,1,1/),float)
prate3=new((/1488,1,1/),float)
fsds3=new((/1488,1,1/),float)
flds3=new((/1488,1,1/),float)

tbot4=new((/1440,1,1/),float)
rh4  =new((/1440,1,1/),float)
wind4=new((/1440,1,1/),float)
pres4=new((/1440,1,1/),float)
prate4=new((/1440,1,1/),float)
fsds4=new((/1440,1,1/),float)
flds4=new((/1440,1,1/),float)

tbot5=new((/1488,1,1/),float)
rh5  =new((/1488,1,1/),float)
wind5=new((/1488,1,1/),float)
pres5=new((/1488,1,1/),float)
prate5=new((/1488,1,1/),float)
fsds5=new((/1488,1,1/),float)
flds5=new((/1488,1,1/),float)

tbot6=new((/1440,1,1/),float)
rh6  =new((/1440,1,1/),float)
wind6=new((/1440,1,1/),float)
pres6=new((/1440,1,1/),float)
prate6=new((/1440,1,1/),float)
fsds6=new((/1440,1,1/),float)
flds6=new((/1440,1,1/),float)

tbot7=new((/1488,1,1/),float)
rh7  =new((/1488,1,1/),float)
wind7=new((/1488,1,1/),float)
pres7=new((/1488,1,1/),float)
prate7=new((/1488,1,1/),float)
fsds7=new((/1488,1,1/),float)
flds7=new((/1488,1,1/),float)

tbot8=new((/1488,1,1/),float)
rh8  =new((/1488,1,1/),float)
wind8=new((/1488,1,1/),float)
pres8=new((/1488,1,1/),float)
prate8=new((/1488,1,1/),float)
fsds8=new((/1488,1,1/),float)
flds8=new((/1488,1,1/),float)

tbot9=new((/1440,1,1/),float)
rh9  =new((/1440,1,1/),float)
wind9=new((/1440,1,1/),float)
pres9=new((/1440,1,1/),float)
prate9=new((/1440,1,1/),float)
fsds9=new((/1440,1,1/),float)
flds9=new((/1440,1,1/),float)

tbot10=new((/1488,1,1/),float)
rh10  =new((/1488,1,1/),float)
wind10=new((/1488,1,1/),float)
pres10=new((/1488,1,1/),float)
prate10=new((/1488,1,1/),float)
fsds10=new((/1488,1,1/),float)
flds10=new((/1488,1,1/),float)

tbot11=new((/1440,1,1/),float)
rh11  =new((/1440,1,1/),float)
wind11=new((/1440,1,1/),float)
pres11=new((/1440,1,1/),float)
prate11=new((/1440,1,1/),float)
fsds11=new((/1440,1,1/),float)
flds11=new((/1440,1,1/),float)

tbot12=new((/1488,1,1/),float)
rh12  =new((/1488,1,1/),float)
wind12=new((/1488,1,1/),float)
pres12=new((/1488,1,1/),float)
prate12=new((/1488,1,1/),float)
fsds12=new((/1488,1,1/),float)
flds12=new((/1488,1,1/),float)

z1=new((/1488,1,1/),float)
z2=new((/1344,1,1/),float)
z3=new((/1440,1,1/),float)

prate_t200b=new((/17520,1,1/),float)


delim     = ","

  tbot = stringtofloat(str_get_field(lines(0:),1,delim))
  rh = stringtofloat(str_get_field(lines(0:),2,delim))
  wind = stringtofloat(str_get_field(lines(0:),3,delim))
  pres = stringtofloat(str_get_field(lines(0:),4,delim))
  prate = stringtofloat(str_get_field(lines(0:),5,delim))
  fsds = stringtofloat(str_get_field(lines(0:),6,delim))
  flds = stringtofloat(str_get_field(lines(0:),7,delim))
print(rh(0:7))



ncdf1  = addfile("/home/ash/clm4.0/clm4run-out/2013maduoqiangpochang/2013-01.nc" ,"c")  ; open output netCDF file
ncdf2  = addfile("/home/ash/clm4.0/clm4run-out/2013maduoqiangpochang/2013-02.nc" ,"c")  ; open output netCDF file
ncdf3  = addfile("/home/ash/clm4.0/clm4run-out/2013maduoqiangpochang/2013-03.nc" ,"c")  ; open output netCDF file
ncdf4  = addfile("/home/ash/clm4.0/clm4run-out/2013maduoqiangpochang/2013-04.nc" ,"c")  ; open output netCDF file
ncdf5  = addfile("/home/ash/clm4.0/clm4run-out/2013maduoqiangpochang/2013-05.nc" ,"c")  ; open output netCDF file
ncdf6  = addfile("/home/ash/clm4.0/clm4run-out/2013maduoqiangpochang/2013-06.nc" ,"c")  ; open output netCDF file
ncdf7  = addfile("/home/ash/clm4.0/clm4run-out/2013maduoqiangpochang/2013-07.nc" ,"c")  ; open output netCDF file
ncdf8  = addfile("/home/ash/clm4.0/clm4run-out/2013maduoqiangpochang/2013-08.nc" ,"c")  ; open output netCDF file
ncdf9  = addfile("/home/ash/clm4.0/clm4run-out/2013maduoqiangpochang/2013-09.nc" ,"c")  ; open output netCDF file
ncdf10  = addfile("/home/ash/clm4.0/clm4run-out/2013maduoqiangpochang/2013-10.nc" ,"c")  ; open output netCDF file
ncdf11  = addfile("/home/ash/clm4.0/clm4run-out/2013maduoqiangpochang/2013-11.nc" ,"c")  ; open output netCDF file
ncdf12  = addfile("/home/ash/clm4.0/clm4run-out/2013maduoqiangpochang/2013-12.nc" ,"c")  ; open output netCDF file

        fAtt               = True            ; assign file attributes
        fAtt@title         = "Maduo atmospheric data at 0.5-hour resolution"  
        fAtt@Conventions   = "None"   
        fAtt@creation_date = systemfunc ("date")        

fileattdef( ncdf1, fAtt )            ; copy file attributes
fileattdef( ncdf2, fAtt )
fileattdef( ncdf3, fAtt )
fileattdef( ncdf4, fAtt )
fileattdef( ncdf5, fAtt )
fileattdef( ncdf6, fAtt )
fileattdef( ncdf7, fAtt )
fileattdef( ncdf8, fAtt )
fileattdef( ncdf9, fAtt )
fileattdef( ncdf10, fAtt )
fileattdef( ncdf11, fAtt )
fileattdef( ncdf12, fAtt )

; make time and UNLIMITED dimension      ; recommended  for most applications
filedimdef(ncdf1,"time",-1,True)
filedimdef(ncdf2,"time",-1,True)
filedimdef(ncdf3,"time",-1,True)
filedimdef(ncdf4,"time",-1,True)
filedimdef(ncdf5,"time",-1,True)
filedimdef(ncdf6,"time",-1,True)
filedimdef(ncdf7,"time",-1,True)
filedimdef(ncdf8,"time",-1,True)
filedimdef(ncdf9,"time",-1,True)
filedimdef(ncdf10,"time",-1,True)
filedimdef(ncdf11,"time",-1,True)
filedimdef(ncdf12,"time",-1,True)

;print(ts(13248,:))
do i=0,1487,1
tbot1(i,:,:)=tbot(i)
rh1(i,:,:)  =rh(i)
wind1(i,:,:)=wind(i)
pres1(i,:,:)=pres(i)
prate1(i,:,:)=prate(i)
fsds1(i,:,:)=fsds(i)
flds1(i,:,:)=flds(i)
end do
tbot131=tbot1(0:1487,:,:)
rh131=rh1(0:1487,:,:)
wind131=wind1(0:1487,:,:)
pres131=pres1(0:1487,:,:)*100
prate131=prate1(0:1487,:,:)/1800
;prate131=prate_t200b(0:1487,:,:)
fsds131=fsds1(0:1487,:,:)
flds131=flds1(0:1487,:,:)
print(tbot1(0:7,:,:))

do i=0,1343,1
tbot2(i,:,:)=tbot(i)
rh2(i,:,:)  =rh(i)
wind2(i,:,:)=wind(i)
pres2(i,:,:)=pres(i)
prate2(i,:,:)=prate(i)
fsds2(i,:,:)=fsds(i)
flds2(i,:,:)=flds(i)
end do
tbot21=tbot2(0:1343,:,:)
rh21=rh2(0:1343,:,:)
wind21=wind2(0:1343,:,:)
pres21=pres2(0:1343,:,:)*100
prate21=prate2(0:1343,:,:)/1800
;prate21=prate_t200b(1488:2831,:,:)

fsds21=fsds2(0:1343,:,:)
flds21=flds2(0:1343,:,:)

do i=0,1487,1
tbot3(i,:,:)=tbot(i)
rh3(i,:,:)  =rh(i)
wind3(i,:,:)=wind(i)
pres3(i,:,:)=pres(i)
prate3(i,:,:)=prate(i)
fsds3(i,:,:)=fsds(i)
flds3(i,:,:)=flds(i)
end do
tbot31=tbot3(0:1487,:,:)
rh31=rh3(0:1487,:,:)
wind31=wind3(0:1487,:,:)
pres31=pres3(0:1487,:,:)*100
prate31=prate3(0:1487,:,:)/1800
;prate31=prate_t200b(2832:4319,:,:)
fsds31=fsds3(0:1487,:,:)
flds31=flds3(0:1487,:,:)

do i=0,1439,1
tbot4(i,:,:)=tbot(i)
rh4(i,:,:)  =rh(i)
wind4(i,:,:)=wind(i)
pres4(i,:,:)=pres(i)
prate4(i,:,:)=prate(i)
fsds4(i,:,:)=fsds(i)
flds4(i,:,:)=flds(i)
end do
tbot41=tbot4(0:1439,:,:)
rh41=rh4(0:1439,:,:)
wind41=wind4(0:1439,:,:)
pres41=pres4(0:1439,:,:)*100
prate41=prate4(0:1439,:,:)/1800

fsds41=fsds4(0:1439,:,:)
flds41=flds4(0:1439,:,:)

do i=0,1487,1
tbot5(i,:,:)=tbot(i)
rh5(i,:,:)  =rh(i)
wind5(i,:,:)=wind(i)
pres5(i,:,:)=pres(i)
prate5(i,:,:)=prate(i)+prate(i)+prate(i)
fsds5(i,:,:)=fsds(i)
flds5(i,:,:)=flds(i)
end do
tbot51=tbot5(0:1487,:,:)
rh51=rh5(0:1487,:,:)
wind51=wind5(0:1487,:,:)
pres51=pres5(0:1487,:,:)*100
prate51=prate5(0:1487,:,:)/1800
fsds51=fsds5(0:1487,:,:)
flds51=flds5(0:1487,:,:)


do i=0,1439,1
tbot6(i,:,:)=tbot(i)
rh6(i,:,:)  =rh(i)
wind6(i,:,:)=wind(i)
pres6(i,:,:)=pres(i)
prate6(i,:,:)=prate(i)
fsds6(i,:,:)=fsds(i)
flds6(i,:,:)=flds(i)
end do
tbot61=tbot6(0:1439,:,:)
rh61=rh6(0:1439,:,:)
wind61=wind6(0:1439,:,:)
pres61=pres6(0:1439,:,:)*100
prate61=prate6(0:1439,:,:)/1800
fsds61=fsds6(0:1439,:,:)
flds61=flds6(0:1439,:,:)

do i=0,1487,1
tbot7(i,:,:)=tbot(i)
rh7(i,:,:)  =rh(i)
wind7(i,:,:)=wind(i)
pres7(i,:,:)=pres(i)
prate7(i,:,:)=prate(i)
fsds7(i,:,:)=fsds(i)
flds7(i,:,:)=flds(i)
end do
tbot71=tbot7(0:1487,:,:)
rh71=rh7(0:1487,:,:)
wind71=wind7(0:1487,:,:)
pres71=pres7(0:1487,:,:)*100
prate71=prate7(0:1487,:,:)/1800
fsds71=fsds7(0:1487,:,:)
flds71=flds7(0:1487,:,:)

do i=0,1487,1
tbot8(i,:,:)=tbot(i)
rh8(i,:,:)  =rh(i)
wind8(i,:,:)=wind(i)
pres8(i,:,:)=pres(i)
prate8(i,:,:)=prate(i)
fsds8(i,:,:)=fsds(i)
flds8(i,:,:)=flds(i)
end do
tbot81=tbot8(0:1487,:,:)
rh81=rh8(0:1487,:,:)
wind81=wind8(0:1487,:,:)
pres81=pres8(0:1487,:,:)*100.
prate81=prate8(0:1487,:,:)/1800
fsds81=fsds8(0:1487,:,:)
flds81=flds8(0:1487,:,:)

do i=0,1439,1
tbot9(i,:,:)=tbot(i)
rh9(i,:,:)  =rh(i)
wind9(i,:,:)=wind(i)
pres9(i,:,:)=pres(i)
prate9(i,:,:)=prate(i)
fsds9(i,:,:)=fsds(i)
flds9(i,:,:)=flds(i)
end do
tbot91=tbot9(0:1439,:,:)
rh91=rh9(0:1439,:,:)
wind91=wind9(0:1439,:,:)
pres91=pres9(0:1439,:,:)*100
prate91=prate9(0:1439,:,:)/1800
fsds91=fsds9(0:1439,:,:)
flds91=flds9(0:1439,:,:)

do i=0,1487,1
tbot10(i,:,:)=tbot(i)
rh10(i,:,:)  =rh(i)
wind10(i,:,:)=wind(i)
pres10(i,:,:)=pres(i)
prate10(i,:,:)=prate(i)
fsds10(i,:,:)=fsds(i)
flds10(i,:,:)=flds(i)
end do
tbot101=tbot10(0:1487,:,:)
rh101=rh10(0:1487,:,:)
wind101=wind10(0:1487,:,:)
pres101=pres10(0:1487,:,:)*100
prate101=prate10(0:1487,:,:)/1800
fsds101=fsds10(0:1487,:,:)
flds101=flds10(0:1487,:,:)

do i=0,1439,1
tbot11(i,:,:)=tbot(i)
rh11(i,:,:)  =rh(i)
wind11(i,:,:)=wind(i)
pres11(i,:,:)=pres(i)
prate11(i,:,:)=prate(i)
fsds11(i,:,:)=fsds(i)
flds11(i,:,:)=flds(i)
end do
tbot111=tbot11(0:1439,:,:)
rh111=rh11(0:1439,:,:)
wind111=wind11(0:1439,:,:)
pres111=pres11(0:1439,:,:)*100
prate111=prate11(0:1439,:,:)/1800
;prate111=prate_t200b(14592:16031,:,:)
fsds111=fsds11(0:1439,:,:)
flds111=flds11(0:1439,:,:)

do i=0,1487,1
tbot12(i,:,:)=tbot(i)
rh12(i,:,:)  =rh(i)
wind12(i,:,:)=wind(i)
pres12(i,:,:)=pres(i)
prate12(i,:,:)=prate(i)
fsds12(i,:,:)=fsds(i)
flds12(i,:,:)=flds(i)
end do
tbot121=tbot12(0:1487,:,:)
rh121=rh12(0:1487,:,:)
wind121=wind12(0:1487,:,:)
pres121=pres12(0:1487,:,:)*100
prate121=prate12(0:1487,:,:)/1800


fsds121=fsds12(0:1487,:,:)
flds121=flds12(0:1487,:,:)


lon1=new((/1,1/),float)
lon1!0="lat"
lon1!1="lon"
lon1=96.4
;lon1=102.5
;lon1=102
lon1@long_name="longitude"
lon1@units="degrees E"
lon1@mode="time-invariant"

lat1=new((/1,1/),float)
lat1!0="lat"
lat1!1="lon"
lat1=35
;lat1=34
;lat1=33.5
lat1@long_name="latitude"
lat1@units="degrees N"
lat1@mode="time-invariant"

a=34.95
;a=33.5
;a=32.5
a!0="scalar"
a@long_name="southern edge in atmospheric data"
a@units="degrees N"
a@mode="time-invariant"

b=35.05
;b=34.5
b!0="scalar"
b@long_name="northern edge in atmospheric data"
b@units="degrees N"
b@mode="time-invariant"

c=96.45
;c=103
c!0="scalar"
c@long_name="eastern edge in atmospheric data"
c@units="degrees E"
c@mode="time-invariant"

d=96.35
;d=102
;d=101
d!0="scalar"
d@long_name="western edge in atmospheric data"
d@units="degrees E"
d@mode="time-invariant"

time1=fspan(0,30.97917,1488)
time1!0="time"
time1@long_name="observation time"
time1@units="days since 2013-01-01 00:00:00"
time1@calender="noleap"
;print(time1)

time2=fspan(0,27.97949,1344)
time2!0="time"
time2@long_name="observation time"
time2@units="days since 2013-02-01 00:00:00"
time2@calender="noleap"

time3=fspan(0,30.97917,1488)
time3!0="time"
time3@long_name="observation time"
time3@units="days since 2013-03-01 00:00:00"
time3@calender="noleap"

time4=fspan(0,29.97949,1440)
time4!0="time"
time4@long_name="observation time"
time4@units="days since 2013-04-01 00:00:00"
time4@calender="noleap"

time5=fspan(0,30.97917,1488)
time5!0="time"
time5@long_name="observation time"
time5@units="days since 2013-05-01 00:00:00"
time5@calender="noleap"

time6=fspan(0,29.97949,1440)
time6!0="time"
time6@long_name="observation time"
time6@units="days since 2013-06-01 00:00:00"
time6@calender="noleap"

time7=fspan(0,30.97917,1488)
time7!0="time"
time7@long_name="observation time"
time7@units="days since 2013-07-01 00:00:00"
time7@calender="noleap"

time8=fspan(0,30.97917,1488)
time8!0="time"
time8@long_name="observation time"
time8@units="days since 2013-08-01 00:00:00"
time8@calender="noleap"

time9=fspan(0,29.97949,1440)
time9!0="time"
time9@long_name="observation time"
time9@units="days since 2013-09-01 00:00:00"
time9@calender="noleap"

time10=fspan(0,30.97917,1488)
time10!0="time"
time10@long_name="observation time"
time10@units="days since 2013-10-01 00:00:00"
time10@calender="noleap"

time11=fspan(0,29.97949,1440)
time11!0="time"
time11@long_name="observation time"
time11@units="days since 2013-11-01 00:00:00"
time11@calender="noleap"

time12=fspan(0,29.97949,1440)
time12!0="time"
time12@long_name="observation time"
time12@units="days since 2013-12-01 00:00:00"
time12@calender="noleap"

do i=0,1487,1
z1(i,:,:)=2.0
end do

do i=0,1343,1
z2(i,:,:)=2.0
end do

do i=0,1439,1
z3(i,:,:)=2.0
end do

tbot131@long_name="temperature at the lowest atm level (TBOT)"
tbot131@units="K"
tbot131@mode="time-dependent"
tbot131!0="time"
tbot131&time=time1
tbot131!1="lat"
tbot131!2="lon"
ncdf1->TBOT=tbot131
rh131@long_name="specific humidity at the lowest atm level (RH)"
rh131@units="kg/kg"
rh131@mode="time-dependent"
rh131!0="time"
rh131&time=time1
rh131!1="lat"
rh131!2="lon"
ncdf1->RH=rh131
wind131@long_name="wind at the lowest atm level (WIND)"
wind131@units="m/s"
wind131@mode="time-dependent"
wind131!0="time"
wind131&time=time1
wind131!1="lat"
wind131!2="lon"
ncdf1->WIND=wind131
pres131@long_name="surface pressure at the lowest atm level (PSRF)"
pres131@units="Pa"
pres131@mode="time-dependent"
pres131!0="time"
pres131&time=time1
pres131!1="lat"
pres131!2="lon"
ncdf1->PSRF=pres131
prate131@long_name="preciptation(PRECT)"
prate131@units="mm/s"
prate131@mode="time-dependent"
prate131!0="time"
prate131&time=time1
prate131!1="lat"
prate131!2="lon"
ncdf1->PRECTmms=prate131
fsds131@long_name="incident solar(FSDS)"
fsds131@units="W/m2"
fsds131@mode="time-dependent"
fsds131!0="time"
fsds131&time=time1
fsds131!1="lat"
fsds131!2="lon"
ncdf1->FSDS=fsds131
flds131@long_name="incident longwave radiation(FLDS)"
flds131@units="W/m2"
flds131@mode="time-dependent"
flds131!0="time"
flds131&time=time1
flds131!1="lat"
flds131!2="lon"
ncdf1->FLDS=flds131
z1@long_name="observational height"
z1@units="m"
z1@mode="time-invariant"
z1!0="time"
z1&time=time1
z1!1="lat"
z1!2="lon"
ncdf1->ZBOT = z1
ncdf1->LONGXY = lon1
ncdf1->LATIXY = lat1
ncdf1->EDGES = a
ncdf1->EDGEN = b
ncdf1->EDGEE = c
ncdf1->EDGEW = d


tbot21@long_name="temperature at the lowest atm level (TBOT)"
tbot21@units="K"
tbot21@mode="time-dependent"
tbot21!0="time"
tbot21&time=time2
tbot21!1="lat"
tbot21!2="lon"
ncdf2->TBOT=tbot21
rh21@long_name="specific humidity at the lowest atm level (RH)"
rh21@units="kg/kg"
rh21@mode="time-dependent"
rh21!0="time"
rh21&time=time2
rh21!1="lat"
rh21!2="lon"
ncdf2->RH=rh21
wind21@long_name="wind at the lowest atm level (WIND)"
wind21@units="m/s"
wind21@mode="time-dependent"
wind21!0="time"
wind21&time=time2
wind21!1="lat"
wind21!2="lon"
ncdf2->WIND=wind21
pres21@long_name="surface pressure at the lowest atm level (PSRF)"
pres21@units="Pa"
pres21@mode="time-dependent"
pres21!0="time"
pres21&time=time2
pres21!1="lat"
pres21!2="lon"
ncdf2->PSRF=pres21
prate21@long_name="preciptation(PRECT)"
prate21@units="mm/s"
prate21@mode="time-dependent"
prate21!0="time"
prate21&time=time2
prate21!1="lat"
prate21!2="lon"
ncdf2->PRECTmms=prate21
fsds21@long_name="incident solar(FSDS)"
fsds21@units="W/m2"
fsds21@mode="time-dependent"
fsds21!0="time"
fsds21&time=time2
fsds21!1="lat"
fsds21!2="lon"
ncdf2->FSDS=fsds21
flds21@long_name="incident longwave radiation(FLDS)"
flds21@units="W/m2"
flds21@mode="time-dependent"
flds21!0="time"
flds21&time=time2
flds21!1="lat"
flds21!2="lon"
ncdf2->FLDS=flds21
z2@long_name="observational height"
z2@units="m"
z2@mode="time-invariant"
z2!0="time"
z2&time=time2
z2!1="lat"
z2!2="lon"
ncdf2->ZBOT = z2
ncdf2->LONGXY = lon1
ncdf2->LATIXY = lat1
ncdf2->EDGES = a
ncdf2->EDGEN = b
ncdf2->EDGEE = c
ncdf2->EDGEW = d

tbot31@long_name="temperature at the lowest atm level (TBOT)"
tbot31@units="K"
tbot31@mode="time-dependent"
tbot31!0="time"
tbot31&time=time3
tbot31!1="lat"
tbot31!2="lon"
ncdf3->TBOT=tbot31
rh31@long_name="specific humidity at the lowest atm level (RH)"
rh31@units="kg/kg"
rh31@mode="time-dependent"
rh31!0="time"
rh31&time=time3
rh31!1="lat"
rh31!2="lon"
ncdf3->RH=rh31
wind31@long_name="wind at the lowest atm level (WIND)"
wind31@units="m/s"
wind31@mode="time-dependent"
wind31!0="time"
wind31&time=time3
wind31!1="lat"
wind31!2="lon"
ncdf3->WIND=wind31
pres31@long_name="surface pressure at the lowest atm level (PSRF)"
pres31@units="Pa"
pres31@mode="time-dependent"
pres31!0="time"
pres31&time=time3
pres31!1="lat"
pres31!2="lon"
ncdf3->PSRF=pres31
prate31@long_name="preciptation(PRECT)"
prate31@units="mm/s"
prate31@mode="time-dependent"
prate31!0="time"
prate31&time=time3
prate31!1="lat"
prate31!2="lon"
ncdf3->PRECTmms=prate31
fsds31@long_name="incident solar(FSDS)"
fsds31@units="W/m2"
fsds31@mode="time-dependent"
fsds31!0="time"
fsds31&time=time3
fsds31!1="lat"
fsds31!2="lon"
ncdf3->FSDS=fsds31
flds31@long_name="incident longwave radiation(FLDS)"
flds31@units="W/m2"
flds31@mode="time-dependent"
flds31!0="time"
flds31&time=time3
flds31!1="lat"
flds31!2="lon"
ncdf3->FLDS=flds31
z1@long_name="observational height"
z1@units="m"
z1@mode="time-invariant"
z1!0="time"
z1&time=time3
z1!1="lat"
z1!2="lon"
ncdf3->ZBOT = z1
ncdf3->LONGXY = lon1
ncdf3->LATIXY = lat1
ncdf3->EDGES = a
ncdf3->EDGEN = b
ncdf3->EDGEE = c
ncdf3->EDGEW = d

tbot41@long_name="temperature at the lowest atm level (TBOT)"
tbot41@units="K"
tbot41@mode="time-dependent"
tbot41!0="time"
tbot41&time=time4
tbot41!1="lat"
tbot41!2="lon"
ncdf4->TBOT=tbot41
rh41@long_name="specific humidity at the lowest atm level (RH)"
rh41@units="kg/kg"
rh41@mode="time-dependent"
rh41!0="time"
rh41&time=time4
rh41!1="lat"
rh41!2="lon"
ncdf4->RH=rh41
wind41@long_name="wind at the lowest atm level (WIND)"
wind41@units="m/s"
wind41@mode="time-dependent"
wind41!0="time"
wind41&time=time4
wind41!1="lat"
wind41!2="lon"
ncdf4->WIND=wind41
pres41@long_name="surface pressure at the lowest atm level (PSRF)"
pres41@units="Pa"
pres41@mode="time-dependent"
pres41!0="time"
pres41&time=time4
pres41!1="lat"
pres41!2="lon"
ncdf4->PSRF=pres41
prate41@long_name="preciptation(PRECT)"
prate41@units="mm/s"
prate41@mode="time-dependent"
prate41!0="time"
prate41&time=time4
prate41!1="lat"
prate41!2="lon"
ncdf4->PRECTmms=prate41
fsds41@long_name="incident solar(FSDS)"
fsds41@units="W/m2"
fsds41@mode="time-dependent"
fsds41!0="time"
fsds41&time=time4
fsds41!1="lat"
fsds41!2="lon"
ncdf4->FSDS=fsds41
flds41@long_name="incident longwave radiation(FLDS)"
flds41@units="W/m2"
flds41@mode="time-dependent"
flds41!0="time"
flds41&time=time4
flds41!1="lat"
flds41!2="lon"
ncdf4->FLDS=flds41
z3@long_name="observational height"
z3@units="m"
z3@mode="time-invariant"
z3!0="time"
z3&time=time4
z3!1="lat"
z3!2="lon"
ncdf4->ZBOT = z3
ncdf4->LONGXY = lon1
ncdf4->LATIXY = lat1
ncdf4->EDGES = a
ncdf4->EDGEN = b
ncdf4->EDGEE = c
ncdf4->EDGEW = d

tbot51@long_name="temperature at the lowest atm level (TBOT)"
tbot51@units="K"
tbot51@mode="time-dependent"
tbot51!0="time"
tbot51&time=time5
tbot51!1="lat"
tbot51!2="lon"
ncdf5->TBOT=tbot51
rh51@long_name="specific humidity at the lowest atm level (RH)"
rh51@units="kg/kg"
rh51@mode="time-dependent"
rh51!0="time"
rh51&time=time5
rh51!1="lat"
rh51!2="lon"
ncdf5->RH=rh51
wind51@long_name="wind at the lowest atm level (WIND)"
wind51@units="m/s"
wind51@mode="time-dependent"
wind51!0="time"
wind51&time=time5
wind51!1="lat"
wind51!2="lon"
ncdf5->WIND=wind51
pres51@long_name="surface pressure at the lowest atm level (PSRF)"
pres51@units="Pa"
pres51@mode="time-dependent"
pres51!0="time"
pres51&time=time5
pres51!1="lat"
pres51!2="lon"
ncdf5->PSRF=pres51
prate51@long_name="preciptation(PRECT)"
prate51@units="mm/s"
prate51@mode="time-dependent"
prate51!0="time"
prate51&time=time5
prate51!1="lat"
prate51!2="lon"
ncdf5->PRECTmms=prate51
fsds51@long_name="incident solar(FSDS)"
fsds51@units="W/m2"
fsds51@mode="time-dependent"
fsds51!0="time"
fsds51&time=time5
fsds51!1="lat"
fsds51!2="lon"
ncdf5->FSDS=fsds51
flds51@long_name="incident longwave radiation(FLDS)"
flds51@units="W/m2"
flds51@mode="time-dependent"
flds51!0="time"
flds51&time=time5
flds51!1="lat"
flds51!2="lon"
ncdf5->FLDS=flds51
z1@long_name="observational height"
z1@units="m"
z1@mode="time-invariant"
z1!0="time"
z1&time=time5
z1!1="lat"
z1!2="lon"
ncdf5->ZBOT = z1
ncdf5->LONGXY = lon1
ncdf5->LATIXY = lat1
ncdf5->EDGES = a
ncdf5->EDGEN = b
ncdf5->EDGEE = c
ncdf5->EDGEW = d

tbot61@long_name="temperature at the lowest atm level (TBOT)"
tbot61@units="K"
tbot61@mode="time-dependent"
tbot61!0="time"
tbot61&time=time6
tbot61!1="lat"
tbot61!2="lon"
ncdf6->TBOT=tbot61
rh61@long_name="specific humidity at the lowest atm level (RH)"
rh61@units="kg/kg"
rh61@mode="time-dependent"
rh61!0="time"
rh61&time=time6
rh61!1="lat"
rh61!2="lon"
ncdf6->RH=rh61
wind61@long_name="wind at the lowest atm level (WIND)"
wind61@units="m/s"
wind61@mode="time-dependent"
wind61!0="time"
wind61&time=time6
wind61!1="lat"
wind61!2="lon"
ncdf6->WIND=wind61
pres61@long_name="surface pressure at the lowest atm level (PSRF)"
pres61@units="Pa"
pres61@mode="time-dependent"
pres61!0="time"
pres61&time=time6
pres61!1="lat"
pres61!2="lon"
ncdf6->PSRF=pres61
prate61@long_name="preciptation(PRECT)"
prate61@units="mm/s"
prate61@mode="time-dependent"
prate61!0="time"
prate61&time=time6
prate61!1="lat"
prate61!2="lon"
ncdf6->PRECTmms=prate61
fsds61@long_name="incident solar(FSDS)"
fsds61@units="W/m2"
fsds61@mode="time-dependent"
fsds61!0="time"
fsds61&time=time6
fsds61!1="lat"
fsds61!2="lon"
ncdf6->FSDS=fsds61
flds61@long_name="incident longwave radiation(FLDS)"
flds61@units="W/m2"
flds61@mode="time-dependent"
flds61!0="time"
flds61&time=time6
flds61!1="lat"
flds61!2="lon"
ncdf6->FLDS=flds61
z3@long_name="observational height"
z3@units="m"
z3@mode="time-invariant"
z3!0="time"
z3&time=time6
z3!1="lat"
z3!2="lon"
ncdf6->ZBOT = z3
ncdf6->LONGXY = lon1
ncdf6->LATIXY = lat1
ncdf6->EDGES = a
ncdf6->EDGEN = b
ncdf6->EDGEE = c
ncdf6->EDGEW = d

tbot71@long_name="temperature at the lowest atm level (TBOT)"
tbot71@units="K"
tbot71@mode="time-dependent"
tbot71!0="time"
tbot71&time=time7
tbot71!1="lat"
tbot71!2="lon"
ncdf7->TBOT=tbot71
rh71@long_name="specific humidity at the lowest atm level (RH)"
rh71@units="kg/kg"
rh71@mode="time-dependent"
rh71!0="time"
rh71&time=time7
rh71!1="lat"
rh71!2="lon"
ncdf7->RH=rh71
wind71@long_name="wind at the lowest atm level (WIND)"
wind71@units="m/s"
wind71@mode="time-dependent"
wind71!0="time"
wind71&time=time7
wind71!1="lat"
wind71!2="lon"
ncdf7->WIND=wind71
pres71@long_name="surface pressure at the lowest atm level (PSRF)"
pres71@units="Pa"
pres71@mode="time-dependent"
pres71!0="time"
pres71&time=time7
pres71!1="lat"
pres71!2="lon"
ncdf7->PSRF=pres71
prate71@long_name="preciptation(PRECT)"
prate71@units="mm/s"
prate71@mode="time-dependent"
prate71!0="time"
prate71&time=time7
prate71!1="lat"
prate71!2="lon"
ncdf7->PRECTmms=prate71
fsds71@long_name="incident solar(FSDS)"
fsds71@units="W/m2"
fsds71@mode="time-dependent"
fsds71!0="time"
fsds71&time=time7
fsds71!1="lat"
fsds71!2="lon"
ncdf7->FSDS=fsds71
flds71@long_name="incident longwave radiation(FLDS)"
flds71@units="W/m2"
flds71@mode="time-dependent"
flds71!0="time"
flds71&time=time7
flds71!1="lat"
flds71!2="lon"
ncdf7->FLDS=flds71
z1@long_name="observational height"
z1@units="m"
z1@mode="time-invariant"
z1!0="time"
z1&time=time7
z1!1="lat"
z1!2="lon"
ncdf7->ZBOT = z1
ncdf7->LONGXY = lon1
ncdf7->LATIXY = lat1
ncdf7->EDGES = a
ncdf7->EDGEN = b
ncdf7->EDGEE = c
ncdf7->EDGEW = d

tbot81@long_name="temperature at the lowest atm level (TBOT)"
tbot81@units="K"
tbot81@mode="time-dependent"
tbot81!0="time"
tbot81&time=time8
tbot81!1="lat"
tbot81!2="lon"
ncdf8->TBOT=tbot81
rh81@long_name="specific humidity at the lowest atm level (RH)"
rh81@units="kg/kg"
rh81@mode="time-dependent"
rh81!0="time"
rh81&time=time8
rh81!1="lat"
rh81!2="lon"
ncdf8->RH=rh81
wind81@long_name="wind at the lowest atm level (WIND)"
wind81@units="m/s"
wind81@mode="time-dependent"
wind81!0="time"
wind81&time=time8
wind81!1="lat"
wind81!2="lon"
ncdf8->WIND=wind81
pres81@long_name="surface pressure at the lowest atm level (PSRF)"
pres81@units="Pa"
pres81@mode="time-dependent"
pres81!0="time"
pres81&time=time8
pres81!1="lat"
pres81!2="lon"
ncdf8->PSRF=pres81
prate81@long_name="preciptation(PRECT)"
prate81@units="mm/s"
prate81@mode="time-dependent"
prate81!0="time"
prate81&time=time8
prate81!1="lat"
prate81!2="lon"
ncdf8->PRECTmms=prate81
fsds81@long_name="incident solar(FSDS)"
fsds81@units="W/m2"
fsds81@mode="time-dependent"
fsds81!0="time"
fsds81&time=time8
fsds81!1="lat"
fsds81!2="lon"
ncdf8->FSDS=fsds81
flds81@long_name="incident longwave radiation(FLDS)"
flds81@units="W/m2"
flds81@mode="time-dependent"
flds81!0="time"
flds81&time=time8
flds81!1="lat"
flds81!2="lon"
ncdf8->FLDS=flds81
z1@long_name="observational height"
z1@units="m"
z1@mode="time-invariant"
z1!0="time"
z1&time=time8
z1!1="lat"
z1!2="lon"
ncdf8->ZBOT = z1
ncdf8->LONGXY = lon1
ncdf8->LATIXY = lat1
ncdf8->EDGES = a
ncdf8->EDGEN = b
ncdf8->EDGEE = c
ncdf8->EDGEW = d

tbot91@long_name="temperature at the lowest atm level (TBOT)"
tbot91@units="K"
tbot91@mode="time-dependent"
tbot91!0="time"
tbot91&time=time9
tbot91!1="lat"
tbot91!2="lon"
ncdf9->TBOT=tbot91
rh91@long_name="specific humidity at the lowest atm level (RH)"
rh91@units="kg/kg"
rh91@mode="time-dependent"
rh91!0="time"
rh91&time=time9
rh91!1="lat"
rh91!2="lon"
ncdf9->RH=rh91
wind91@long_name="wind at the lowest atm level (WIND)"
wind91@units="m/s"
wind91@mode="time-dependent"
wind91!0="time"
wind91&time=time9
wind91!1="lat"
wind91!2="lon"
ncdf9->WIND=wind91
pres91@long_name="surface pressure at the lowest atm level (PSRF)"
pres91@units="Pa"
pres91@mode="time-dependent"
pres91!0="time"
pres91&time=time9
pres91!1="lat"
pres91!2="lon"
ncdf9->PSRF=pres91
prate91@long_name="preciptation(PRECT)"
prate91@units="mm/s"
prate91@mode="time-dependent"
prate91!0="time"
prate91&time=time9
prate91!1="lat"
prate91!2="lon"
ncdf9->PRECTmms=prate91
fsds91@long_name="incident solar(FSDS)"
fsds91@units="W/m2"
fsds91@mode="time-dependent"
fsds91!0="time"
fsds91&time=time9
fsds91!1="lat"
fsds91!2="lon"
ncdf9->FSDS=fsds91
flds91@long_name="incident longwave radiation(FLDS)"
flds91@units="W/m2"
flds91@mode="time-dependent"
flds91!0="time"
flds91&time=time9
flds91!1="lat"
flds91!2="lon"
ncdf9->FLDS=flds91
z3@long_name="observational height"
z3@units="m"
z3@mode="time-invariant"
z3!0="time"
z3&time=time9
z3!1="lat"
z3!2="lon"
ncdf9->ZBOT = z3
ncdf9->LONGXY = lon1
ncdf9->LATIXY = lat1
ncdf9->EDGES = a
ncdf9->EDGEN = b
ncdf9->EDGEE = c
ncdf9->EDGEW = d

tbot101@long_name="temperature at the lowest atm level (TBOT)"
tbot101@units="K"
tbot101@mode="time-dependent"
tbot101!0="time"
tbot101&time=time10
tbot101!1="lat"
tbot101!2="lon"
ncdf10->TBOT=tbot101
rh101@long_name="specific humidity at the lowest atm level (RH)"
rh101@units="kg/kg"
rh101@mode="time-dependent"
rh101!0="time"
rh101&time=time10
rh101!1="lat"
rh101!2="lon"
ncdf10->RH=rh101
wind101@long_name="wind at the lowest atm level (WIND)"
wind101@units="m/s"
wind101@mode="time-dependent"
wind101!0="time"
wind101&time=time10
wind101!1="lat"
wind101!2="lon"
ncdf10->WIND=wind101
pres101@long_name="surface pressure at the lowest atm level (PSRF)"
pres101@units="Pa"
pres101@mode="time-dependent"
pres101!0="time"
pres101&time=time10
pres101!1="lat"
pres101!2="lon"
ncdf10->PSRF=pres101
prate101@long_name="preciptation(PRECT)"
prate101@units="mm/s"
prate101@mode="time-dependent"
prate101!0="time"
prate101&time=time10
prate101!1="lat"
prate101!2="lon"
ncdf10->PRECTmms=prate101
fsds101@long_name="incident solar(FSDS)"
fsds101@units="W/m2"
fsds101@mode="time-dependent"
fsds101!0="time"
fsds101&time=time10
fsds101!1="lat"
fsds101!2="lon"
ncdf10->FSDS=fsds101
flds101@long_name="incident longwave radiation(FLDS)"
flds101@units="W/m2"
flds101@mode="time-dependent"
flds101!0="time"
flds101&time=time10
flds101!1="lat"
flds101!2="lon"
ncdf10->FLDS=flds101
z1@long_name="observational height"
z1@units="m"
z1@mode="time-invariant"
z1!0="time"
z1&time=time10
z1!1="lat"
z1!2="lon"
ncdf10->ZBOT = z1
ncdf10->LONGXY = lon1
ncdf10->LATIXY = lat1
ncdf10->EDGES = a
ncdf10->EDGEN = b
ncdf10->EDGEE = c
ncdf10->EDGEW = d

tbot111@long_name="temperature at the lowest atm level (TBOT)"
tbot111@units="K"
tbot111@mode="time-dependent"
tbot111!0="time"
tbot111&time=time11
tbot111!1="lat"
tbot111!2="lon"
ncdf11->TBOT=tbot111
rh111@long_name="specific humidity at the lowest atm level (RH)"
rh111@units="kg/kg"
rh111@mode="time-dependent"
rh111!0="time"
rh111&time=time11
rh111!1="lat"
rh111!2="lon"
ncdf11->RH=rh111
wind111@long_name="wind at the lowest atm level (WIND)"
wind111@units="m/s"
wind111@mode="time-dependent"
wind111!0="time"
wind111&time=time11
wind111!1="lat"
wind111!2="lon"
ncdf11->WIND=wind111
pres111@long_name="surface pressure at the lowest atm level (PSRF)"
pres111@units="Pa"
pres111@mode="time-dependent"
pres111!0="time"
pres111&time=time11
pres111!1="lat"
pres111!2="lon"
ncdf11->PSRF=pres111
prate111@long_name="preciptation(PRECT)"
prate111@units="mm/s"
prate111@mode="time-dependent"
prate111!0="time"
prate111&time=time11
prate111!1="lat"
prate111!2="lon"
ncdf11->PRECTmms=prate111
fsds111@long_name="incident solar(FSDS)"
fsds111@units="W/m2"
fsds111@mode="time-dependent"
fsds111!0="time"
fsds111&time=time11
fsds111!1="lat"
fsds111!2="lon"
ncdf11->FSDS=fsds111
flds111@long_name="incident longwave radiation(FLDS)"
flds111@units="W/m2"
flds111@mode="time-dependent"
flds111!0="time"
flds111&time=time11
flds111!1="lat"
flds111!2="lon"
ncdf11->FLDS=flds111
z3@long_name="observational height"
z3@units="m"
z3@mode="time-invariant"
z3!0="time"
z3&time=time11
z3!1="lat"
z3!2="lon"
ncdf11->ZBOT = z3
ncdf11->LONGXY = lon1
ncdf11->LATIXY = lat1
ncdf11->EDGES = a
ncdf11->EDGEN = b
ncdf11->EDGEE = c
ncdf11->EDGEW = d

tbot121@long_name="temperature at the lowest atm level (TBOT)"
tbot121@units="K"
tbot121@mode="time-dependent"
tbot121!0="time"
tbot121&time=time12
tbot121!1="lat"
tbot121!2="lon"
ncdf12->TBOT=tbot121
rh121@long_name="specific humidity at the lowest atm level (RH)"
rh121@units="kg/kg"
rh121@mode="time-dependent"
rh121!0="time"
rh121&time=time12
rh121!1="lat"
rh121!2="lon"
ncdf12->RH=rh121
wind121@long_name="wind at the lowest atm level (WIND)"
wind121@units="m/s"
wind121@mode="time-dependent"
wind121!0="time"
wind121&time=time12
wind121!1="lat"
wind121!2="lon"
ncdf12->WIND=wind121
pres121@long_name="surface pressure at the lowest atm level (PSRF)"
pres121@units="Pa"
pres121@mode="time-dependent"
pres121!0="time"
pres121&time=time12
pres121!1="lat"
pres121!2="lon"
ncdf12->PSRF=pres121
prate121@long_name="preciptation(PRECT)"
prate121@units="mm/s"
prate121@mode="time-dependent"
prate121!0="time"
prate121&time=time12
prate121!1="lat"
prate121!2="lon"
ncdf12->PRECTmms=prate121
fsds121@long_name="incident solar(FSDS)"
fsds121@units="W/m2"
fsds121@mode="time-dependent"
fsds121!0="time"
fsds121&time=time12
fsds121!1="lat"
fsds121!2="lon"
ncdf12->FSDS=fsds121
flds121@long_name="incident longwave radiation(FLDS)"
flds121@units="W/m2"
flds121@mode="time-dependent"
flds121!0="time"
flds121&time=time12
flds121!1="lat"
flds121!2="lon"
ncdf12->FLDS=flds121
z3@long_name="observational height"
z3@units="m"
z3@mode="time-invariant"
z3!0="time"
z3&time=time12
z3!1="lat"
z3!2="lon"
ncdf12->ZBOT = z3
ncdf12->LONGXY = lon1
ncdf12->LATIXY = lat1
ncdf12->EDGES = a
ncdf12->EDGEN = b
ncdf12->EDGEE = c
ncdf12->EDGEW = d

end

jdy

有木有高人回答啊~~~~之前区域的完全能跑  没有问题。。但是单点老提示这个