Muutokset MLAMIA varten ryhmiteltynä aliohjamakirjastojen mukaan: CLIM: ---- clim/G3COLD.f,v: G=9.80665 clim/G3COLD.f,v: FLD(IPS+JL)=100000. ???? clim/G3COLD.f,v: FLD(IPS+JL)=100000. ????? clim/PCNCART.f,v: G =9.80665 !ACC OF EARTH clim/PCORFIT.f,v: G=9.80665 clim/PCORFIT.f,v: A=6371000. clim/PCSNOW.f,v: G=9.80665 !ACC OF EARTH clim/PCTSBLN.f,v: G =9.80665 ! clim/PCTSTAR.f,v: G=9.80665 clim/TOPPO.f,v: orostd(i,j) = sqrt(orostd(i,j))*9.80665 clim/TOPPO.f,v: 500 ORO(I,J) = ORO(I,J)*9.80665 clim/TOPPO.f,v: mean = 1000000. clim/TOPPO.f,v: mean = 1000000. clim/CRDDRS.f,v: IDDR20(IF0+211)=1013 clim/INTUSN.f,v: ZREA=6371000. GCOD: ---- gcod/GRIBEX.f,v:C (6378.160km,6356.775km,f=1/297.0) oletukset pallosta, ei tehty mitään GRDY: ---- grdy/STATIS.f,v: ZRG = 1.0/9.8066 grdy/STATIS.f,v: ZRG = 1./9.8066 grdy/STATIS.f,v: ZRG = 1./9.8066 grdy/STATIS.f,v: ZRGASD = 287.04 grdy/STATIS.f,v: ZRGASD = 287.04 grdy/CHFRAME.f,v: if(ppsz(jx,jy).eq.rmdi) ppsz(jx,jy)=101300. grdy/DIFINI.f,v: ZLNPRS=ALOG( 101320. ) !CRM UNDEFINED grdy/DIFINI.f,v: RPRERS = 101320. grdy/DIFINI.f,v: ZLNPRS = ALOG( 101320. ) grdy/SLDYNM.f,v: ETAHALF(JK) = AHYB(JK)/101326.0 + BHYB(JK) grdy/SLDYNM.f,v: ETAHALF(K) = AHYB(K)/101326.0 + BHYB(K) grdy/SLDYNM.f,v: PPX(JX,JY,JK) = PLNPSZ(JX,JY) - log(100000.) grdy/SLDYNM.f,v: PLNPSP(JX,JY) = PLNPSP(JX,JY) + log(100000.) grdy/SLDYNM.f,v: @Subtract 1000 hPa to increase accuracy grdy/SLDYNM.f,v: PPX(JX,JY,JK) = PLNPSZ(JX,JY) - log(100000.) grdy/SLDYNM.f,v: PLNPSP(JX,JY) = PLNPSP(JX,JY) + log(100000.) grdy/SLDYNM.f,v: ZRGASD = 287.04 grdy/GEMINI.f,v: INTL = NEXPHL - INTF*1000000 - INTS*1000 ????? grdy/GEMINI.f,v: INTF = NEXPHL/1000000 ? grdy/EDCOMP.f,v: ZDEK = ZDAK/101325.0 + ZDBK grdy/GETDAT.f,v: HYBI(JK) = AHYB(JK)/101325.0 + BHYB(JK) grdy/DELNMI.f,v: ZRGASD = 287.04 grdy/DYN.f,v: ZRGASD = 287.04 grdy/EXPADJ.f,v: ZRGASD = 287.04 grdy/IMPADJ.f,v: ZRGASD = 287.04 grdy/IMPINI.f,v: ZRGASD = 287.04 grdy/SLDYN.f,v: ZRGASD = 287.04 grdy/SLEXPA.f,v: ZRGASD = 287.04 grdy/DIFINI.f,v: DATA ZPI,ZR/3.1415927,6371220./ grdy/OMCOMP.f,v: DATA AEARTH/6371000./ GRW1: ---- grw1/AS2DDR.f,v: RLEVHL( 1)= 101325.000000000 grw1/BLOCK1.f,v: KDEF(13) = NDTBHL / 1000000 ????????? grw1/ACCU.f,v:C REDUCED ACCURACY FOR FIELDS ABOVE 10000 PA grw1/ACCU.f,v: IF(PRESS.GT.10000.)THEN NASU: ---- nasu/HORINT.f,v: G=9.8066 nasu/TQ2DIA.f,v: G=9.8066 nasu/TQ2DIA.f,v: RD=287.049 nasu/TSRED.f,v: PARAMETER (G=9.8) nasu/VERINT.f,v: G=9.8066 nasu/VERINT.f,v: RD=287.049 nasu/FGINT.f,v: REFPRE= 101325.00 nasu/AOFOBS.f,v: NCONV=100000 ???????? nasu/AOFOBS.f,v: NCONV=100000 ???? nasu/VV2MET.f,v: DATA KOD/25,100,300,700,1500,3000,6000,15000,30000,100000/ ????? nasu/VV2MET.f,v: UT = 100000 ???????? nasu/VV2MET.f,v: DATA KOD/25,100,300,700,1500,3000,6000,15000,30000,100000/ nasu/VV2MET.f,v: UT = 100000 ???? nasu/NASURF.f,v: &, A=6378160.0 nasu/RELHUM.f,v: C -- 5418.118466 = 0.622*2.5E+6/287. PHYS: ---- phys/CONPHYS.f,v: GRAVIT = 9.80665 phys/SUCST.f,v: RV=1000.*R/RMV ???? phys/SUCST.f,v:C RA=6371229. phys/SUCST.f,v:C RG=9.80665 phys/CLOUDCV.f,v: zp0=101325. phys/POSTQV.f,v: ZHYBI(JK) = PAHYB(JK)/101325.0 + PBHYB(JK) phys/RADIA.f,v: zpref=101300. phys/RADIA.f,v: cpsref=100000. phys/SUPHEC.f,v: VP00=101320 phys/VDIFF.f,v: zp0 = 101325. phys/CUMASTR.f,v: ZMF2=ZHFL/MAX(ZDH,1.E5*ZDQMIN) phys/INICONDS.f,v: hp0 = 1.e5 phys/INICONS.f,v: hp0 = 1.e5 phys/PBLHGT.f,v: zp0=1.e5 phys/VCBR.f,v: zp0=1.e5 phys/VDIFNL.f,v: zp0=1.e5 phys/CLOUD.f,v: cw(jl,jk)=gravit*50./(dpf(jl,jk)*1000.) phys/CONDCV.f,v: zpcr=10000. phys/CONTUN.f,v: HGTMIX = 1000. ?? phys/ESAT.inc,v: + ,zc5p=10000. phys/ESAT.inc,v: + zc5p=10000. phys/STPHYS.f,v: zph(11) = phstat_l(11,16)*anorm*1000. ??? phys/STPHYS.f,v: zph(12) = phstat_l(12,16)*anorm*1000. ??? PRPO: ---- prpo/INTPRE.f,v: fxtr(jhor) = fxtr(jhor) / 9.80665 prpo/INTPRE.f,v: FXTR(JHOR) = FXTR(JHOR) / 9.80665 prpo/INTPRE.f,v: fxtr(jhor) = fxtr(jhor) / 9.80665 prpo/POSTPP.f,v: FSTORE(JHOR) = FSTORE(JHOR)/9.80665 prpo/POSTPP.f,v: DATA REFPRE /101350.0/ prpo/POSTPP.f,v: ZDHIGH=100000. prpo/POSTPP.f,v: ZDMEAN=100000. prpo/POSTPP.f,v: ZDP=ABS(AF(JK)+BF(JK)*100000.-ZPHIGH) prpo/POSTPP.f,v: ZDP=ABS(AF(JK)+BF(JK)*100000.-ZPMEAN) prpo/POSTPP.f,v: ZDP=ABS(AF(JLEV)+BF(JLEV)*100000.-ZPHIGH) prpo/POSTPP.f,v: ZDP=ABS(AF(JLEV)+BF(JLEV)*100000.-ZPMEAN) prpo/TQUV.f,v: G=9.8066 prpo/TQUV.f,v: RD=287.049 prpo/VINETA.f,v: DATA REFPRE /101325.0/ prpo/VINETA.f,v: ZPREF=100000. prpo/FLTRORO.f,v: PARAMETER (JPLF=1000000) ???? prpo/FLTRORO.f,v: PARAMETER (JPLF=1000000) ???? prpo/INITSF.f,v: isle( 2) = 010001 ! Surface (geopotential) height (m) SCRIPTS: ------- scripts/Forecast,v :SIP0=101320.0 scripts/Forecast,v : sip0=101320.0, scripts/Forecast,v :SIP0=101320.0 scripts/Forecast,v: alevmlp=50000.,100000., scripts/Anaveri,v: MCURRL=100000, ???? scripts/Compare,v: verlev='100000.,85000.,70000.,50000.,25000.' scripts/Postpp,v: 50000.0,70000.0,85000.0,92500.0,100000.0, VARI: ---- vari/GB2LPA.f,v: ZLAC=6371000.*(1.+ABS(SIN(PTHE*ZRHO))) vari/GB2LPC.f,v: ZLAC=6371000.*(1.+ABS(SIN(PTHE*ZRHO))) vari/GB2PLA.f,v: ZLAC=6371000.*(1.+ABS(SIN(PTHE*ZRHO))) vari/GB2PLC.f,v: ZLAC=6371000.*(1.+ABS(SIN(PTHE*ZRHO))) vari/GB2PPA.f,v: ZSCALE=-ZLAC2*ZLAC1*6371000.**2 vari/GB2PPC.f,v: ZSCALE=-ZLAC2*ZLAC1*6371000.**2