limit daytime irradiance

source/mom/tracer.F

@@ -390,7 +390,7 @@
      &           *1.2
 #  endif
             dayfrac = min( 1., -tan(tlat(i,jrow)/radian)*tan(declin))
-            dayfrac = max(1e-12, acos(max(-1., dayfrac))/pi)
+            dayfrac = max(1.e-12, acos(max(-1., dayfrac))/pi)
 #  if defined O_embm
 ! swr is in W/m2 = kg/s3, dnswr is in g/s3 = mW/m2
 ! from Andreas: swr is 24h-avg.
@@ -399,7 +399,7 @@
             swr = 200.
 #  endif
 #  if defined O_opem
-     &            /dayfrac ! with cmo, swr is mean DAYTIME irradiance
+     &            /MAX(dayfrac, 1.e-3) ! with cmo, swr is mean DAYTIME irradiance
             glb = swr
             expo_c = 0.0
             expo_p = 0.0