This is the minimum required update to use 2020.01 FMS & FMScoupler

- Adding domain position arguments to diag_axis_init calls
dpsarmie · Jun 24, 2020 · a7ab529 · a7ab529
1 parent ac231db
commit a7ab529
Showing 1 changed file with 5 additions and 5 deletions.
diff --git a/src/framework/MOM_diag_mediator.F90 b/src/framework/MOM_diag_mediator.F90
@@ -33,7 +33,7 @@ module MOM_diag_mediator
 use diag_axis_mod, only : get_diag_axis_name
 use diag_data_mod, only : null_axis_id
 use diag_manager_mod, only : diag_manager_init, diag_manager_end
-use diag_manager_mod, only : send_data, diag_axis_init, diag_field_add_attribute
+use diag_manager_mod, only : send_data, diag_axis_init, EAST, NORTH, diag_field_add_attribute
 ! The following module is needed for PGI since the following line does not compile with PGI 6.5.0
 ! was: use diag_manager_mod, only : register_diag_field_fms=>register_diag_field
 use MOM_diag_manager_wrapper, only : register_diag_field_fms
@@ -365,14 +365,14 @@ subroutine set_axes_info(G, GV, US, param_file, diag_cs, set_vertical)
   ! Horizontal axes for the native grids
   if (G%symmetric) then
     id_xq = diag_axis_init('xq', G%gridLonB(G%isgB:G%iegB), G%x_axis_units, 'x', &
-              'q point nominal longitude', Domain2=G%Domain%mpp_domain)
+              'q point nominal longitude', Domain2=G%Domain%mpp_domain, domain_position=EAST)
     id_yq = diag_axis_init('yq', G%gridLatB(G%jsgB:G%jegB), G%y_axis_units, 'y', &
-              'q point nominal latitude', Domain2=G%Domain%mpp_domain)
+              'q point nominal latitude', Domain2=G%Domain%mpp_domain, domain_position=NORTH)
   else
     id_xq = diag_axis_init('xq', G%gridLonB(G%isg:G%ieg), G%x_axis_units, 'x', &
-              'q point nominal longitude', Domain2=G%Domain%mpp_domain)
+              'q point nominal longitude', Domain2=G%Domain%mpp_domain, domain_position=EAST)
     id_yq = diag_axis_init('yq', G%gridLatB(G%jsg:G%jeg), G%y_axis_units, 'y', &
-              'q point nominal latitude', Domain2=G%Domain%mpp_domain)
+              'q point nominal latitude', Domain2=G%Domain%mpp_domain, domain_position=NORTH)
   endif
   id_xh = diag_axis_init('xh', G%gridLonT(G%isg:G%ieg), G%x_axis_units, 'x', &
               'h point nominal longitude', Domain2=G%Domain%mpp_domain)