Code: Select all
ic=0
DO i=1,NT(ng)
IF (LtracerCLM(i,ng)) THEN
ic=ic+1
CALL get_3dfld (ng, iNLM, idTclm(i), CLM(ng)%ncid, &
& 1, CLM(ng), update(1), &
& LBi, UBi, LBj, UBj, 1, N(ng), 2, 1, &
& GRID(ng) % rmask, &
& CLIMA(ng) % tclmG(:,:,:,:,ic))
IF (exit_flag.ne.NoError) RETURN
END IF
END DO
when ifield=idTclm(3)=0;
It seems idTclm(3) is not initialized and something like idTclm(iNO3_) should be add to varinfo.dat and fennel_var.h.