Changeset 2718 for LMDZ5/trunk/libf/phylmd/surf_ocean_mod.F90

Timestamp:

Nov 29, 2016, 4:03:12 PM (8 years ago)

Author:

fhourdin

Message:

Correction de bug : (:) -> (1:knon) ou (1:klon) suivant les tableaux

File:

• LMDZ5/trunk/libf/phylmd/surf_ocean_mod.F90 (modified) (7 diffs)

LMDZ5/trunk/libf/phylmd/surf_ocean_mod.F90

@@ -104 +104 @@
 !
 !******************************************************************************
-    radsol(:) = 0.0
+    radsol(1:knon) = 0.0
     radsol(1:knon) = swnet(1:knon) + lwnet(1:knon)
 
@@ -118 +118 @@
 ! Si on suit les formulations par exemple de Tessel, on 
 ! a z0h=0.4*nu/u*, z0q=0.62*nu/u*, d'ou f_z0qh_oce=0.62/0.4=1.55
-       cdragq(:)=cdragh(:)*                                      &
-       log(z1lay(:)/z0h(:))/log(z1lay(:)/(f_z0qh_oce*z0h(:)))
+       cdragq(1:knon)=cdragh(1:knon)*                                      &
+       log(z1lay(1:knon)/z0h(1:knon))/log(z1lay(1:knon)/(f_z0qh_oce*z0h(1:knon)))
     ELSE
-       cdragq(:)=cdragh(:)
+       cdragq(1:knon)=cdragh(1:knon)
     ENDIF
 
@@ -169 +169 @@
 !******************************************************************************
     IF (type_ocean.NE.'slab') THEN
-        lmt_bils(:)=0.
+        lmt_bils(1:knon)=0.
         DO i=1,knon
            lmt_bils(knindex(i))=(swnet(i)+lwnet(i)+fluxsens(i)+fluxlat(i)) &
@@ -189 +189 @@
 !--ad-hoc correction for model radiative balance tuning
 !--now outside alboc_cd routine
-       alb_eau(:) = fmagic*alb_eau(:) + pmagic
-       alb_eau=MIN(MAX(alb_eau,0.0),1.0)
+       alb_eau(1:klon) = fmagic*alb_eau(1:klon) + pmagic
+       alb_eau(1:klon)=MIN(MAX(alb_eau(1:klon),0.0),1.0)
 !
     ELSE
@@ -197 +197 @@
 !--ad-hoc correction for model radiative balance tuning
 !--now outside alboc routine
-       alb_eau(:) = fmagic*alb_eau(:) + pmagic
-       alb_eau=MIN(MAX(alb_eau(i),0.04),0.60)
+       alb_eau(1:klon) = fmagic*alb_eau(1:klon) + pmagic
+       alb_eau(1:klon)=MIN(MAX(alb_eau(1:klon),0.04),0.60)
 !
     ENDIF
@@ -209 +209 @@
 !IM 09122015 next line corresponds to the old way of doing in LMDZ5A/IPSLCM5A versions 
 !albedo for diffuse radiation is taken the same as for direct radiation
-     alb_dif_new=alb_dir_new
+     alb_dif_new(1:knon,:)=alb_dir_new(1:knon,:)
 !IM 09122015 end
 !
@@ -219 +219 @@
 !
 !--ad-hoc correction for model radiative balance tuning
-    alb_dir_new(:,:) = fmagic*alb_dir_new(:,:) + pmagic
-    alb_dir_new=MIN(MAX(alb_dir_new,0.0),1.0)
-    alb_dif_new=MIN(MAX(alb_dif_new,0.0),1.0)
+    alb_dir_new(1:knon,:) = fmagic*alb_dir_new(1:knon,:) + pmagic
+    alb_dir_new(1:knon,:)=MIN(MAX(alb_dir_new(1:knon,:),0.0),1.0)
+    alb_dif_new(1:knon,:)=MIN(MAX(alb_dif_new(1:knon,:),0.0),1.0)
 !
 ENDIF