I don’t understand if it is running through different PFTs why does it run for three times (sf_stom) ? JULES has 9 PFTs / 2 PFT (C3 and C4). None of them should make sf_stom run for three times within one timestep? Do you know how to figure out what PFT is it currently using? Thanks for your time and help. Very much appreciated.
Hi Wenyao:
Now, I can’t enter your directory
/home/users/gwy1998/JULES/r24333_pmodel/src/science/surface/
I think you need to give me the execute permission for that directory.
It should match the permissions of the other subdirectories of
/home/users/gwy1998/JULES/r24333_pmodel/src/science/
Patrick
Hi again Wenyao:
I used MOSRS to checkout a copy of your most-recently-checked-in JULES branch: r24333_pmodel
In line 880 of your file:
r24333_pmodel/src/science/surface/physiol_jls_mod.F90
there is a loop over the PFTs:
DO n = 1,npft
This loop calls in-serial (not in-parallel) the sf_stom( ) subroutine on line 1041, with the PFT index n for many of the I/O variables of that subroutine.
So, it should call that subroutine npft=9 times.
You can use several more print statements to figure out why it’s only getting 3 times (instead of 9 times) from line 880 of physiol_jls_mod.f90 to line 820 of sf_stom_jls_mod.F90 (where it has print*,'iparb',ipar(l)).
Patrick