Difference between revisions of "Solution to UKCA Chemistry and Aerosol Tutorial 5 Task 5.2"
From UKCA
(Blanked the page) |
|||
| Line 1: | Line 1: | ||
| − | [[UKCA Chemistry and Aerosol Tutorials | Back to UKCA Chemistry and Aerosol Tutorials]] |
||
| − | |||
| − | [[UKCA Chemistry and Aerosol Tutorial 5#Solution to Task_5.2: Add these two new tracers to UKCA | Back to the adding new chemical tracers tutorial]] |
||
| − | |||
| − | ==Task== |
||
| − | |||
| − | You were asked to |
||
| − | |||
| − | <blockquote> |
||
| − | Make the required code changes so that your '''ALICE''' and '''BOB''' tracers are now specified in the UKCA CheST/StratTrop scheme. You should set the conversion factor for each of these to 1.0. |
||
| − | </blockquote> |
||
| − | |||
| − | ==Solution== |
||
| − | |||
| − | The solution to this task involves making changes to 4 files in the UKCA sub-directory of your source code. These files are |
||
| − | |||
| − | ukca_setd1defs.F90 |
||
| − | ukca_chem_strattrop.F90 |
||
| − | ukca_constants.F90 |
||
| − | ukca_cspecies.F90 |
||
| − | |||
| − | The following changes need to be made |
||
| − | |||
| − | ===Changes to ukca_setd1defs.F90=== |
||
| − | |||
| − | The value of '''n_chem_tracers''' need to be increased from 71 to '''73'''. |
||
| − | |||
| − | The '''nm_spec''' array has been edited to add in the ALICE and BOB tracers (changes in <span style="color:red">red</span>): |
||
| − | |||
| − | nm_spec(1:n_all_tracers) = (/ & |
||
| − | 'O3 ','NO ','NO3 ','NO2 ','N2O5 ', & |
||
| − | 'HO2NO2 ','HONO2 ','H2O2 ','CH4 ','CO ', & !10 |
||
| − | 'HCHO ','MeOOH ','HONO ','C2H6 ','EtOOH ', & |
||
| − | 'MeCHO ','PAN ','C3H8 ','n-PrOOH ','i-PrOOH ', & !20 |
||
| − | 'EtCHO ','Me2CO ','MeCOCH2OOH','PPAN ','MeONO2 ', & |
||
| − | 'O3_S ','C5H8 ','ISOOH ','ISON ','MACR ', & !30 |
||
| − | 'MACROOH ','MPAN ','HACET ','MGLY ','NALD ', & |
||
| − | 'HCOOH ','MeCO3H ','MeCO2H ','H2O ','ISO2 ', & !40 |
||
| − | 'Cl ','ClO ','Cl2O2 ','OClO ','Br ', & |
||
| − | 'BrO ','BrCl ','BrONO2 ','N2O ','HCl ', & !50 |
||
| − | 'HOCl ','HBr ','HOBr ','ClONO2 ','CFCl3 ', & |
||
| − | 'CF2Cl2 ','MeBr ','N ','O(3P) ','MACRO2 ', & !60 |
||
| − | 'MeCl ','CF2ClBr ','CCl4 ','<span style="color:red">ALICE</span> ','<span style="color:red">BOB</span> ', & |
||
| − | 'MeCCl3 ','CF3Br ','H2OS ','CH2Br2 ','H2 ', & !70 |
||
| − | 'DMS ','SO2 ','H2SO4 ','MSA ','DMSO ', & |
||
| − | 'NH3 ','CS2 ','COS ','H2S ','H ', & !80 |
||
| − | 'OH ','HO2 ','MeOO ','EtOO ','MeCO3 ', & |
||
| − | 'n-PrOO ','i-PrOO ','EtCO3 ','MeCOCH2OO ','MeOH ', & !90 |
||
| − | 'Monoterp ','Sec_Org ','SESQUITERP','SO3 ','AROM ', & |
||
| − | 'O(3P)_S ','O(1D)_S ','NO2 ','BrO ','HCl ', & !100 |
||
| − | 'ND_Nuc_SOL','Nuc_SOL_SU','ND_Ait_SOL','Ait_SOL_SU','Ait_SOL_BC', & |
||
| − | 'Ait_SOL_OC','ND_Acc_SOL','Acc_SOL_SU','Acc_SOL_BC','Acc_SOL_OC', & !110 |
||
| − | 'Acc_SOL_SS','Acc_SOL_DU','ND_Cor_SOL','Cor_SOL_SU','Cor_SOL_BC', & |
||
| − | 'Cor_SOL_OC','Cor_SOL_SS','Cor_SOL_DU','ND_Ait_INS','Ait_INS_BC', & !120 |
||
| − | 'Ait_INS_OC','ND_Acc_INS','Acc_INS_DU','ND_Cor_INS','Cor_INS_Du', & |
||
| − | 'Nuc_SOL_OC','Ait_SOL_SS','Nuc_SOL_OZ','Ait_SOL_OZ','Acc_SOL_OZ', & !130 |
||
| − | 'Cor_SOL_OZ','Nuc_SOL_NH','Ait_SOL_NH','Acc_SOL_NH','Cor_SOL_NH', & |
||
| − | 'Nuc_SOL_NT','Ait_SOL_NT','Acc_SOL_NT','Cor_SOL_NT','XXX ', & !140 |
||
| − | 'Anth_Prec ','Bio_Prec ','Anth_Cond ','Bio_Cond ','XXX ', & |
||
| − | 'XXX ','XXX ','XXX ','PASSIVE O3','AGE OF AIR' & !150 |
||
| − | /) |
||
| − | |||
| − | ===Changes to ukca_chem_strattrop.F90=== |
||
| − | |||
| − | The following lines have been added to the end of the '''chch_defs_strattrop_chem''' specification: |
||
| − | |||
| − | chch_t( 76,'ALICE ', 1,'TR ',' ', 0, 0, 0), & ! 76 |
||
| − | chch_t( 77,'BOB ', 1,'TR ',' ', 0, 0, 0) & ! 77 |
||
| − | |||
| − | (also, a comma has been added after the final ")" on line 75, specifying N2). The size of the <tt>chch_defs_strattrop_chem</tt> array has been increased from 75 to '''77'''. |
||
| − | |||
| − | ===Changes to ukca_constants.F90=== |
||
| − | |||
| − | The following lines have been added to the UKCA_CONSTANTS module: |
||
| − | |||
| − | ! UKCA Tutorial tracers |
||
| − | REAL, PARAMETER :: C_ALICE = 1.0000 |
||
| − | REAL, PARAMETER :: C_BOB = 1.0000 |
||
| − | |||
| − | ===Changes to ukca_cspecies.F90=== |
||
| − | |||
| − | The following lines have been added to the UKCA_CALC_CSPECIES subroutine which is held within the UKCA_CSPECIES module |
||
| − | |||
| − | ! UKCA Tutorial Tracers |
||
| − | WHERE (advt == 'ALICE ') c_species = C_ALICE |
||
| − | WHERE (advt == 'BOB ') c_species = C_BOB |
||
| − | |||
| − | ===Hand-edit to increase values of JPCTR and JPSPEC=== |
||
| − | |||
| − | A hand-edit is required to increase the values of '''JPCTR''' and '''JPSPEC'''. This should contain the following: |
||
| − | |||
| − | ed CNTLATM<<\EOF |
||
| − | /JPCTR/ |
||
| − | d |
||
| − | i |
||
| − | JPCTR = 73, |
||
| − | . |
||
| − | /JPSPEC/ |
||
| − | d |
||
| − | i |
||
| − | JPSPEC = 77, |
||
| − | . |
||
| − | w |
||
| − | q |
||
| − | EOF |
||
| − | |||
| − | An example can be found at |
||
| − | |||
| − | /home/ukca/hand_edits/VN8.2/Tutorial/Task4.2_incr_JPvals.ed |
||
| − | |||
| − | ==Output== |
||
| − | |||
| − | If you open the '''pb''' file in your '''[[UKCA & UMUI Tutorials: Things to know before you start#Archiving|archive]]''' directory you will find that it still contains the fields |
||
| − | |||
| − | 0 : 192 145 85 1 Stash code = 34001 |
||
| − | 1 : 192 145 85 1 Stash code = 34064 |
||
| − | 2 : 192 145 85 1 Stash code = 34065 |
||
| − | |||
| − | If you open 34064 and 34065 you will find that they still only contain zeros. This will change in later tasks as we add-in emissions and reactions. |
||
| − | |||
| − | '''Sample output''' from this job can be found in |
||
| − | |||
| − | /work/n02/n02/ukca/Tutorial/sample_output/Task4.2/ |
||
| − | |||
| − | on HECToR, and in |
||
| − | |||
| − | /projects/ukca/Tutorial/sample_ouput/Task4.2/ |
||
| − | |||
| − | on MONSooN. |
||
| − | |||
| − | ==Worked Solution== |
||
| − | |||
| − | There is a worked solution to this problem in the UMUI Tutorial experiment. This is job '''d''': ''Tutorial: solution to Task 4.2 - adding new chemical tracers to UKCA''. |
||
| − | |||
| − | The code changes can be viewed by using the following FCM command |
||
| − | |||
| − | fcm diff -g fcm:um_br/dev/luke/vn8.2_UKCA_Tutorial_Solns@12143 fcm:um_br/dev/luke/vn8.2_UKCA_Tutorial_Solns@12148 |
||
| − | |||
| − | This gives the following (non-graphical) output: |
||
| − | |||
| − | Index: src/atmosphere/UKCA/ukca_setd1defs.F90 |
||
| − | =================================================================== |
||
| − | --- src/atmosphere/UKCA/ukca_setd1defs.F90 (revision 12143) |
||
| − | +++ src/atmosphere/UKCA/ukca_setd1defs.F90 (revision 12148) |
||
| − | @@ -256,7 +256,7 @@ |
||
| − | (/'NO ','CH4 ','CO ','HCHO ', & |
||
| − | 'C2H6 ','C3H8 ','Me2CO ','MeCHO ', & |
||
| − | 'C5H8 ','NO_aircrft'/) |
||
| − | - n_chem_tracers = 71 ! No chem tracers |
||
| − | + n_chem_tracers = 73 ! No chem tracers |
||
| − | nr_therm = 220 ! thermal reactions |
||
| − | nr_phot = 55 ! photolytic (ATA) |
||
| − | |||
| − | @@ -409,7 +409,7 @@ |
||
| − | 'BrO ','BrCl ','BrONO2 ','N2O ','HCl ', & !50 |
||
| − | 'HOCl ','HBr ','HOBr ','ClONO2 ','CFCl3 ', & |
||
| − | 'CF2Cl2 ','MeBr ','N ','O(3P) ','MACRO2 ', & !60 |
||
| − | - 'MeCl ','CF2ClBr ','CCl4 ','CF2ClCFCl2','CHF2Cl ', & |
||
| − | + 'MeCl ','CF2ClBr ','CCl4 ','ALICE ','BOB ', & |
||
| − | 'MeCCl3 ','CF3Br ','H2OS ','CH2Br2 ','H2 ', & !70 |
||
| − | 'DMS ','SO2 ','H2SO4 ','MSA ','DMSO ', & |
||
| − | 'NH3 ','CS2 ','COS ','H2S ','H ', & !80 |
||
| − | Index: src/atmosphere/UKCA/ukca_constants.F90 |
||
| − | =================================================================== |
||
| − | --- src/atmosphere/UKCA/ukca_constants.F90 (revision 12143) |
||
| − | +++ src/atmosphere/UKCA/ukca_constants.F90 (revision 12148) |
||
| − | @@ -267,6 +267,11 @@ |
||
| − | REAL, PARAMETER :: C_ISOSVOC2 = 2.3473 ! as C5H8 |
||
| − | REAL, PARAMETER :: C_ISOSOA = 4.4874 ! 130.0 |
||
| − | |||
| − | +! UKCA Tutorial tracers |
||
| − | + REAL, PARAMETER :: C_ALICE = 1.0000 |
||
| − | + REAL, PARAMETER :: C_BOB = 1.0000 |
||
| − | + |
||
| − | + |
||
| − | ! molecular masses in g/mol of emitted species, |
||
| − | ! for budget calculations |
||
| − | |||
| − | Index: src/atmosphere/UKCA/ukca_chem_strattrop.F90 |
||
| − | =================================================================== |
||
| − | --- src/atmosphere/UKCA/ukca_chem_strattrop.F90 (revision 12143) |
||
| − | +++ src/atmosphere/UKCA/ukca_chem_strattrop.F90 (revision 12148) |
||
| − | @@ -65,7 +65,7 @@ |
||
| − | |||
| − | |||
| − | ! ATA NLA CheST Chemistry v1.2 |
||
| − | -TYPE(CHCH_T), DIMENSION( 75), PUBLIC :: chch_defs_strattrop_chem=(/ & |
||
| − | +TYPE(CHCH_T), DIMENSION( 77), PUBLIC :: chch_defs_strattrop_chem=(/ & |
||
| − | chch_t( 1,'O(3P) ', 1,'TR ','Ox ', 0, 0, 0), & ! 1 |
||
| − | chch_t( 2,'O(1D) ', 1,'SS ','Ox ', 0, 0, 0), & ! 2 |
||
| − | chch_t( 3,'O3 ', 1,'TR ','Ox ', 1, 0, 0), & ! 3 DD: 1, |
||
| − | @@ -140,7 +140,9 @@ |
||
| − | chch_t( 72,'MeOH ', 1,'TR ',' ', 1, 1, 0), & ! 72 DD:36,WD:29, |
||
| − | chch_t( 73,'CO2 ', 1,'CT ',' ', 0, 0, 0), & ! 73 |
||
| − | chch_t( 74,'O2 ', 1,'CT ',' ', 0, 0, 0), & ! 74 |
||
| − | -chch_t( 75,'N2 ', 1,'CT ',' ', 0, 0, 0) & ! 75 |
||
| − | +chch_t( 75,'N2 ', 1,'CT ',' ', 0, 0, 0), & ! 75 |
||
| − | +chch_t( 76,'ALICE ', 1,'TR ',' ', 0, 0, 0), & ! 76 |
||
| − | +chch_t( 77,'BOB ', 1,'TR ',' ', 0, 0, 0) & ! 77 |
||
| − | /) |
||
| − | |||
| − | TYPE(CHCH_T), DIMENSION( 87), PUBLIC :: chch_defs_strattrop_aer=(/ & |
||
| − | Index: src/atmosphere/UKCA/ukca_cspecies.F90 |
||
| − | =================================================================== |
||
| − | --- src/atmosphere/UKCA/ukca_cspecies.F90 (revision 12143) |
||
| − | +++ src/atmosphere/UKCA/ukca_cspecies.F90 (revision 12148) |
||
| − | @@ -270,6 +270,9 @@ |
||
| − | WHERE (advt == 'ORGNIT ') c_species = c_orgnit |
||
| − | WHERE (advt == 'PASSIVE O3') c_species = 1.0 |
||
| − | WHERE (advt == 'AGE OF AIR') c_species = 1.0 |
||
| − | +! UKCA Tutorial Tracers |
||
| − | + WHERE (advt == 'ALICE ') c_species = C_ALICE |
||
| − | + WHERE (advt == 'BOB ') c_species = C_BOB |
||
| − | |||
| − | ! non-advected tracers |
||
| − | c_na_species=0.0 |
||
| − | |||
| − | |||
| − | |||
| − | ---- |
||
| − | ''Written by [[User:Nla27 | Luke Abraham]] 2014'' |
||
