Ignore:
Timestamp:
Oct 21, 2016, 7:58:08 PM (8 years ago)
Author:
fhourdin
Message:

Rustine sur ocean_albedo pour le 1D
Correction sur la version conservative de Mellor et Yamada

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ5/trunk/libf/phylmd/ocean_albedo.F90

    r2677 r2680  
    8080!--initialisations-------------
    8181!
     82
     83IF (knon==0) RETURN ! A verifier pourquoi on en a besoin...
     84
    8285alb_dir_new(:,:) = 0.
    8386alb_dif_new(:,:) = 0.
     
    182185  ZUE=0.676               ! equivalent u_unif for diffuse incidence
    183186  ZUE2=SQRT(1.0-(1.0-ZUE**2)/ZREFM**2)
     187print*,'knon',knon
     188stop
    184189  ZRR0(1:knon)=0.50*(((ZUE2-ZREFM*ZUE)/(ZUE2+ZREFM*ZUE))**2 +((ZUE-ZREFM*ZUE2)/(ZUE+ZREFM*ZUE2))**2)
    185190  ZRRR(1:knon)=0.50*(((ZUE2-1.34*ZUE)/(ZUE2+1.34*ZUE))**2 +((ZUE-1.34*ZUE2)/(ZUE+1.34*ZUE2))**2)
Note: See TracChangeset for help on using the changeset viewer.