#################### PRE-PROCESSING ##################### #***************************************************** #1# /clima_nfs/nmmb_regPEJA/fix * #***************************************************** * ovde se radi preprocesing za fix polja (topo, landuse, seamask, albedo, vegetacija, vertikalna koordinata, co2 ...) * pusti se runfixed <--- uzima podatke iz geodata, racuna, a zatim upisuje u dir output (u nasem slucaju outnmmb) * lmimjm.inc <--- ovde odredjujemo (tlm0d, tph0d,wbd, sbd, pt,ptsgm, dlmd, dphd) * modelgrid.inc <--- isto kao i lmimjm.inc, definisemo oblast integracije, rezoluciju ... * ako nesto menjamo u modelu(rezoluciju, oblast integracije i slicno), to prvo prepravimo u prethodna dva inc-a #***************************************************** #2# /data/clima/oper/nmme4/ll2tll.ecmwf * #***************************************************** * ovde se nalazi ECM.ctl fajl * ecmwf2nmmbUPP.f <--- program koji: - cita grib (globalna polja), dekodira podatke - pravi fajlove(llspl, llsmst, llstmp, llgsst, llgsno, llgcic) u direktorijumu outnmmb - u njih upisuje (gsmst, gstmp, gsst, gsnow, gcice) - u njemu definisemo oblast, p nivoe koje zelimo, duzinu integracije... #***************************************************** #3# /data/clima/nmmb_regPEJA/vrbl * #***************************************************** * Ukoliko radimo sa GFS uslovima, prvo treba videti da li su prekopirane fajle iz direktorijuma GFS - fajlovi gfs.t00z... treba da budu i u GFS i u vrbl direktorijumu - u GFS su pre bili bocni i pocetni uslovi, sada radimo sa ECMWF ! * llgrid_midres.inc <--- pocetni grid(globalni model) - u njemu su definisane granice integracije, rezolucija modela, deltat i broj vertikalnih nivoa globalnog modela * startovenjem shella runprepreg.sh(cita iz dir output i upisuje u isti taj), pokrecu se 3 stvari: - runalbedo - runvegfrac - run_midres * run_midres <--- ovaj shell vrsi kompajliranje programa allprep_midres.f i njegovo pokretanje #********** allprep_midres.f start ********* * allprep_midres.f <--- sluzi za interpolaciju pocetnih polja i granicnih vrednosti sa globalne na modelsku mrezu tacaka * program grads2anec.f (koji se nalazi u /data/clima/oper/nmme4/ll2tll.ecmwf) vec napravi neka polja i upise ih u direktorijum output. Pusta se sa yymmhh00 (u nasem slucaju grads2anec_ECM_UPP.f) * allprep_midres.f, pored ostalih polja cita neka od ovih koje je napravio grads2anec.f, racuna sa podacima tj. vrsi interpolaciju globalnih polja u modelsku mrezu tacaka i pravi dodatne datoteke u outputu (boco.$ihr, llspl.$ihr) - imi, jmi <--- modelske tacke - imll, jmll <--- podaci iz ECMWF - gblprep <--- sabrutina za interpolaciju modelskih polja na p nivoe, zove sabrutine gtlli i p2hyb (2 puta zove, prvi put za analize, drugi put za granicne uslove) - ihrend, ihboco, LM(nalaze se u rutini gblprep) !!! <--- rucno se menja kraj integracionog perioda, vremenski interval za bocne granicne uslove i broj vertikalnih nivoa - p2hyb <--- prebacuje globalna polja(ECWMF) iz p nivoa u hibridne(ovde smo definisali p nivoe koje mi zelimo) #********** allprep_midres.f end ********* * ako smo menjali rezoluciju, oblast integracije i slicno, moramo i ovde to da promenimo u lmimjm.inc i modelgrid.inc (promenimo u direktorijumu fix, a zatim prekopiramo ove inc-ove u vrbl) * iskljucili smo runcnv_midres, a koristimo grads2anec.f #################### STARTOVANJE MODELA ####################### #***************************************************** #4# /data/clima/nmmb_regPEJA/nmmbrun * #***************************************************** * configure_file.nml <--- setup file za model - definisana rezolucija, broj vertikalnih nivoa, vreme starta modela, period integracije, frekvencija izlaznih podataka, parametrizacione seme... * bsub_cron <--- sluzi za pustanje run-a, napravi lokalni direktoirijum RUN i u njega iskopira sve sto mu je potrebno za rad. - /data/clima/nmmb_regPEJA/nmmbrun/run <--- output modela se ovde upisuje, fajlovi nmmb_hst_01_bin_HHHHh_00m_00.00s #################### POST-PROCESSING ####################### #***************************************************** #5# /data/clima/nmmb_regPEJA/post * #***************************************************** * runpost <--- shell koji kompajlira, a zatim pusta program postall.f - u ovom shellu treba podesiti DATADIR i SRCDIR 1) DATADIR=/data/clima/nmmb_regPEJA/nmmbrun/run <-- direktorijum iz koga cita podatke 2) SRCDIR=/data/clima/nmmb_regPEJA/post * postall.f <--- prebacuje iz hibridnog u p koordinatni sistem, cita nmmb_hst_01_bin_HHHHh_00m_00.00s fajlove iz run direktorijuma i isto u njemu, uspisuje poutHHH fajlove * lmimjm.inc <--- treba da bude isti kao i u direktorijumima fix i vrbl (prekopirati iz fix-a) * outll.inc <--- definisan broj tacaka u oba pravca, rezolucija i granice * vrbldata <--- definisani: period za koji radimo postproc, promenljive, p nivoi i konturni interval 1) ihrstp - prvi trenutak prognoze od koga nadalje zelimo da radimo postproc 2) kount - ukupan broj termina koji posmatramo 3) ninc - vremenski interval izmedju termina 4) var - promenljiva (upisuje se u formi KOD-a) 5) lev - p nivo koji na kome posmatramo promenljivu 6 cntr - konturni interval #***************************************************** #6# /data/clima/nmmb_regPEJA/north * #***************************************************** * rungrads <--- shell koji od poutHHH fajlova pravi poutHHH.grads i pout.ctl fajlove * lmimjm.inc <--- treba da bude isti kao i u prethodnim direktorijumima (prekopirati iz fix-a) * vrbldata <--- treba da bude isti kao i u direktorijumu post # OSTALE NAPOMENE: * Problem je bio sa vlagom u tlu. Vlaga koja je osmotrena(tj iz ECMWF) je veca od naseg qsat( max moguca kolicina vlage ... kod ECM je preko 0.65, a kod nas je ispod 0.5) * sst je problematicna za interpolaciju kod nedefinisanih tacaka u moru (sea_mask) /data/clima/model64/togradsnew.lib/togradsB.f !******************************* inc fajle!!!! u wdir-u modela promeniti .inc fajle: 1. lmimjm.inc (oblast modela) 2. llgrid_midres.inc (oblast globalnih polja) 3. modelgrid.inc (oblast modela - isto kao lmimjm.inc) ./prep.sh kopira .inc fajle tamo gde treba. !********************************* # cvetex et al. :)