-
Notifications
You must be signed in to change notification settings - Fork 2
Integration of tropical PEATCLSM in GEOSldas (CLSM and CatchmentCN) - Round 1 #1
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.
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
base: develop
Are you sure you want to change the base?
Conversation
…s the first changes for CN - catchmentCN.F90 is set-up so that GEOSldas compiles
|
||
! MB: FOXY only called here, ZBAR unchanged | ||
CALL OXYFAC ( & | ||
NCH, ZBAR, POROS, & |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
ZBAR is input to this function but not calculated yet, see line above that you missed to add.
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Cannot add it to the line above. So I write it here. Add at line 1842:
!ZBAR defined here positive below ground and in meter
ZBAR(CHNO)=SQRT(1.e-20+CATDEF(CHNO)/BF1(CHNO))-BF2(CHNO)
compare with 3184 of https://github.com/mbechtold/PEATCLSM_T/blob/main/catchment.F90
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
suggested bug fix for OXYFAC error message
! comments) | ||
ZBAR1 = catch_calc_zbar( BF1(N), BF2(N), CATDEF(N) ) | ||
SYSOIL = (2.*bf1(N)*amin1(amax1(zbar1,0.),PEATCLSM_ZBARMAX_4_SYSOIL) + 2.*bf1(N)*bf2(N))/1000. | ||
SYSOIL = (2.*bf1(N)*amin1(amax1(zbar1,0.),0.45) + 2.*bf1(N)*bf2(N))/1000. |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
This should become and "if statement" where you have the three porosity cases, I assume.
There are 4 files 'added' and 4 files 'adapted'.
The 4 added files: 3 of them _original are the original files and 1 of them _SA is not finished yet with my adaptations.
The 4 adapted files: catchment.F90, lsm_routines.F90 and catch_constants.f90 include my final adaptations and catchmentCN.F90 is adapted so that is allows GEOSldas to compile with my adaptation in the other 3 files but is still quite similar to the original one.