Modified calculation of latin1 for grid template 1 in wind rotation

LarissaReames-NOAA · Feb 26, 2024 · 804ea27 · 804ea27
1 parent 492eed9
commit 804ea27
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 2 deletions.
diff --git a/sorc/chgres_cube.fd/atm_input_data.F90 b/sorc/chgres_cube.fd/atm_input_data.F90
@@ -3020,7 +3020,7 @@ subroutine read_winds(u,v,localpet,octet_23,rlevs,lugb,pdt_num)
 
    elseif (gfld%igdtnum == 1) then ! grid definition template number - non-E stagger rotated lat/lon grid
 
-     latin1 = real(float(gfld%igdtmpl(20))/1.0E6, kind=esmf_kind_r4)
+     latin1 = real(float(gfld%igdtmpl(20))/1.0E6, kind=esmf_kind_r4) + 90.0_esmf_kind_r4
      lov = real(float(gfld%igdtmpl(21))/1.0E6, kind=esmf_kind_r4)
 
      print*, "- CALL CALCALPHA_ROTLATLON with center lat,lon = ",latin1,lov

diff --git a/sorc/chgres_cube.fd/model_grid.F90 b/sorc/chgres_cube.fd/model_grid.F90
@@ -1557,9 +1557,9 @@ subroutine gdt_to_gds(igdtnum, igdstmpl, igdtlen, kgds, ni, nj, res)
      kgds(7)=nint(float(igdstmpl(20))/float(iscale)*1000.)       ! octs 18-20,
                                                                  ! Lat of cent of rotation
      kgds(8)=nint(float(igdstmpl(21))/float(iscale)*1000.)       ! octs 21-23,
-     print*, "INPUT LAT, LON CENTER ", kgds(7), kgds(8) 
                                                                  ! Lon of cent of rotation
      kgds(7) = kgds(7) + 90000.0
+     print*, "INPUT LAT, LON CENTER ", kgds(7), kgds(8)
 
      DPR = 180.0/3.1415926
      CLATR=COS((float(kgds(4))/1000.0)/DPR)