Skip to content

Commit

Permalink
debug mode3
Browse files Browse the repository at this point in the history
  • Loading branch information
AnningCheng-NOAA committed Nov 26, 2024
1 parent b33ccf9 commit 21c0af0
Showing 1 changed file with 4 additions and 4 deletions.
8 changes: 4 additions & 4 deletions physics/MP/Morrison_Gettelman/aerinterp.F90
Original file line number Diff line number Diff line change
Expand Up @@ -225,7 +225,7 @@ SUBROUTINE read_aerdataf_dl ( me, master, iflip, idate, FHOUR, errmsg, errflg)
write(myr,'(i4.4)') jdat(1)
write(mn,'(i2.2)') jdat(2)
write(dy,'(i2.2)') jdat(3)
fname=trim("merra2_"//myr//mn//dy//".nc")
fname="merra2_"//myr//mn//dy//".nc"
! rjday is the minutes in a day
rjday = jdat(5)*60+jdat(6)+jdat(7)/60.
!
Expand All @@ -250,7 +250,7 @@ SUBROUTINE read_aerdataf_dl ( me, master, iflip, idate, FHOUR, errmsg, errflg)
n2sv=1
t2sv=1440.
call fdnx_fname (jdat(1), jdat(2),jdat(3),fname)
call read_netfaer_dl(trim(fname), n2sv, iflip, 2, errmsg, errflg)
call read_netfaer_dl(fname, n2sv, iflip, 2, errmsg, errflg)
end if
END SUBROUTINE read_aerdataf_dl
!**********************************************************************
Expand Down Expand Up @@ -301,7 +301,6 @@ SUBROUTINE aerinterpol_dl( me,master,nthrds,npts,IDATE,FHOUR,iflip, jindx1,jindx
call w3doxdat(jdat,jdow,jdoy,jday)
! rjday is the minutes in a day
rjday = jdat(5)*60+jdat(6)+jdat(7)/60.
write(*,*)"BBBBB",jdat(2),jdat(3), jdat(5)
if(rjday >= t2sv .or. jdat(3).ne.n1sv) then !!need to either to read in a record or open a new file
#ifdef DEBUG
if (me == master) write(*,*)"read in a new MERRA2 record", n2sv+1
Expand All @@ -326,6 +325,7 @@ SUBROUTINE aerinterpol_dl( me,master,nthrds,npts,IDATE,FHOUR,iflip, jindx1,jindx
write(mn,'(i2.2)') jdat(2)
write(dy,'(i2.2)') jdat(3)
fname="merra2_"//myr//mn//dy//".nc"
write(*,*)"BBBB", fname
call read_netfaer_dl(fname,n2sv, iflip, 2, errmsg, errflg)
write(*,*)"AAAAAND",fname, t1sv, rjday, t2sv
else if (rjday >= t2sv) then
Expand All @@ -346,7 +346,7 @@ SUBROUTINE aerinterpol_dl( me,master,nthrds,npts,IDATE,FHOUR,iflip, jindx1,jindx
n2sv=1
t2sv=1440.
call fdnx_fname (jdat(1), jdat(2),jdat(3),fname)
call read_netfaer_dl(trim(fname), n2sv, iflip, 2, errmsg, errflg)
call read_netfaer_dl(fname, n2sv, iflip, 2, errmsg, errflg)
write(*,*)"AAAAANF",fname, t1sv, rjday, t2sv
end if
end if
Expand Down

0 comments on commit 21c0af0

Please sign in to comment.