fix corner lat/lon issue, use getLatBu/getLonBu instead of ulon/vlat

jiandewang · Dec 11, 2017 · 7872ef1 · 7872ef1
1 parent 4253f8c
commit 7872ef1
Showing 1 changed file with 6 additions and 1 deletion.
diff --git a/mom_cap.F90 b/mom_cap.F90
@@ -1154,7 +1154,12 @@ subroutine InitializeRealize(gcomp, importState, exportState, clock, rc)
     enddo
     enddo
 
+#ifdef MOM5_CAP
     call ocean_model_data_get(Ocean_state, Ocean_sfc, 'ulon', ofld, isc, jsc)
+#endif
+#ifdef MOM6_CAP
+    call ocean_model_data_get(Ocean_state, Ocean_sfc, 'geoLonBu', ofld, isc, jsc)
+#endif
     write(tmpstr,*) subname//' ofld xu = ',minval(ofld),maxval(ofld)
     call ESMF_LogWrite(trim(tmpstr), ESMF_LOGMSG_INFO, rc=dbrc)
     call mpp_global_field(Ocean_sfc%domain, ofld, gfld)
@@ -1190,7 +1195,7 @@ subroutine InitializeRealize(gcomp, importState, exportState, clock, rc)
 #endif
 
 #ifdef MOM6_CAP
-    call ocean_model_data_get(Ocean_state, Ocean_sfc, 'vlat', ofld, isc, jsc)
+     call ocean_model_data_get(Ocean_state, Ocean_sfc, 'geoLatBu', ofld, isc, jsc)     
 #endif
 
     write(tmpstr,*) subname//' ofld yu = ',minval(ofld),maxval(ofld)