add new mx->mx weights

DeniseWorthen · Apr 12, 2024 · fbc751c · fbc751c
1 parent f1f9464
commit fbc751c
Showing 1 changed file with 21 additions and 0 deletions.
diff --git a/sorc/cpld_gridgen.fd/gen_fixgrid.F90 b/sorc/cpld_gridgen.fd/gen_fixgrid.F90
@@ -524,7 +524,17 @@ program gen_fixgrid
         fwgt = trim(dirout)//'/'//'tripole.mx025.Ct.to.mx'//trim(res)//'.Ct.neareststod.nc'
         logmsg = 'creating weight file '//trim(fwgt)
         print '(a)',trim(logmsg)
+        call ESMF_RegridWeightGen(srcFile=trim(fsrc),dstFile=trim(fdst), &
+             weightFile=trim(fwgt), regridmethod=method,                 &
+             ignoreDegenerate=.true., unmappedaction=ESMF_UNMAPPEDACTION_IGNORE, rc=rc)
+        if (ESMF_LogFoundError(rcToCheck=rc, msg=ESMF_LOGERR_PASSTHRU, &
+             line=__LINE__, file=__FILE__)) call ESMF_Finalize(endflag=ESMF_END_ABORT)
 
+        method=ESMF_REGRIDMETHOD_BILINEAR
+        fdst = trim(dirout)//'/'//'Ct.mx'//trim(res)//'_SCRIP.nc'
+        fwgt = trim(dirout)//'/'//'tripole.mx025.Ct.to.mx'//trim(res)//'.Ct.bilinear.nc'
+        logmsg = 'creating weight file '//trim(fwgt)
+        print '(a)',trim(logmsg)
         call ESMF_RegridWeightGen(srcFile=trim(fsrc),dstFile=trim(fdst), &
              weightFile=trim(fwgt), regridmethod=method,                 &
              ignoreDegenerate=.true., unmappedaction=ESMF_UNMAPPEDACTION_IGNORE, rc=rc)
@@ -537,6 +547,17 @@ program gen_fixgrid
      end if
   end if
 
+  ! tripole Ct->tripole Bu for CICE are only for CICE IC creation
+  fsrc = trim(dirout)//'/'//'Ct.mx'//trim(res)//'_SCRIP.nc'
+  fdst = trim(dirout)//'/'//'Bu.mx'//trim(res)//'_SCRIP.nc'
+  fwgt = trim(dirout)//'/'//'tripole.mx'//trim(res)//'.Ct.to.Bu.bilinear.nc'
+  logmsg = 'creating weight file '//trim(fwgt)
+  print '(a)',trim(logmsg)
+  call ESMF_RegridWeightGen(srcFile=trim(fsrc),dstFile=trim(fdst), &
+       weightFile=trim(fwgt), regridmethod=method,                 &
+       ignoreDegenerate=.true., unmappedaction=ESMF_UNMAPPEDACTION_IGNORE, rc=rc)
+  if (ESMF_LogFoundError(rcToCheck=rc, msg=ESMF_LOGERR_PASSTHRU, &
+       line=__LINE__, file=__FILE__)) call ESMF_Finalize(endflag=ESMF_END_ABORT)
   !---------------------------------------------------------------------
   !
   !---------------------------------------------------------------------