Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Jump to bottom

[PULL REQUEST] More stable stratospheric H2O boundary condition #689

Merged
merged 1 commit into from Apr 21, 2021
Merged

[PULL REQUEST] More stable stratospheric H2O boundary condition #689

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
1 commit
Select commit
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Jump to
Jump to file
Failed to load files.
Diff view
Diff view
15 changes: 10 additions & 5 deletions GeosCore/input_mod.F90
View file
Open in desktop
Expand Up @@ -1961,6 +1961,14 @@ SUBROUTINE READ_EMISSIONS_MENU( Input_Opt, RC )
ENDIF
READ( SUBSTRS(1:N), * ) Input_Opt%LSETH2O

! Use static strat H2O boundary condition?
CALL SPLIT_ONE_LINE( SUBSTRS, N, 1, 'LSTATICH2OBC', RC )
IF ( RC /= GC_SUCCESS ) THEN
CALL GC_Error( ErrMsg, RC, ThisLoc )
RETURN
ENDIF
READ( SUBSTRS(1:N), * ) Input_Opt%LStaticH2OBC

! Separator line
CALL SPLIT_ONE_LINE( SUBSTRS, N, 1, 'separator 4', RC )
IF ( RC /= GC_SUCCESS ) THEN
Expand All @@ -1972,11 +1980,6 @@ SUBROUTINE READ_EMISSIONS_MENU( Input_Opt, RC )
! Error check logical flags
!=================================================================

! Turn off full-chem only switches
IF ( .not. Input_Opt%ITS_A_FULLCHEM_SIM ) THEN
Input_Opt%LSETH2O = .FALSE.
ENDIF

! Return success
RC = GC_SUCCESS

Expand All @@ -1994,6 +1997,8 @@ SUBROUTINE READ_EMISSIONS_MENU( Input_Opt, RC )
Input_Opt%LCH4EMIS
WRITE( 6, 100 ) 'Set initial strat H2O? : ', &
Input_Opt%LSETH2O
WRITE( 6, 100 ) 'Use robust strat H2O BC? : ', &
Input_Opt%LStaticH2OBC
ENDIF

! FORMAT statements
Expand Down
8 changes: 6 additions & 2 deletions GeosCore/ucx_mod.F90
View file
Open in desktop
Expand Up @@ -3678,8 +3678,12 @@ SUBROUTINE SET_H2O_TRAC( SETSTRAT, Input_Opt, State_Chm, State_Grid, &
IF ( ABS(State_Grid%YMid(I,J)) <= 30 ) THEN
! Level with minimum temperature,
! use MASK to screen for >10 hPa
LEVCPT = MINLOC( State_Met%T(I,J,:), DIM=1, &
MASK=(State_Met%PMID(I,J,:) >= 10) )
If (Input_Opt%LStaticH2OBC) Then
LEVCPT = MINLOC(ABS(State_Met%PMID(I,J,:) - 70), DIM=1)
Else
LEVCPT = MINLOC( State_Met%T(I,J,:), DIM=1, &
MASK=(State_Met%PMID(I,J,:) >= 10) )
End If
ELSE
LEVCPT = -1
ENDIF
Expand Down
2 changes: 2 additions & 0 deletions Headers/input_opt_mod.F90
View file
Open in desktop
Expand Up @@ -147,6 +147,7 @@ MODULE Input_Opt_Mod
LOGICAL :: LCH4EMIS
LOGICAL :: LCH4SBC
LOGICAL :: LSETH2O
LOGICAL :: LStaticH2OBC
LOGICAL :: LHCodedOrgHal
LOGICAL :: LCMIP6OrgHal
LOGICAL :: DoLightNOx ! Shadow for LightNOX extension
Expand Down Expand Up @@ -643,6 +644,7 @@ SUBROUTINE Set_Input_Opt( am_I_Root, Input_Opt, RC )
Input_Opt%LCH4EMIS = .FALSE.
Input_Opt%LCH4SBC = .FALSE.
Input_Opt%LSETH2O = .FALSE.
Input_Opt%LStaticH2OBC = .FALSE.
Input_Opt%LHCodedOrgHal = .FALSE.
Input_Opt%LCMIP6OrgHal = .FALSE.
Input_Opt%DoLightNOx = .FALSE.
Expand Down
1 change: 1 addition & 0 deletions run/GCHP/input.geos.templates/input.geos.TransportTracers
View file
Open in desktop
Expand Up @@ -53,6 +53,7 @@ HEMCO Input file : HEMCO_Config.rc
Switches for UCX :---
=> Use CH4 emissions? : F
=> Set init. strat. H2O: F
=> Use static H2O BC : T
------------------------+------------------------------------------------------
%%% AEROSOL MENU %%% :
Online SULFATE AEROSOLS : F
Expand Down
1 change: 1 addition & 0 deletions run/GCHP/input.geos.templates/input.geos.fullchem
View file
Open in desktop
Expand Up @@ -235,6 +235,7 @@ HEMCO Input file : HEMCO_Config.rc
Switches for UCX :---
=> Use CH4 emissions? : F
=> Set init. strat. H2O: T
=> Use static H2O BC : T
------------------------+------------------------------------------------------
%%% AEROSOL MENU %%% :
Online SULFATE AEROSOLS : T
Expand Down