Changeset 203 for trunk/src/cor30a.pro

Timestamp:

02/11/14 18:41:24 (10 years ago)

Author:

pinsard

Message:

fix thanks to coding rules

File:

• trunk/src/cor30a.pro (modified) (6 diffs)

trunk/src/cor30a.pro

@@ -44 +44 @@
 ; ========
 ;
-; ::
+; .. parsed-literal::
 ;
 ;                  u  us  ts  ta   qs   qa   Qsw IRd   r  pbl  Ps   zu   zt  zq lat
@@ -50 +50 @@
 ;    x=cor30a(5.5,0,28.7,27.2,24.2,18.5,141.,419.,0.,600.,1010.,15.,15.,15.,0.,1,1,5.,1.)
 ;
-; Result with these sample values with Matlab code::
+; Result with these sample values with Matlab code:
+;
+; .. parsed-literal::
 ;
 ;      8.64830      101.640    0.0352910  2.17780e-05  0.000115000  0.000115000
@@ -57 +59 @@
 ;  0.000941350   0.00107910   0.00107910     0.780060
 ;
-; Result obtained with this idl routine::
-;
+; Result obtained with this idl routine:
+;
+; .. parsed-literal::
+
 ;      8.64829      101.640    0.0352913  2.17780e-05  0.000115000  0.000115000
 ;     -29.5802     0.175432   -0.0423667 -0.000205610     0.250948  0.000351304
@@ -223 +227 @@
   hlb=-rhoa*LLe*usr*qsr
   qout=Rnl+hsb+hlb
-  dels=Rns*(.065+11*tkt-6.6e-5/tkt*(1-exp(-tkt/8.0e-4)))        ; Eq.16 Shortwave
+  dels=Rns*(.065+11*tkt-6.6e-5/tkt*(1-exp(-tkt/8.0e-4)))      ; Eq.16 Shortwave
   qcol=qout-dels
-  alq=Al*qcol+be*hlb*cpw/LLe                                    ; Eq. 7 Buoy flux water
+  alq=Al*qcol+be*hlb*cpw/LLe                     ; Eq. 7 Buoy flux water
 
 ;  if (alq gt 0) then begin
-;    xlamx=6./(1+(bigc*alq/usr^4)^.75)^.333                     ; Eq 13 Saunders
-;    tkt=xlamx*visw/(sqrt(rhoa/rhow)*usr)                       ;Eq.11 Sub. thk
+;    xlamx=6./(1+(bigc*alq/usr^4)^.75)^.333             ; Eq 13 Saunders
+;    tkt=xlamx*visw/(sqrt(rhoa/rhow)*usr)             ;Eq.11 Sub. thk
 ;  endif else begin
 ;    xlamx=6.0
-;  tkt=min([.01,xlamx*visw/(sqrt(rhoa/rhow)*usr)])              ;Eq.11 Sub. thk
+;  tkt=min([.01,xlamx*visw/(sqrt(rhoa/rhow)*usr)])         ;Eq.11 Sub. thk
 ;    tkt=(xlamx*visw/(sqrt(rhoa/rhow)*usr) < .01)
 ;
@@ -253 +257 @@
 
 ;****************   rain heat flux ********
-dwat=2.11e-5*((t+tdk)/tdk)^1.94                         ;! water vapour diffusivity
-dtmp=(1.+3.309e-3*t-1.44e-6*t*t)*0.02411/(rhoa*cpa)     ;!heat diffusivity
-alfac= 1/(1+(wetc*LLe*dwat)/(cpa*dtmp))                 ;! wet bulb factor
+dwat=2.11e-5*((t+tdk)/tdk)^1.94             ;! water vapour diffusivity
+dtmp=(1.+3.309e-3*t-1.44e-6*t*t)*0.02411/(rhoa*cpa)      ;!heat diffusivity
+alfac= 1/(1+(wetc*LLe*dwat)/(cpa*dtmp))               ;! wet bulb factor
 RF= rain*alfac*cpw*((ts-t-dter*jcool)+(Qs-Q-dqer*jcool)*LLe/cpa)/3600.
 
@@ -262 +266 @@
 
 ;****************   Webb et al. correection  ************
-;wbar=1.61*hlb/LLe/(1+1.61*Q)/rhoa+hsb/rhoa/cpa/ta      ;formulation in hlb already includes webb
+;wbar=1.61*hlb/LLe/(1+1.61*Q)/rhoa+hsb/rhoa/cpa/ta     ;formulation in hlb already includes webb
 ;hl_webb=rhoa*wbar*Q*LLe
 ;**************   compute transfer coeffs relative to ut @meas. ht **********