Switch to method used by RWG to output weight matrix.

Author: oehmke

src/Infrastructure/Field/interface/ESMF_Field_C.F90

@@ -1356,10 +1356,21 @@ subroutine f_esmf_regridstorefile(srcField, dstField, fileName, &
     endif
 
     if (filemode_local == ESMF_FILEMODE_BASIC) then
+#ifndef OLDWAY
+       call ESMF_OutputSimpleWeightFile (fileName, localFactorList, localFactorIndexList, &
+            title="Weight file created by the ESMPy Regrid Class", &
+            method=regridmethod, &
+            largeFileFlag=l_largeFileFlag, &
+            ! TODO: Support this  netcdf4FileFlag, &
+            rc=localrc)
+       if (ESMF_LogFoundError(localrc, ESMF_ERR_PASSTHRU, &
+            ESMF_CONTEXT, rcToReturn=rc)) return       
+#else       
       call ESMF_SparseMatrixWrite(localFactorList, localFactorIndexList, &
                                   fileName, rc=localrc)
       if (ESMF_LogFoundError(localrc, ESMF_ERR_PASSTHRU, &
-          ESMF_CONTEXT, rcToReturn=rc)) return
+           ESMF_CONTEXT, rcToReturn=rc)) return
+#endif      
     elseif (filemode_local == ESMF_FILEMODE_WITHAUX) then
       ! query field for geom type
       call ESMF_FieldGet(srcField, geomType=srcgt, typekind=srctk, rc=localrc)