Difference between revisions of "CH4 simulation"

From Geos-chem
Jump to: navigation, search
(Error reading GWETTOP fields)
Line 41: Line 41:
 
Melissa Payer discovered a file read error in the current CH4 simulation.  We will fix this issue in [[GEOS-Chem v9-01-03]].
 
Melissa Payer discovered a file read error in the current CH4 simulation.  We will fix this issue in [[GEOS-Chem v9-01-03]].
  
'''''[mailto mpayer@seas.harvard.edu Melissa Payer] wrote:'''''
+
'''''[mailto:mpayer@seas.harvard.edu Melissa Payer] wrote:'''''
  
 
:I got the offline CH4 simulation rice emissions to work. It turns out that XTAU is being reset [to a the wrong value] when GWETTOP is read.  To fix this, I added the following code in RICE_EMIS prior to setting the soil wetness FILENAME:
 
:I got the offline CH4 simulation rice emissions to work. It turns out that XTAU is being reset [to a the wrong value] when GWETTOP is read.  To fix this, I added the following code in RICE_EMIS prior to setting the soil wetness FILENAME:
  
          ! Determine year for soil wetness files (mpayer, 1/31/12)
+
      ! Determine year for soil wetness files (mpayer, 1/31/12)
          IF ( GET_YEAR() < 2004 ) THEN
+
      IF ( GET_YEAR() < 2004 ) THEN
            CYEAR  = '2004'
+
        CYEAR  = '2004'
            CYEAR2 = 2004
+
        CYEAR2 = 2004
          ELSEIF (GET_YEAR() > 2009 ) THEN
+
      ELSEIF (GET_YEAR() > 2009 ) THEN
            CYEAR  = '2009'
+
        CYEAR  = '2009'
            CYEAR2 = 2009
+
        CYEAR2 = 2009
          ELSE
+
      ELSE
            WRITE( CYEAR, '(i4)' ) GET_YEAR()
+
        WRITE( CYEAR, '(i4)' ) GET_YEAR()
            CYEAR2 = GET_YEAR()
+
        CYEAR2 = GET_YEAR()
          ENDIF
+
      ENDIF
  
 
:I also changed:
 
:I also changed:
  
          XTAU = GET_TAU0( 1, 1, GET_YEAR() )
+
      XTAU = GET_TAU0( 1, 1, GET_YEAR() )
  
 
:to:
 
:to:
  
          XTAU = GET_TAU0( 1, 1, CYEAR2 )
+
      XTAU = GET_TAU0( 1, 1, CYEAR2 )
  
 
:For both annual and monthly mean GEWTTOP. That seems to work and it reads the file for the proper simulation year.  
 
:For both annual and monthly mean GEWTTOP. That seems to work and it reads the file for the proper simulation year.  
  
 
--[[User:Bmy|Bob Y.]] 10:29, 1 February 2012 (EST)
 
--[[User:Bmy|Bob Y.]] 10:29, 1 February 2012 (EST)

Revision as of 15:30, 1 February 2012

NOTE: Page under construction

This page contains information about the methane (CH4) simulation in GEOS-Chem.

Overview

The CH4 simulation in GEOS-Chem ...

Authors and collaborators

CH4 simulation user groups

User Group Personnel Projects
Harvard University Christopher Pickett-Heaps
Kevin Wecht
Methane sources in the Arctic
U. Wollongong Chris Chan Miller Methane sources in tropical Australia

Validation

References

Known issues

Error reading GWETTOP fields

Melissa Payer discovered a file read error in the current CH4 simulation. We will fix this issue in GEOS-Chem v9-01-03.

Melissa Payer wrote:

I got the offline CH4 simulation rice emissions to work. It turns out that XTAU is being reset [to a the wrong value] when GWETTOP is read. To fix this, I added the following code in RICE_EMIS prior to setting the soil wetness FILENAME:
     ! Determine year for soil wetness files (mpayer, 1/31/12)
     IF ( GET_YEAR() < 2004 ) THEN
        CYEAR  = '2004'
        CYEAR2 = 2004
     ELSEIF (GET_YEAR() > 2009 ) THEN
        CYEAR  = '2009'
        CYEAR2 = 2009
     ELSE
        WRITE( CYEAR, '(i4)' ) GET_YEAR()
        CYEAR2 = GET_YEAR()
     ENDIF
I also changed:
     XTAU = GET_TAU0( 1, 1, GET_YEAR() )
to:
     XTAU = GET_TAU0( 1, 1, CYEAR2 )
For both annual and monthly mean GEWTTOP. That seems to work and it reads the file for the proper simulation year.

--Bob Y. 10:29, 1 February 2012 (EST)