geoschem / geos-chem

GEOS-Chem "Science Codebase" repository. Contains GEOS-Chem science routines, run directory generation scripts, and interface code. This repository is used as a submodule within the GCClassic and GCHP wrappers, as well as in other modeling contexts (external ESMs).
http://geos-chem.org
Other
164 stars 154 forks source link

[QUESTION] The model value of NH4 wet deposition #1518

Closed LuoBiaodll closed 1 year ago

LuoBiaodll commented 1 year ago

Ask a question about GEOS-Chem:

I simulated 1-month and got wetlossls and wetlossconv diagnostics file in GCHP 13.2.1. I found that the unit is kg/s, how can i convert it into kg/m2? I got kg/m2 by summing the wetloss vertical layers, multiplying by the simulation time and dividing by the grid area. However, I found that the conversion result is quite different from the actual. In addition, I found the Budgetwetdep diagnostic value is negative when the wetloss summing value is positive. Is that right?

Could you give me some suggestions? Thanks very much for your time and help.

yantosca commented 1 year ago

Thanks for writing @LuoBiaodll. Each diagnostic file has an AREA variable containing the surface area. You can use that to convert to kg/m2.

Also the budget diagnostic is not a true budget but just the change in concentration after each operation. Also keep in mind that re-evaporation is treated as a negative deposition.

yantosca commented 1 year ago

Tagging @Jourdan-He and @SaptSinha for reference.

LuoBiaodll commented 1 year ago

Thanks for your reply. I have checked the diagnostic file, in which i didn't find AREA varaible. I also checked the HISTORY.rc. For budget, may be it's difference between initial and final (final minus innitial)?

yantosca commented 1 year ago

@LuoBiaodll: are you using GCHP? If so then the AREA would not be in the diagnostic files (it is only for GCClassic). You can also use the area values (m^2) stored in Met_AREAM2 in the StateMet collection.

LuoBiaodll commented 1 year ago

Yes, i am using GCHP. I will turn on 'StateMet' collection. Thank you very much!

yantosca commented 1 year ago

Thanks @LuoBiaodll. I'll close out this issue now.