Comments (8)
Thanks for writing @ast222340. We do not have much experience with Matlab, so we would not be able to assist you with that. However, GCPy (our Python package) has the functionality to read and plot data on on GCHP cubed-sphere/stretched grids. We would recommend that you look into using GCPy for your analysis. See: https://gcpy.readthedocs.io.
I think that the grid in the first figure would be more proper for the entire country of India (even though you have some border countries included).
from gchp.
Thank you @yantosca.I was late replying because you tagged them.I have started the STRETCHED GRID simulation of the first figure as you said.I am using "GEOSChem.Restart.fullchem.20190701_0000z.c180.nc4" file for STRETCHED GRID which is available in ExtData as GEOSCHEM_RESTARTS.My reason is that I use C180 because comparing with 05x0625.I am facing many errors.
Integer*4 Resource Parameter: GCHPchem_REFERENCE_TIME:500
Character Resource Parameter: GCHPchem_INTERNAL_CHECKPOINT_FILE:Restarts/gcchem_internal_checkpoint
Character Resource Parameter: GCHPchem_INTERNAL_RESTART_FILE:gchp_restart.nc4
Character Resource Parameter: MAPL_ENABLE_BOOTSTRAP:YES
pe=00011 FAIL at line=06123 MAPL_Generic.F90
pe=00011 FAIL at line=01609 MAPL_Generic.F90
pe=00015 FAIL at line=06123 MAPL_Generic.F90
I just showed a variable in that restart file.
SPC_N2O
Size: 180x1080x72x1
Dimensions: lon,lat,lev,time
Datatype: single
I don't understand why the model is not running.Please let me know where I am doing wrong.
GCHP.log.txt
from gchp.
@ast222340, please read the stretched grid section of GCHP ReadTheDocs. It states that there is a stretched grid restart file available for a test case. The test case is not at c180 so it sounds like you are trying to use a regular grid c180 restart file for your stretched grid run. This will fail. You need to make a restart file that matches your grid resolution and stretch grid parameters.
Please also make sure that if you are using the standard simulation you update setCommonRunSettings.sh
to allow missing variables in the restart file. You can find the setting for that by scrolling through the configuration section of the file and reading the comments. Ideally your restart file would be generated to match the version you are running with. You can do this by creating your stretched grid restart file from one of the restarts that is included in the Restarts
subdirectory in the run directory you create from your source code.
Lastly, the latest ReadTheDocs has an updated debugging page. Please see https://gchp.readthedocs.io/en/stable/user-guide/debugging.html.
from gchp.
Thank you @lizziel.As I read GCHP ReadTheDocs,I realized that need to convert the restart file mentioned above to stretched grid using regridding method(like GEOSChem.Restart.20190701_0000z.c180.s3_21N_82E.nc').Then transfer this new file to the Restarts subdirectory in the run directory.Did I say it right?please correct me if i am wrong.
When STRETCH_GRID=OFF,then model started running fine.I ran the first 1 months for testing.Then a file was created on the Restarts subdirectory as name "gcchem_internal_checkpoint".I renamed the checkpoint file (GEOSChem.Restart.fullchem.20190801_0000z.c180.nc4) to long run for 1 year.This process is right.If this is wrong,please tell me otherwise my analysis will be wrong.
This restart file is provided for testing purposes.Can't i do the simulation without taking restart file just for asking?if possible,where to change.it is required for long simulations because if the model stops midway.
setCommonRunSettings.sh.txt
from gchp.
Hi @ast222340, that is correct that you need to create a restart file for your stretched grid simulation. When you transfer it to the Restarts subdirectory make sure that you set the symbolic link to it before running. See the GCHP ReadTheDocs section called "Configure run directory" within the stretched grid page for how to do this.
That's great that you have the model working for C180 with stretched grid off. You can avoid having to manually rename the gcchem_internal_checkpoint
file by using a run script from subdirectory runScriptSamples
.
I renamed the checkpoint file (GEOSChem.Restart.fullchem.20190801_0000z.c180.nc4) to long run for 1 year.
Could you clarify what you mean?
This restart file is provided for testing purposes.Can't i do the simulation without taking restart file just for asking?if possible,where to change.it is required for long simulations because if the model stops midway.
I'm not sure I understand your question. Can you clarify?
from gchp.
Thank You @lizziel.
I renamed the checkpoint file (GEOSChem.Restart.fullchem.20190801_0000z.c180.nc4) to long run for 1 year.
---Example,I have a Restart file "GEOSChem.Restart.fullchem.20190801_0000z.c180.nc4".First i ran it for a month for testing purposes.A gcchem_internal_checkpoint is created in the Restart subdirectory at the end of the run.I was saying that using this check file I can run for one year.
This restart file is provided for testing purposes.Can't i do the simulation without taking restart file just for asking?if possible,where to change.it is required for long simulations because if the model stops midway
GEOS Chem classic model does not required a restart file without test case,it has an option to create it automatically via "HISTORY.rc" file.
==============================================================================
%%%%% THE Restart COLLECTION %%%%%
GEOS-Chem restart file fields
Available for all simulations
==============================================================================
Restart.filename: './Restarts/GEOSChem.Restart.%y4%m2%d2_%h2%n2z.nc4',
Restart.frequency: 00000001 000000
Restart.duration: 00000001 000000
Restart.mode: 'instantaneous'
Restart.fields: 'SpeciesRst_?ALL? ',
If this model had such an option then we would not need to provide any restart file initially.
I used "MATLAB" for all of my research work.I would like to know how to plot variables in Matlab, such as SpeciesConcVV_CO.'Cubed-sphere' grid plots are something I have never done.If you could explain the process to me, let me know what method to use in this situation, or provide some sample Matlab code.
from gchp.
Hi @ast222340, the code in HISTORY.rc you are referring to is the output restart file for GEOS-Chem Classic. You need an initial restart file for GEOS-Chem Classic. Like GCHP, a restart file is provided in the run directory.
For both GEOS-Chem Classic and GCHP we recommend using the provided restart file to do a "spin-up", meaning a simulation run prior to your science run start date. This better captures the state of the atmosphere at the date you wish to start simulation the atmosphere. We also recommend taking seasonality into account. For example, If you have a July 2019 restart file that comes with the run directory and you want to start your simulation on January 1st 2018, then you would rename your GCHP restart file to July 1st 2017 and run for 6 months. The output restart file would be the restart file to use for your science run. This assumes a 6 month spin-up is okay for the research you are doing. If you are looking at the stratosphere then you may want to spin up the model for longer.
The GEOS-Chem Support Team does not advise on use of MATLAB for analyzing GEOS-Chem data. We provide the python package GCPy which includes examples for plotting. See https://github.com/geoschem/gcpy.
from gchp.
I am closing this issue because the original questions about stretched grid are resolved.
from gchp.
Related Issues (20)
- Feature request: Update the docs/requirements.yml to avoid security issues w/ old python versions (used for ReadTheDocs) HOT 1
- Problem about compilation of GCHP HOT 11
- Wrong units for SLP and TROPP in preprocessed GEOS-IT ExtData.rc HOT 2
- Constant value for all grid boxes in Passive Tracer monthly mean diagnostic HOT 1
- Single column/grid subset mode? HOT 6
- Access issue with restart files of GCHP 14.3.2 HOT 10
- Stretched Grid Runs Failing with "Error calling DO_WETDEP" HOT 7
- Create Run Directory Problem with 14.2.2 and 14.2.3 container HOT 4
- Using gcchem_internal_checkpoint for grid-stretching simulation but fails with 'Factories not equal' error. HOT 15
- GCHP simulation stopped after 7 month simulation (total set time for 1 yr) HOT 9
- Wrong surface type field when using raw meteorology file HOT 7
- [SUBMODULE UPDATES] MAPL fix for restarting with stretched grid checkpoints HOT 1
- GCHP carbon simulations with CH4 take much longer than simulations without CH4 HOT 14
- Using a new inventory in ExtData.rc HOT 8
- GCHP 14.3.1 out of memory when writing checkpoint files HOT 12
- Error with CMIP6 file HOT 26
- ERROR during 'Initialize' stage of the gridded component 'EXTDATA' HOT 6
- Transport Budget Diagnostics at surface level for GCHP14.3.1 C180 simulation are all zero HOT 7
- error running carbon simulation HOT 16
- Segmentation fault - invalid memory reference in mapl_capmod_MOD_run_model HOT 2
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from gchp.