[RegCNET] interpolation of re-analysis data and SIGMAtoP.f for ICBC
Moetasim
mashfaq at purdue.edu
Sun Apr 1 06:40:07 CEST 2007
For sigma2p problem, add
integer iy,jx,kz,np
parameter(iy=??,jx=??,kx=18)
in the code after include .... line. Set iy and jx according to your domain.
For Geopotential, the use of icbc.f is not straighforward. The ERA40
geopotential is read in era6hour.f but we do not interpolate it to the regcm
grid. It is used for someother variables calulation in ICBC. If you are good in
fortran, you may try looking in era6hour.f and getera40 subroutine in icbc.f .
Moet
Quoting ywyang <ywyang at lzb.ac.cn>:
> Hi,everyone,
> I want to compare the ERA40 data with the PLEV.yyyymmddhh at 700hPa and
> 500hPa.It seems that the ICBC.yyyymmddhh is the interpolation of re-analysis
> data to the model grid,so can I use the ICBC files to compare with the PLEV
> files?
> However,the ICBC files have no the variable of geopotential height.If I
> want get the variable,can I use the icbc.f to interpolate the re-anlysis
> data?And which lines should be changed?Does anyone provide any suggestion ?
>
> Ps:when I do SIGMAtoP.f for ICBC.yyyymmddhh,there are some false
> suggestion.The following is my steps :
> 1.in input, copy tools/SIGMAtoP.f .
> 2.vi SIGMAtoP.f,change 'fn' and 'fz'
> data fn/'ICBC2000070100'/
> data fz/'IC_P2000070100'/
> 3.pgf77 -byteswapio -o SIGMAtoP SIGMAtoP.f
> 4.there are the false suggestion:
> PGFTN-S-0017-Unable to open include file:
> ../../PreProc/ICBC/icbc.param (SIGMAtoP.f: 7)
> PGFTN-S-0038-Symbol, jx, has not been explicitly declared
> (SIGMAtoP.f: 13)
> PGFTN-S-0038-Symbol, iy, has not been explicitly declared
> (SIGMAtoP.f: 13)
> PGFTN-S-0038-Symbol, kz, has not been explicitly declared
> (SIGMAtoP.f: 16)
> PGFTN-S-0050-Adjustable or assumed size array, u, is not a dummy
> argument (SIGMAtoP.f: 34)
> PGFTN-S-0050-Adjustable or assumed size array, v, is not a dummy
> argument (SIGMAtoP.f: 34)
> PGFTN-S-0050-Adjustable or assumed size array, t, is not a dummy
> argument (SIGMAtoP.f: 34)
> PGFTN-S-0050-Adjustable or assumed size array, qv, is not a dummy
> argument (SIGMAtoP.f: 34)
> PGFTN-S-0050-Adjustable or assumed size array, w, is not a dummy
> argument (SIGMAtoP.f: 34)
> PGFTN-S-0050-Adjustable or assumed size array, pp, is not a dummy
> argument (SIGMAtoP.f: 34)
> PGFTN-S-0050-Adjustable or assumed size array, ps, is not a dummy
> argument (SIGMAtoP.f: 34)
> PGFTN-S-0050-Adjustable or assumed size array, tgb, is not a dummy
> argument (SIGMAtoP.f: 34)
> PGFTN-S-0050-Adjustable or assumed size array, ht, is not a dummy
> argument (SIGMAtoP.f: 34)
> PGFTN-S-0050-Adjustable or assumed size array, h, is not a dummy
> argument (SIGMAtoP.f: 34)
> PGFTN-S-0050-Adjustable or assumed size array, satbrt, is not a dummy
> argument (SIGMAtoP.f: 34)
> PGFTN-S-0050-Adjustable or assumed size array, slp1, is not a dummy
> argument (SIGMAtoP.f: 34)
> PGFTN-S-0050-Adjustable or assumed size array, slp2, is not a dummy
> argument (SIGMAtoP.f: 34)
> PGFTN-S-0050-Adjustable or assumed size array, up, is not a dummy
> argument (SIGMAtoP.f: 34)
> PGFTN-S-0050-Adjustable or assumed size array, vp, is not a dummy
> argument (SIGMAtoP.f: 34)
> PGFTN-S-0050-Adjustable or assumed size array, tp, is not a dummy
> argument (SIGMAtoP.f: 34)
> PGFTN-S-0050-Adjustable or assumed size array, hp, is not a dummy
> argument (SIGMAtoP.f: 34)
> PGFTN-S-0050-Adjustable or assumed size array, qp, is not a dummy
> argument (SIGMAtoP.f: 34)
> PGFTN-S-0050-Adjustable or assumed size array, sig, is not a dummy
> argument (SIGMAtoP.f: 34)
> PGFTN-S-0034-Syntax error at or near end of line (SIGMAtoP.f: 34)
> PGFTN-S-0038-Symbol, iomega, has not been explicitly declared
> (SIGMAtoP.f: 100)
> PGFTN-F-0008-Error limit exceeded (SIGMAtoP.f: 100)
>
>
> How to solve the problem?Any suggestion is appreciated.Thanking you in
> advance.
>
> Best Regards
>
>
>
> ¡¡¡¡¡¡¡¡¡¡¡¡¡¡¡¡ywyang
> ¡¡¡¡¡¡¡¡¡¡¡¡¡¡¡¡ywyang at lzb.ac.cn
> ¡¡¡¡¡¡¡¡¡¡¡¡¡¡¡¡¡¡¡¡2007-03-31
> _______________________________________________
> RegCNET mailing list
> RegCNET at lists.ictp.it
> https://lists.ictp.it/mailman/listinfo/regcnet
>
>
More information about the RegCNET
mailing list