Bug fix

GEOS-ESM · Dec 20, 2023 · bfb1d72 · bfb1d72
1 parent 00903b5
commit bfb1d72
Show file tree

Hide file tree

Showing 2 changed files with 9 additions and 5 deletions.
diff --git a/GEOSogcm_GridComp/GEOS_OceanBioGeoChemGridComp/ta_iter_SWS.F90 b/GEOSogcm_GridComp/GEOS_OceanBioGeoChemGridComp/ta_iter_SWS.F90
@@ -71,9 +71,9 @@ subroutine ta_iter_SWS(x,fn,df)
 
 !	df = dfn/dx
 
-       if (a2 == 0.) then
-          a2 = tiny(a2)
-       end if
+!       if (a2 == 0.) then
+!          a2 = tiny(a2)
+!       end if
 
        df = ((k1*dic*b) - k1*x*dic*db)/b2 -                        &
              2.0*dic*k12*db/b2 -                                   &

diff --git a/GEOSogcm_GridComp/GEOS_OradBioGridComp/setlte.F90 b/GEOSogcm_GridComp/GEOS_OradBioGridComp/setlte.F90
@@ -1,5 +1,5 @@
-      subroutine setlte(lam,aw,bw,ac,bc,bpic,excdom,exdet,WtoQ,  &
-       wfac)
+      subroutine setlte(rad,pi2,lam,aw,bw,ac,bc,bpic,excdom,exdet,    &
+                        WtoQ,wfac)
 
 !  Sets parameters for ocean irradiance.
 
@@ -15,6 +15,10 @@ subroutine setlte(lam,aw,bw,ac,bc,bpic,excdom,exdet,WtoQ,  &
       data a0,a1,a2,a3 /0.9976,0.2194,5.554E-2,6.7E-3/
       data b0,b1,b2,b3 /5.026,-0.01138,9.552E-6,-2.698E-9/
 
+!  Degrees to radians conversion
+      pi = dacos(-1.0D0)
+      pi2 = pi*2.0
+      rad = 180.0D0/pi
 !   Obtain Light data
       call lidata(lam,aw,bw,ac,bc,bpic)
 !   Quanta conversion