Hi Lars,
Could you check the code in this line?
The array is allocated as allocate(id_lstate_in_pstate(dim_l)).
In the case that all variables are ndims == 2, dim_l is the number of variables. The check here uses wet_pts(3,domain_p), which is the number of wet layers. I assume it only makes sense for 3D variables. I guess a more appropriate check is id_lstate_in_pstate(dim_l).
|
if (id_lstate_in_pstate(wet_pts(3,domain_p)) > dim_state_p) then |
Cheers,
Yumeng
Hi Lars,
Could you check the code in this line?
The array is allocated as
allocate(id_lstate_in_pstate(dim_l)).In the case that all variables are
ndims == 2,dim_lis the number of variables. The check here uses wet_pts(3,domain_p), which is the number of wet layers. I assume it only makes sense for 3D variables. I guess a more appropriate check isid_lstate_in_pstate(dim_l).NEMO-PDAF/MY_SRC_pdaf/init_dim_l_pdaf.F90
Line 140 in 69616c6
Cheers,
Yumeng