Changeset 204 for trunk/src/calc_cloud_vlat.pro

Timestamp:

03/26/14 11:58:56 (10 years ago)

Author:

pinsard

Message:

fix thanks to coding rules; typo

File:

• trunk/src/calc_cloud_vlat.pro (modified) (3 diffs)

trunk/src/calc_cloud_vlat.pro

@@ -7 +7 @@
 ; .. function:: calc_cloud_vlat(jday, swr, xlat)
 ;
-; DESCRITPION
+; DESCRIPTION
 ; ===========
 ;
@@ -19 +19 @@
 ;  INPUT:
 ;       jday: value jday
-;       SWR: vector of  SW rad in Watts/m^2
+;       SWR: vector of SW rad in Watts/m^2
 ;       XLAT: corresponding vector of lat
 ;
@@ -74 +74 @@
 ;B2 = B2+(-38.79 + 2.43*XLAT - 0.034*XLAT^2)*((abs(XLAT) gt 40.) and (abs(XLAT) le 60.))
 ;
-  DAY=jday ; day is day of that year;
-  XDIV=365.15
+DAY=jday ; day is day of that year;
+XDIV=365.15
 ;
 ; CALCULATE PHI
-  PHI = (360./XDIV * (DAY-21.0))*pi/180.
+PHI = (360./XDIV * (DAY-21.0))*pi/180.
 ;
-; CALCULATE Q0; MULTIPY BY .4846 TO CONVERT FROM
+; CALCULATE Q0; MULTIPLY BY .4846 TO CONVERT FROM
 ; LY/DAY TO W/M**2
-  Q0 = A0+A1*cos(PHI)+B1*sin(PHI)+A2*cos(2*PHI)+B2*sin(2*PHI);
-  Q0=Q0*.4846;
+Q0 = A0+A1*cos(PHI)+B1*sin(PHI)+A2*cos(2*PHI)+B2*sin(2*PHI);
+Q0=Q0*.4846;
 ;
-  SNALPHA=sin(XLAT*pi/180)*sin(23.45*sin((DAY-82)*pi/180)*pi/180) + $
-           cos(XLAT*pi/180)*cos(23.45*sin((DAY-82)*pi/180)*pi/180)
-  ALPHA = asin(SNALPHA)*180/pi
+SNALPHA=sin(XLAT*pi/180)*sin(23.45*sin((DAY-82)*pi/180)*pi/180) + $
+cos(XLAT*pi/180)*cos(23.45*sin((DAY-82)*pi/180)*pi/180)
+ALPHA = asin(SNALPHA)*180/pi
 ;
-; SEE MARINE CLIMATE ATLAS OF TROPICAL PACIFC OCEAN, WEARE,STRUB,
+; SEE MARINE CLIMATE ATLAS OF TROPICAL PACIFIC OCEAN, WEARE,STRUB,
 ; AND SAMUEL -- UC,DAVIS;  THIS IS ACTUALLY REEDS 1977
 ; FORMULA
-  CLD = (1.0+0.0019*ALPHA-SWR/Q0)/.62
-  Qclsky = Q0
+CLD = (1.0+0.0019*ALPHA-SWR/Q0)/.62
+Qclsky = Q0
 ;
 indx=where(CLD ge 1.0)