debugging tries for read routine

JessicaMeixner-NOAA · Apr 3, 2024 · 7929998 · 7929998
1 parent 3502b6d
commit 7929998
Showing 1 changed file with 8 additions and 1 deletion.
diff --git a/model/src/w3iopomd.F90 b/model/src/w3iopomd.F90
@@ -1224,11 +1224,18 @@ SUBROUTINE W3IOPON_READ(IOTST, IMOD, filename, ncerr)
         CALL EXTCDE ( 12 )
       END IF
 
+      !check reading: 
+      write(*,*) 'MK,MTH', MK, MTH
+      write(*,*) 'dimensions:', d_nopts_len,d_grdidlen_len, d_namelen_len,d_vsize_len,d_nspec_len
+
+
       ! Read vars with nopts as a dimension.
       ncerr = nf90_inq_varid(fh, VNAME_PTLOC, v_ptloc)
       if (ncerr .ne. 0) return
-      ncerr = nf90_get_var(fh, v_ptloc, PTLOC)
+      ncerr = nf90_get_var(fh, v_ptloc, PTLOC, start = (/ 1, 1/), &
+       count = (/ d_vsize_len, d_nopts_len /)) 
       if (ncerr .ne. 0) return
+      write(*,*) 'PTLOC', PTLOC
       ncerr = nf90_inq_varid(fh, VNAME_PTNME, v_ptnme)
       if (ncerr .ne. 0) return
       !code segfaults reading this, skipping for now to see other issues