Merge pull request #1531 from emkemp/fix/ldt76_smap_e_opl

SMAP_E_OPL fixes and scripts
NASA-LIS · Apr 29, 2024 · f446cb8 · f446cb8
2 parents bb3ac32 + 291b358
commit f446cb8
Show file tree

Hide file tree

Showing 8 changed files with 1,410 additions and 245 deletions.
diff --git a/ldt/SMAP_E_OPL/LDT_smap_e_oplMod.F90 b/ldt/SMAP_E_OPL/LDT_smap_e_oplMod.F90
@@ -274,7 +274,7 @@ subroutine LDT_smap_e_oplRun(n)
     real                    :: teff_02(LDT_rc%lnc(n),LDT_rc%lnr(n))
     real                    :: teff_03(LDT_rc%lnc(n),LDT_rc%lnr(n))
     real                    :: SnowDepth(LDT_rc%lnc(n),LDT_rc%lnr(n))
-    real                    :: TIMEsec(LDT_rc%lnc(n),LDT_rc%lnr(n))
+    real*8                  :: TIMEsec(LDT_rc%lnc(n),LDT_rc%lnr(n))
     real                    :: UTChr(LDT_rc%lnc(n),LDT_rc%lnr(n))
     integer                 :: L1B_dir_len
     integer                 :: doy_pre, doy_curr
@@ -502,6 +502,8 @@ subroutine LDT_smap_e_oplRun(n)
              ! get UTC
              call get_UTC(n,TIMEsec,UTChr)
 
+             !write(LDT_logunit,*)'EMK: UTChr = ', UTChr
+
              ! retrieve
              ierr = LDT_create_subdirs(len_trim(SMAPeOPL%SMoutdir), &
                 trim(SMAPeOPL%SMoutdir))

diff --git a/ldt/SMAP_E_OPL/SMAPL1BTOL1C_ARFS_SUBSET.F90 b/ldt/SMAP_E_OPL/SMAPL1BTOL1C_ARFS_SUBSET.F90
@@ -15,6 +15,7 @@
 !  18 Dec 2021: Yonghwan Kwon; modified for LDT
 !  09 Feb 2023: Eric Kemp; now processes subset of fields, no output of
 !    data to separate binary files.
+!  19 Apr 2024: Eric Kemp; changed pixel times to REAL*8.
 !
 ! DESCRIPTION: RESAMPLE SMAPL1B TB TO AIR FORCE GRID
 ! INPUT : SMAP - L1B Brightness Temperature
@@ -45,12 +46,14 @@ subroutine SMAPL1BRESAMPLE_SUBSET(SMAPFILE,L1B_dir,Orbit,ARFS_TIME,rc)
   integer :: ierr
   integer :: rc
 
-  REAL*4,DIMENSION(:,:),ALLOCATABLE :: TIME_L1B, TBV_COR_L1B
+  REAL*8,DIMENSION(:,:),ALLOCATABLE :: TIME_L1B
+  REAL*4,DIMENSION(:,:),ALLOCATABLE :: TBV_COR_L1B
   REAL*4,DIMENSION(:,:),ALLOCATABLE :: LAT_L1B, LON_L1B, SCNANG_L1B
   REAL*4,DIMENSION(:),ALLOCATABLE :: ANTSCN_L1B
   INTEGER*4,DIMENSION(:,:),ALLOCATABLE :: TBVFLAG_L1B, TBHFLAG_L1B
   REAL*8,DIMENSION(:), ALLOCATABLE :: ARFS_LAT, ARFS_LON
-  REAL*4,DIMENSION(2560,1920) :: ARFS_TIME, ARFS_COR_TBV
+  REAL*8,DIMENSION(2560,1920) :: ARFS_TIME
+  REAL*4,DIMENSION(2560,1920) :: ARFS_COR_TBV
 
   REAL :: T1, T2