; modified for seac4rs ; Lin Zhang, 05/25/2012 ; Modified by Karen Yu 23 April 2013 @create_nested.pro @create_nested_ascii.pro @regridh_unit_area.pro @regridh_sst.pro pro make_ch_wrapper Model = 'GEOS5_47L' Res = 0.25 Tag = '025x03125' preprocess = 0 xr = [-130, -60] yr = [9.7, 60] ;=============================== ; In RETRO_201103/ ;=============================== ; for RETRO VOC emissions retro_2011 = 0 ; done IF ( retro_2011 ) then begin InDir = '/as/data/geos/GEOS_2x2.5/RETRO_201103/' OutDir1 = 'RETRO_201103_global/' OutDir2 = 'RETRO_201103/' nm = 12 for im = 1, nm do begin date = string( 200000L + im, format='(i6)' ) print, 'regrid RETRO ', date infile = date + '.2x25' outfile= date + '.' + Tag if ( ~preprocess ) then $ Regridh_unit_area, infilename = indir + infile, $ OutFilename = outdir1 + outfile, $ OutModelname = Model, OutResolution = Res, diagn = 'ANTHSRCE' Create_Nested, InFileName = OutDir1 + OutFile, $ OutFileName = OutDir2 + OutFile, $ xrange=xr, yrange=yr endfor ENDIF ;=============================== ; In pco_lco_200203/ ;=============================== pco_lco = 0 ; done, need to make sure tracerinfo.dat and diaginfo.dat is not ; in the local directory IF ( pco_lco ) then begin InFileName = 'data/raw/CO.P_L.data' OutDir1 = 'pco_lco_200203_global/' OutDir2 = 'pco_lco_200203/' if ( ~preprocess ) then $ RegridVH_PCO_LCO, InFileName=InFileName, $ OutModelName = Model, OutResolution = Res, OutDir = OutDir1 OutFile = 'COprod.geos5.' + Tag Create_Nested, InFileName = OutDir1 + OutFile, $ OutFileName = OutDir2 + OutFile, $ xrange=xr, yrange=yr OutFile = 'COloss.geos5.' + Tag Create_Nested, InFileName = OutDir1 + OutFile, $ OutFileName = OutDir2 + OutFile, $ xrange=xr, yrange=yr ENDIF ;=============================== ; In stratjv_200203/ ;=============================== stratjv = 1 ; done IF ( stratjv ) then begin InDir = '/as/data/geos/GEOS_2x2.5/stratjv_200203/' OutDir1 = '../stratjv_200203_global/' OutDir2 = './' ; JV-MAP InFile = 'stratjv.geos5.2x25' OutFile = 'stratjv.geos5.' + Tag if ( ~preprocess ) then $ RegridVH_StratJV, InFileName = InDir + InFile, $ OutFileName = OutDir1 + OutFile, $ OutModelName = Model, OutResolution = Res Create_Nested, InFileName = OutDir1 + OutFile, $ OutFileName = OutDir2 + OutFile, $ xrange=xr, yrange=yr ENDIF ;=============================== ; In stratOH_200203/ *** ;=============================== stratOH = 0 ; done IF ( stratOH ) then begin InDir = '/as/data/geos/GEOS_2x2.5/stratOH_200203/' OutDir1 = 'stratOH_200203_global/' OutDir2 = 'stratOH_200203/' ; straOH InFile = 'stratOH.geos5.2x25' OutFile = 'stratOH.geos5.' + Tag if ( ~preprocess ) then $ RegridVH_StratJV, InFileName = InDir + InFile, $ OutFileName = OutDir1 + OutFile, $ OutModelName = Model, OutResolution = Res, DiagN='CHEM-L=$' Create_Nested, InFileName = OutDir1 + OutFile, $ OutFileName = OutDir2 + OutFile, $ xrange=xr, yrange=yr ENDIF ;=============================== ; In pnoy_200106/ ;=============================== ; for NOy upper condition pnoy_200106 = 0 ; done IF ( pnoy_200106 ) then begin Infile = '/home/bmy/archive/data/pnoy_200106/raw/PNOY.data' OutDir1 = 'pnoy_200106_global/' OutDir2 = 'pnoy_200106/' ; 1st file OutFile = 'pnoy_nox_hno3.geos5.' + Tag if ( ~preprocess ) then $ RegridVH_NOy, InFileName = Infile, OutFileName = OutDir1 + OutFile, $ OutModelName = Model, OutResolution = Res Create_Nested, InFileName = OutDir1 + OutFile, $ OutFileName = OutDir2 + OutFile, $ xrange=[-130,-60], yrange=[9.7,60] ENDIF ;=============================== ; In biofuel_200202/ ;=============================== ; for biofuel gases biofuel_2002 = 0 ; don't know where to find this data either. Bob has a .tar.gz file ; should I use it? ; kyu 25 Apr 2013 IF ( biofuel_2002 ) then begin InDir = 'data/1x1_gen/' OutDir1 = 'biofuel_200202_global/' OutDir2 = 'biofuel_200202/' InFile = 'biofuel.generic.1x1' OutFile = 'biofuel.geos.' + Tag if ( ~preprocess ) then $ RegridH_NH3, InFileName = InDir + InFile, $ OutFileName = OutDir1 + OutFile, $ OutModelName = Model, OutResolution = Res Create_Nested, InFileName = OutDir1 + OutFile, $ OutFileName = OutDir2 + OutFile, $ xrange=xr, yrange=yr ENDIF ;=============================== ; In acetone_200108/ ;=============================== ; for JO1D acetone_2001 = 0 ; done IF ( acetone_2001 ) then begin InDir = '/as/data/geos/GEOS_2x2.5/acetone_200108/' OutDir1 = 'acetone_200108_global/' OutDir2 = 'acetone_200108/' ; JO1D InFile = 'JO1D.geos.2x25' OutFile = 'JO1D.geos.' + Tag if ( ~preprocess ) then $ RegridH_JO1D, InFileName = InDir + InFile, $ OutFileName = OutDir1 + OutFile, $ OutModelName = Model, OutResolution = Res Create_Nested, InFileName = OutDir1 + OutFile, $ OutFileName = OutDir2 + OutFile, $ xrange=[-130,-60], yrange=[9.7,60] ; HET-RESP InFile = 'resp.geos.2x25' OutFile = 'resp.geos.' + Tag if ( ~preprocess ) then $ RegridH_JO1D, InFileName = InDir + InFile, $ OutFileName = OutDir1 + OutFile, $ OutModelName = Model, OutResolution = Res, DiagN = 'HET-RESP' Create_Nested, InFileName = OutDir1 + OutFile, $ OutFileName = OutDir2 + OutFile, $ xrange=xr, yrange=yr ENDIF ;=============================== ; In fossil_200104/ ;=============================== ; for GEIA emissions fossil_geia = 0 ; done IF ( fossil_geia ) then begin InFileName='/home/bmy/archive/data/fossil_200104/1x1_gen/merge_nobiofuels.1x1_STD_SASS' OutDir1 = 'fossil_200104_global/' OutDir2 = 'fossil_200104/' OutFile = 'merge_nobiofuels.geos.' + Tag if ( ~preprocess ) then $ RegridH_Anthro, InFileName = InFileName, $ OutModelName = Model, OutResolution = Res, $ OutFileName = OutDir1 + OutFile Create_Nested, InFileName = OutDir1 + OutFile, $ OutFileName = OutDir2 + OutFile, $ xrange=xr, yrange=yr ENDIF ;=============================== ; In C3H8_C2H6_200109/ ;=============================== ; for C3H8 C2H6 gases c3h8_c2h6 = 0 ; done IF ( c3h8_c2h6 ) then begin InDir = '/as/data/geos/GEOS_2x2.5/C3H8_C2H6_200109/' OutDir1 = 'C3H8_C2H6_200109_global/' OutDir2 = 'C3H8_C2H6_200109/' File = 'C3H8_C2H6_ngas.geos.' OutFile = File + Tag if ( ~preprocess ) then $ RegridH_C3H8_C2H6, InFileName = InDir+File+'2x25', $ OutFileName = OutDir1 + OutFile, $ OutModelName = Model, OutResolution = Res Create_Nested, InFileName = OutDir1 + OutFile, $ OutFileName = OutDir2 + OutFile, $ xrange=xr, yrange=yr ENDIF ;=============================== ; In carbon_200909/ ;=============================== ; for T. Bonds OC/BC carbon_200909 = 0 ; done IF ( carbon_200909 ) then begin InDir = '/as/data/geos/GEOS_1x1/carbon_200909/' OutDir1 = 'carbon_200909_global/' OutDir2 = 'carbon_200909/' ; 1st file File = 'BCOC_TBond_biofuel.2000.' InFile = File + 'generic.1x1' OutFile = File + 'geos.' + Tag if ( ~preprocess ) then $ REGRIDH_NH3, INFILENAME = indir + InFile, $ OUTFILENAME = OutDir1 + OutFile, $ OutModelName = Model, OutResolution = Res Create_Nested, InFileName = OutDir1 + OutFile, $ OutFileName = OutDir2 + OutFile, $ xrange=xr, yrange=yr ; 2nd file File = 'BCOC_TBond_fossil.2000.' InFile = File + 'generic.1x1' OutFile = File + 'geos.' + Tag if ( ~preprocess ) then $ REGRIDH_NH3, INFILENAME = indir + InFile, $ OUTFILENAME = OutDir1 + OutFile, $ OutModelName = Model, OutResolution = Res Create_Nested, InFileName = OutDir1 + OutFile, $ OutFileName = OutDir2 + OutFile, $ xrange=xr, yrange=yr ENDIF ;=============================== ; In dust_200605/ ;=============================== dust_200605 = 0 ; done if ( dust_200605 ) then begin InDir = '/as/data/geos/GEOS_2x2.5/dust_200605/' OutDir1 = 'dust_200605_global/' OutDir2 = 'dust_200605/' ; 1st file if ( ~preprocess ) then $ REGRIDH_jv, INFILENAME = Indir + 'GOCART_src_fn.geos.2x25', $ OUTFILENAME = OutDir1 + 'GOCART_src_fn.geos.' + Tag, $ outmodelname = Model, OUTRESOLUTION = Res,diagN='DEAD-2D' Create_Nested, InFileName = OutDir1 + 'GOCART_src_fn.geos.' + Tag, $ OutFileName = OutDir2 + 'GOCART_src_fn.geos.' + Tag, $ xrange=xr, yrange=yr ; 2nd file if ( ~preprocess ) then $ REGRIDH_jv, INFILENAME = Indir + 'dst_tibds.geos.2x25', $ OUTFILENAME = OutDir1 + 'dst_tibds.geos.' + Tag, $ outmodelname = Model, OUTRESOLUTION = Res, diagN='DEAD-2D' Create_Nested, InFileName = OutDir1 + 'dst_tibds.geos.' + Tag, $ OutFileName = OutDir2 + 'dst_tibds.geos.' + Tag, $ xrange=xr, yrange=yr ; 3rd file if ( ~preprocess ) then $ REGRIDH_jv, INFILENAME = Indir + 'dst_tvbds.geos.2x25', $ OUTFILENAME = OutDir1 + 'dst_tvbds.geos.' + Tag, $ outmodelname = Model, OUTRESOLUTION = Res, diagN='DEAD-2D' Create_Nested, InFileName = OutDir1 + 'dst_tvbds.geos.' + Tag, $ OutFileName = OutDir2 + 'dst_tvbds.geos.' + Tag, $ xrange=xr, yrange=yr endif ;=============================== ; In uvalbedo_200111/ ;=============================== uvalbedo_200111 = 0 ; done if ( uvalbedo_200111 ) then begin OutDir1 = 'uvalbedo_200111_global/' OutDir2 = 'uvalbedo_200111/' File = 'uvalbedo.geos.' + Tag if ( ~preprocess ) then $ regridh_uvalbedo, outmodelname = Model, Outresolution=Res, $ OutFileName=OutDir1+File Create_Nested, InFileName = OutDir1 + File, $ OutFileName = OutDir2 + File, $ xrange=xr, yrange=yr endif ;=============================== ; In MEGAN_200510/ ;=============================== MEGAN_200510 = 0 ; done if ( MEGAN_200510 ) then begin preprocess = 1 ; already have files at the resolution we want InDir = '/as/data/geos/GEOS_NATIVE/MEGAN_201211/' ; OutDir1 = 'MEGAN_201211_global/' OutDir2 = 'MEGAN_201211/' ; 1st file if ( ~preprocess ) then $ regridh_biomass,infilename = InDir + 'MEGAN_AEF_ISOP.geos.1x1', $ outfilename = OutDir1 + 'MEGAN_AEF_ISOP.geos.' + Tag, $ outmodelname = Model, outresolution = Res, diagN='BIOGSRCE' Create_Nested, InFileName = InDir + 'CLM4_PFT.geos.' + Tag, $ OutFileName = OutDir2 + 'CLM4_PFT.geos.' + Tag, $ xrange=xr, yrange=yr ; 2nd file if ( ~preprocess ) then $ regridh_biomass,infilename = InDir + 'MEGAN_AEF_MBO.geos.1x1', $ outfilename = OutDir1 + 'MEGAN_AEF_MBO.geos.' + Tag, $ outmodelname = Model, outresolution = Res, diagN='BIOGSRCE' Create_Nested, InFileName = InDir + 'MEGAN2.1_EF_APIN_ug.m2.h_geos.' + Tag, $ OutFileName = OutDir2 + 'MEGAN2.1_EF_APIN_ug.m2.h_geos.' + Tag, $ xrange=xr, yrange=yr ; 3rd file if ( ~preprocess ) then $ regridh_biomass,infilename = InDir + 'MEGAN_AEF_MTP.geos.1x1', $ outfilename = OutDir1 + 'MEGAN_AEF_MTP.geos.' + Tag, $ outmodelname = Model, outresolution = Res, diagN='BIOGSRCE' Create_Nested, InFileName = InDir + 'MEGAN2.1_EF_BPIN_ug.m2.h_geos.' + Tag, $ OutFileName = OutDir2 + 'MEGAN2.1_EF_BPIN_ug.m2.h_geos.' + Tag, $ xrange=xr, yrange=yr ; 4th file if ( ~preprocess ) then $ regridh_biomass,infilename = InDir + 'MEGAN_AEF_OVC.geos.1x1', $ outfilename = OutDir1 + 'MEGAN_AEF_OVC.geos.' + Tag, $ outmodelname = Model, outresolution = Res, diagN='BIOGSRCE' Create_Nested, InFileName = InDir + 'MEGAN2.1_EF_CARE_ug.m2.h_geos.' + Tag, $ OutFileName = OutDir2 + 'MEGAN2.1_EF_CARE_ug.m2.h_geos.' + Tag, $ xrange=xr, yrange=yr ; 5th file if ( ~preprocess ) then $ regridh_biomass,infilename = InDir + 'MEGAN_AEF_ISOP.updated_200708.geos.1x1', $ outfilename = OutDir1 + 'MEGAN_AEF_ISOP.updated_200708.geos.' + Tag, $ outmodelname = Model, outresolution = Res, diagN='BIOGSRCE' Create_Nested, InFileName = InDir + 'MEGAN2.1_EF_ISOP_ug.m2.h_geos.' + Tag, $ OutFileName = OutDir2 + 'MEGAN2.1_EF_ISOP_ug.m2.h_geos.' + Tag, $ xrange=xr, yrange=yr ; 6th file Create_Nested, InFileName = InDir + 'MEGAN2.1_EF_LIMO_ug.m2.h_geos.' + Tag, $ OutFileName = OutDir2 + 'MEGAN2.1_EF_LIMO_ug.m2.h_geos.' + Tag, $ xrange=xr, yrange=yr ; 7th file Create_Nested, InFileName = InDir + 'MEGAN2.1_EF_MBO_ug.m2.h_geos.' + Tag, $ OutFileName = OutDir2 + 'MEGAN2.1_EF_MBO_ug.m2.h_geos.' + Tag, $ xrange=xr, yrange=yr ; 8th file Create_Nested, InFileName = InDir + 'MEGAN2.1_EF_MYRC_ug.m2.h_geos.' + Tag, $ OutFileName = OutDir2 + 'MEGAN2.1_EF_MYRC_ug.m2.h_geos.' + Tag, $ xrange=xr, yrange=yr ; 9th file Create_Nested, InFileName = InDir + 'MEGAN2.1_EF_OCIM_ug.m2.h_geos.' + Tag, $ OutFileName = OutDir2 + 'MEGAN2.1_EF_OCIM_ug.m2.h_geos.' + Tag, $ xrange=xr, yrange=yr ; 10th file Create_Nested, InFileName = InDir + 'MEGAN2.1_EF_SABI_ug.m2.h_geos.' + Tag, $ OutFileName = OutDir2 + 'MEGAN2.1_EF_SABI_ug.m2.h_geos.' + Tag, $ xrange=xr, yrange=yr preprocess = 0 endif ;=============================== ; In TOMS_201203/ ;=============================== TOMS_201203 = 0 IF ( TOMS_201203 ) then begin InDir = '/as/data/geos/GEOS_1x1/TOMS_201203/' OutDir1 = 'TOMS_201203_global/' OutDir2 = 'TOMS_201203/' ; loop for years syear = 2006 eyear = 2010 nyr = eyear - syear + 1 For iy = syear, eyear do begin File = 'TOMS_O3col_' + string(iy,format='(i4)') + '.geos.' print, 'processing ', file if ( ~preprocess ) then $ Regrid_toms_sbuv_seac4rs, INFILENAME = Indir + File + '1x1', $ OUTFILENAME = OutDir1 + File + Tag, $ outmodelname= Model, OUTRESOLUTION = Res Create_Nested, InFileName = OutDir1 + File + Tag, $ OutFileName = OutDir2 + File + Tag, $ xrange=xr, yrange=yr Endfor Endif ;=============================== ; In sulfate_sim_200508/ ;=============================== ; for sulfate simulation sulfate_sim = 0 ; done IF ( sulfate_sim ) then begin preprocess = 0 InDir = '/as/data/geos/GEOS_2x2.5/sulfate_sim_200508/' OutDir1 = 'sulfate_sim_200508_global/' OutDir2 = 'sulfate_sim_200508/' DMS = 0 IF ( DMS ) then begin ;------------------------------ ; DMS seawater emissions ;------------------------------ if ( ~preprocess ) then $ REGRIDH_UNIT_AREA, infilename = indir + 'DMS_seawater.geos.2x25', $ outfilename = outdir1 + 'DMS_seawater.geos.' + Tag, $ outmodelname=Model, OUTRESOLUTION = Res, $ diagn = 'BIOGSRCE' Create_Nested, InFileName = OutDir1 + 'DMS_seawater.geos.' + Tag, $ OutFileName = OutDir2 + 'DMS_seawater.geos.' + Tag, $ xrange=xr, yrange=yr ENDIF NH3 = 0 IF ( NH3 ) then begin ;------------------------------ ; NH3 emissions ;------------------------------ if ( ~preprocess ) then $ REGRIDH_NH3, INFILENAME = Indir + 'NH3_anthsrce.geos.2x25', $ OUTFILENAME = OutDir1 + 'NH3_anthsrce.geos.' + Tag, $ outmodelname= Model, OUTRESOLUTION = Res Create_Nested, InFileName = OutDir1 + 'NH3_anthsrce.geos.' + Tag, $ OutFileName = OutDir2 + 'NH3_anthsrce.geos.' + Tag, $ xrange=xr, yrange=yr if ( ~preprocess ) then $ REGRIDH_NH3, INFILENAME = Indir + 'NH3_biofuel.geos.2x25', $ OUTFILENAME = OutDir1 + 'NH3_biofuel.geos.' + Tag, $ outmodelname= Model, OUTRESOLUTION = Res Create_Nested, InFileName = OutDir1 + 'NH3_biofuel.geos.' + Tag, $ OutFileName = OutDir2 + 'NH3_biofuel.geos.' + Tag, $ xrange=xr, yrange=yr if ( ~preprocess ) then $ REGRIDH_NH3, INFILENAME = Indir + 'NH3_natusrce.geos.2x25', $ OUTFILENAME = OutDir1 + 'NH3_natusrce.geos.' + Tag, $ outmodelname = Model, OUTRESOLUTION = Res Create_Nested, InFileName = OutDir1 + 'NH3_natusrce.geos.' + Tag, $ OutFileName = OutDir2 + 'NH3_natusrce.geos.' + Tag, $ xrange=xr, yrange=yr if ( ~preprocess ) then $ REGRIDH_NH3, INFILENAME = Indir + 'NH3emis.monthly.geos.2x25', $ OUTFILENAME = OutDir1 + 'NH3emis.monthly.geos.' + Tag, $ outmodelname= Model, OUTRESOLUTION = Res Create_Nested, InFileName = OutDir1 + 'NH3emis.monthly.geos.' + Tag, $ OutFileName = OutDir2 + 'NH3emis.monthly.geos.' + Tag, $ xrange=xr, yrange=yr ENDIF SST = 0 IF ( SST ) then begin ;------------------------------ ; SST and SOx ;------------------------------ if ( ~preprocess ) then $ RegridH_SST, InFileName = Indir + 'SST.geos.2x25', $ OutFileName = OutDir1 + 'SST.geos.' + Tag, $ OutModelName = Model, OutResolution = Res, $ DiagN='DAO-FLDS', Tracer=1305 Create_Nested, InFileName = OutDir1 + 'SST.geos.' + Tag, $ OutFileName = OutDir2 + 'SST.geos.' + Tag, $ xrange=xr, yrange=yr if ( ~preprocess ) then $ REGRIDH_UNIT_AREA, infilename = indir + 'shipSOx.geos.2x25', $ outfilename = outdir1 + 'shipSOx.geos.' + Tag, $ outmodelname = Model, OUTRESOLUTION = Res, diagn = 'SO2-SHIP', tracer = 26 Create_Nested, InFileName = OutDir1 + 'shipSOx.geos.' + Tag, $ OutFileName = OutDir2 + 'shipSOx.geos.' + Tag, $ xrange=xr, yrange=yr ENDIF ENDIF ;=============================== ; In aircraft_NOx_200202 and sulfate_sim_200508/ ;=============================== ; aircraft data aircraft = 0 NOx = 1 SOx = 0 IF ( aircraft ) then begin ; XRange = [70, 140] & YRange = [15, 55] Months = [ 'jan', 'feb', 'mar', 'apr', 'may', 'jun', $ 'jul', 'aug', 'sep', 'oct', 'nov', 'dec' ] ;==================== ; NOx ;==================== IF ( NOx ) then begin ; InDir = '/home/lzh/seac4rs/nested/GEOS_025x03125/aircraft_NOx_200202/' InDir = 'aircraft_200202_global/' OutDir = 'aircraft_NOx_200202/' For M = 0, n_elements(Months)-1 do begin RegridH_aircraft_data, /NOx, OutModelName=Model, OutResolution=Res CREATE_NESTED_ASCII, InFilename = Indir + 'air' + Months[M] + '.025x03125', $ OutFilename = Outdir + 'air' + Months[M] + '.025x03125', $ InResolution =Res, InModelname=Model, $ Format='(2i4, a)', $ xrange = xr, yrange = yr Endfor ENDIF IF ( SOx ) then begin ;==================== ; SOx ;==================== InDir = 'sulfate_sim_200508_global/' OutDir = 'sulfate_sim_200508/' For M = 0, n_elements(Months)-1 do begin ; RegridH_aircraft_data, /SOx, OutModelName=Model, OutResolution=Res CREATE_NESTED_ASCII, InFilename = Indir + 'aircraft.025x03125.1992.' + Months[M], $ OutFilename = Outdir + 'aircraft.025x03125.1992.' + Months[M], $ InResolution =Res, InModelname=Model, $ Header=2, Format='(2i4, a)', $ xrange = xrange, yrange = yrange Endfor ENDIF ENDIF ;=============================== ; In soil_NOx_200203 ;=============================== soil_NOx = 0 IF ( soil_NOx ) then begin ENDIF end