Changeset 1295

Timestamp:

Jun 13, 2014, 5:28:30 PM (11 years ago)

Author:

milmd

Message:

LMDZ.GENERIC. Diagnostics are done globally. They slightly differ from previous version because of some shorcuts. Correct a bug when using meanOLR option.

Location:

trunk/LMDZ.GENERIC/libf/phystd

Files:

• comgeomfi_h.F90 (modified) (1 diff)
• inifis.F (modified) (2 diffs)
• physiq.F90 (modified) (26 diffs)
• planetwide_mod.F90 (modified) (2 diffs)

trunk/LMDZ.GENERIC/libf/phystd/comgeomfi_h.F90

@@ -5 +5 @@
 
        REAL,ALLOCATABLE,DIMENSION(:) :: long,lati,area
-       REAL :: totarea
+       REAL :: totarea, totarea_planet
 
        end module comgeomfi_h

trunk/LMDZ.GENERIC/libf/phystd/inifis.F

@@ -7 +7 @@
       use datafile_mod, only: datadir
       use comdiurn_h, only: sinlat, coslat, sinlon, coslon
-      use comgeomfi_h, only: long, lati, area, totarea
+      use comgeomfi_h, only: long, lati, area, totarea, totarea_planet
       use comsoil_h, only: ini_comsoil_h
       use control_mod, only: ecritphy
+      use planetwide_mod, only: planetwide_sumval
 
 !=======================================================================
@@ -677 +678 @@
       CALL SCOPY(ngrid,parea,1,area,1)
       totarea=SSUM(ngrid,area,1)
+      call planetwide_sumval(area,totarea_planet)
 
       !! those are defined in comdiurn_h.F90

trunk/LMDZ.GENERIC/libf/phystd/physiq.F90

@@ -19 +19 @@
       use comsaison_h, only: mu0, fract, dist_star, declin
       use comsoil_h, only: nsoilmx, layer, mlayer, inertiedat
-      USE comgeomfi_h, only: long, lati, area, totarea
+      USE comgeomfi_h, only: long, lati, area, totarea, totarea_planet
       USE tracer_h, only: noms, mmol, radius, rho_q, qext, &
                           alpha_lift, alpha_devil, qextrhor, &
@@ -26 +26 @@
       use control_mod, only: ecritphy, iphysiq, nday
       use phyredem, only: physdem0, physdem1
+      use planetwide_mod, only: planetwide_minval,planetwide_maxval,planetwide_sumval
+      use mod_phys_lmdz_para, only : is_master
 
       implicit none
@@ -338 +340 @@
       real madjdE(ngrid), lscaledE(ngrid),madjdEz(ngrid,nlayermx), lscaledEz(ngrid,nlayermx)
 !JL12 conservation test for mean flow kinetic energy has been disabled temporarily
+      real dtmoist_max,dtmoist_min
       
-      real dItot, dVtot
+      real dItot, dItot_tmp, dVtot, dVtot_tmp
 
 !     included by BC for evaporation     
@@ -352 +355 @@
 !     included by RW to test water conservation (by routine)
       real h2otot
-      real dWtot, dWtots
+      real dWtot, dWtot_tmp, dWtots, dWtots_tmp
       real h2o_surf_all
       logical watertest
@@ -914 +917 @@
 ! test energy conservation
          if(enertest)then
-            dEtotSW  = cpp*SUM(massarea(:,:)*zdtsw(:,:))/totarea
-            dEtotLW  = cpp*SUM(massarea(:,:)*zdtlw(:,:))/totarea
-            dEtotsSW = SUM(fluxsurf_sw(:)*(1.-albedo(:))*area(:))/totarea !JL13 carefull, albedo can have changed since the last time we called corrk
-            dEtotsLW = SUM((fluxsurf_lw(:)*emis(:)-zplanck(:))*area(:))/totarea
+            call planetwide_sumval(cpp*massarea(:,:)*zdtsw(:,:)/totarea_planet,dEtotSW)
+            call planetwide_sumval(cpp*massarea(:,:)*zdtlw(:,:)/totarea_planet,dEtotLW)
+            call planetwide_sumval(fluxsurf_sw(:)*(1.-albedo(:))*area(:)/totarea_planet,dEtotsSW) !JL13 carefull, albedo can have changed since the last time we called corrk
+            call planetwide_sumval((fluxsurf_lw(:)*emis(:)-zplanck(:))*area(:)/totarea_planet,dEtotsLW)
             dEzRadsw(:,:)=cpp*mass(:,:)*zdtsw(:,:)
             dEzRadlw(:,:)=cpp*mass(:,:)*zdtlw(:,:)
-            print*,'---------------------------------------------------------------'
-            print*,'In corrk SW atmospheric heating       =',dEtotSW,' W m-2'
-            print*,'In corrk LW atmospheric heating       =',dEtotLW,' W m-2'
-            print*,'atmospheric net rad heating (SW+LW)   =',dEtotLW+dEtotSW,' W m-2'
-            print*,'In corrk SW surface heating           =',dEtotsSW,' W m-2'
-            print*,'In corrk LW surface heating           =',dEtotsLW,' W m-2'
-            print*,'surface net rad heating (SW+LW)       =',dEtotsLW+dEtotsSW,' W m-2'
+            if (is_master) then
+                print*,'---------------------------------------------------------------'
+                print*,'In corrk SW atmospheric heating       =',dEtotSW,' W m-2'
+                print*,'In corrk LW atmospheric heating       =',dEtotLW,' W m-2'
+                print*,'atmospheric net rad heating (SW+LW)   =',dEtotLW+dEtotSW,' W m-2'
+                print*,'In corrk SW surface heating           =',dEtotsSW,' W m-2'
+                print*,'In corrk LW surface heating           =',dEtotsLW,' W m-2'
+                print*,'surface net rad heating (SW+LW)       =',dEtotsLW+dEtotsSW,' W m-2'
+            endif
          endif
 !-------------------------
@@ -989 +994 @@
                dEzdiff(ig,1)= dEzdiff(ig,1)+ sensibFlux(ig)! subtract flux to the ground
             enddo
-            dEtot = SUM(dEdiff(:)*area(:))/totarea
+            call planetwide_sumval(dEdiff(:)*area(:)/totarea_planet,dEtot)
             dEdiffs(:)=capcal(:)*zdtsdif(:)-zflubid(:)-sensibFlux(:)
-            dEtots = SUM(dEdiffs(:)*area(:))/totarea
-            AtmToSurf_TurbFlux=SUM(sensibFlux(:)*area(:))/totarea
-            if (UseTurbDiff) then
-               print*,'In TurbDiff sensible flux (atm=>surf) =',AtmToSurf_TurbFlux,' W m-2'
-               print*,'In TurbDiff non-cons atm nrj change   =',dEtot,' W m-2'
-               print*,'In TurbDiff (correc rad+latent heat) surf nrj change =',dEtots,' W m-2'
-            else
-               print*,'In vdifc sensible flux (atm=>surf)    =',AtmToSurf_TurbFlux,' W m-2'
-               print*,'In vdifc non-cons atm nrj change      =',dEtot,' W m-2'
-               print*,'In vdifc (correc rad+latent heat) surf nrj change =',dEtots,' W m-2'
-            end if
+            call planetwide_sumval(dEdiffs(:)*area(:)/totarea_planet,dEtots)
+            call planetwide_sumval(sensibFlux(:)*area(:)/totarea_planet,AtmToSurf_TurbFlux)
+            if (is_master) then
+             if (UseTurbDiff) then
+                print*,'In TurbDiff sensible flux (atm=>surf) =',AtmToSurf_TurbFlux,' W m-2'
+                print*,'In TurbDiff non-cons atm nrj change   =',dEtot,' W m-2'
+                print*,'In TurbDiff (correc rad+latent heat) surf nrj change =',dEtots,' W m-2'
+             else
+                print*,'In vdifc sensible flux (atm=>surf)    =',AtmToSurf_TurbFlux,' W m-2'
+                print*,'In vdifc non-cons atm nrj change      =',dEtot,' W m-2'
+                print*,'In vdifc (correc rad+latent heat) surf nrj change =',dEtots,' W m-2'
+             end if
+            endif ! of if (is_master)
 ! JL12 note that the black body radiative flux emitted by the surface has been updated by the implicit scheme
 !    but not given back elsewhere
@@ -1010 +1017 @@
          ! test water conservation
          if(watertest.and.water)then
-            dWtot = SUM(massarea(:,:)*zdqdif(:,:,igcm_h2o_vap))*ptimestep/totarea
-            dWtots =  SUM(zdqsdif(:,igcm_h2o_vap)*area(:))*ptimestep/totarea
+            call planetwide_sumval(massarea(:,:)*zdqdif(:,:,igcm_h2o_vap)*ptimestep/totarea_planet,dWtot_tmp)
+            call planetwide_sumval(zdqsdif(:,igcm_h2o_vap)*area(:)*ptimestep/totarea_planet,dWtots_tmp)
             do ig = 1, ngrid
                vdifcncons(ig)=SUM(mass(ig,:)*zdqdif(ig,:,igcm_h2o_vap))
             Enddo
-            dWtot = dWtot + SUM(massarea(:,:)*zdqdif(:,:,igcm_h2o_ice))*ptimestep/totarea
-            dWtots = dWtots + SUM(zdqsdif(:,igcm_h2o_ice)*area(:))*ptimestep/totarea
+            call planetwide_sumval(massarea(:,:)*zdqdif(:,:,igcm_h2o_ice)*ptimestep/totarea_planet,dWtot)
+            call planetwide_sumval(zdqsdif(:,igcm_h2o_ice)*area(:)*ptimestep/totarea_planet,dWtots)
+            dWtot = dWtot + dWtot_tmp
+            dWtots = dWtots + dWtots_tmp
             do ig = 1, ngrid
                vdifcncons(ig)=vdifcncons(ig) + SUM(mass(ig,:)*zdqdif(ig,:,igcm_h2o_ice))
             Enddo           
-            nconsMAX=MAXVAL(vdifcncons(:))
-
-            print*,'---------------------------------------------------------------'
-            print*,'In difv atmospheric water change        =',dWtot,' kg m-2'
-            print*,'In difv surface water change            =',dWtots,' kg m-2'
-            print*,'In difv non-cons factor                 =',dWtot+dWtots,' kg m-2'
-            print*,'In difv MAX non-cons factor             =',nconsMAX,' kg m-2 s-1'
+            call planetwide_maxval(vdifcncons(:),nconsMAX)
+
+            if (is_master) then
+                print*,'---------------------------------------------------------------'
+                print*,'In difv atmospheric water change        =',dWtot,' kg m-2'
+                print*,'In difv surface water change            =',dWtots,' kg m-2'
+                print*,'In difv non-cons factor                 =',dWtot+dWtots,' kg m-2'
+                print*,'In difv MAX non-cons factor             =',nconsMAX,' kg m-2 s-1'
+            endif
 
          endif
@@ -1073 +1084 @@
          ! test energy conservation
          if(enertest)then
-            dEtot=cpp*SUM(massarea(:,:)*zdtadj(:,:))/totarea
-          print*,'In convadj atmospheric energy change  =',dEtot,' W m-2'
+            call planetwide_sumval(cpp*massarea(:,:)*zdtadj(:,:)/totarea_planet,dEtot)
+            if (is_master) print*,'In convadj atmospheric energy change  =',dEtot,' W m-2'
          endif
          !-------------------------
@@ -1081 +1092 @@
          ! test water conservation
          if(watertest)then
-            dWtot = SUM(massarea(:,:)*zdqadj(:,:,igcm_h2o_vap))*ptimestep/totarea
+            call planetwide_sumval(massarea(:,:)*zdqadj(:,:,igcm_h2o_vap)*ptimestep/totarea_planet,dWtot_tmp)
             do ig = 1, ngrid
                cadjncons(ig)=SUM(mass(ig,:)*zdqadj(ig,:,igcm_h2o_vap))
             Enddo
-            dWtot = dWtot + SUM(massarea(:,:)*zdqadj(:,:,igcm_h2o_ice))*ptimestep/totarea
+            call planetwide_sumval(massarea(:,:)*zdqadj(:,:,igcm_h2o_ice)*ptimestep/totarea_planet,dWtot)
+            dWtot = dWtot + dWtot_tmp
             do ig = 1, ngrid
                cadjncons(ig)=cadjncons(ig) + SUM(mass(ig,:)*zdqadj(ig,:,igcm_h2o_ice))
             Enddo           
-            nconsMAX=MAXVAL(cadjncons(:))
-
-            print*,'In convadj atmospheric water change     =',dWtot,' kg m-2'
-            print*,'In convadj MAX non-cons factor          =',nconsMAX,' kg m-2 s-1'
+            call planetwide_maxval(cadjncons(:),nconsMAX)
+
+            if (is_master) then
+                print*,'In convadj atmospheric water change     =',dWtot,' kg m-2'
+                print*,'In convadj MAX non-cons factor          =',nconsMAX,' kg m-2 s-1'
+            endif
          endif
          !-------------------------
@@ -1126 +1140 @@
          ! test energy conservation
          if(enertest)then
-            dEtot = cpp*SUM(massarea(:,:)*zdtc(:,:))/totarea
-            dEtots = SUM(capcal(:)*zdtsurfc(:)*area(:))/totarea
-            print*,'In co2cloud atmospheric energy change   =',dEtot,' W m-2'
-            print*,'In co2cloud surface energy change       =',dEtots,' W m-2'
+            call planetwide_sumval(cpp*massarea(:,:)*zdtc(:,:)/totarea_planet,dEtot)
+            call planetwide_sumval(capcal(:)*zdtsurfc(:)*area(:)/totarea_planet,dEtots)
+            if (is_master) then
+                print*,'In co2cloud atmospheric energy change   =',dEtot,' W m-2'
+                print*,'In co2cloud surface energy change       =',dEtots,' W m-2'
+            endif
          endif
          !-------------------------
@@ -1169 +1185 @@
                ! test energy conservation
                if(enertest)then
-                  dEtot=cpp*SUM(massarea(:,:)*dtmoist(:,:))/totarea
+                  call planetwide_sumval(cpp*massarea(:,:)*dtmoist(:,:)/totarea_planet,dEtot)
+                  call planetwide_maxval(dtmoist(:,:),dtmoist_max)
+                  call planetwide_minval(dtmoist(:,:),dtmoist_min)
                   madjdEz(:,:)=cpp*mass(:,:)*dtmoist(:,:)
                   do ig=1,ngrid
                      madjdE(ig) = cpp*SUM(mass(:,:)*dtmoist(:,:))
                   enddo
-                  print*,'In moistadj atmospheric energy change   =',dEtot,' W m-2'
-                  print*,'In moistadj MAX atmospheric energy change   =',MAXVAL(dtmoist(:,:))*ptimestep,'K/step'
-                  print*,'In moistadj MIN atmospheric energy change   =',MINVAL(dtmoist(:,:))*ptimestep,'K/step'
+                  if (is_master) then
+                        print*,'In moistadj atmospheric energy change   =',dEtot,' W m-2'
+                        print*,'In moistadj MAX atmospheric energy change   =',dtmoist_max*ptimestep,'K/step'
+                        print*,'In moistadj MIN atmospheric energy change   =',dtmoist_min*ptimestep,'K/step'
+                  endif
                   
                 ! test energy conservation
-                  dWtot = SUM(massarea(:,:)*dqmoist(:,:,igcm_h2o_vap))*ptimestep/totarea
-                  dWtot = dWtot + SUM(massarea(:,:)*dqmoist(:,:,igcm_h2o_ice))*ptimestep/totarea
-                  print*,'In moistadj atmospheric water change    =',dWtot,' kg m-2'
+                  call planetwide_sumval(massarea(:,:)*dqmoist(:,:,igcm_h2o_vap)*ptimestep/totarea_planet+      &
+                        massarea(:,:)*dqmoist(:,:,igcm_h2o_ice)*ptimestep/totarea_planet,dWtot)
+                  if (is_master) print*,'In moistadj atmospheric water change    =',dWtot,' kg m-2'
                endif
                !-------------------------
@@ -1202 +1222 @@
                      lscaledE(ig) = cpp*SUM(mass(:,:)*dtlscale(:,:))
                   enddo
-                  dEtot=cpp*SUM(massarea(:,:)*(dtlscale(:,:)))/totarea
+                  call planetwide_sumval(cpp*massarea(:,:)*dtlscale(:,:)/totarea_planet,dEtot)
 !                 if(isnan(dEtot)) then ! NB: isnan() is not a standard function...
 !                    print*,'Nan in largescale, abort'
 !                     STOP
 !                 endif
-                  print*,'In largescale atmospheric energy change =',dEtot,' W m-2'
+                  if (is_master) print*,'In largescale atmospheric energy change =',dEtot,' W m-2'
 
                ! test water conservation
-                  dWtot = SUM(massarea(:,:)*dqvaplscale(:,:))*ptimestep/totarea
-                  dWtot = dWtot + SUM(massarea(:,:)*dqcldlscale(:,:))*ptimestep/totarea
-                  print*,'In largescale atmospheric water change  =',dWtot,' kg m-2'
+                  call planetwide_sumval(massarea(:,:)*dqvaplscale(:,:)*ptimestep/totarea_planet+       &
+                        SUM(massarea(:,:)*dqcldlscale(:,:))*ptimestep/totarea_planet,dWtot)
+                  if (is_master) print*,'In largescale atmospheric water change  =',dWtot,' kg m-2'
                endif
                !-------------------------
@@ -1250 +1270 @@
                ! test energy conservation
                if(enertest)then
-                  dEtot=cpp*SUM(massarea(:,:)*zdtrain(:,:))/totarea
-                 print*,'In rain atmospheric T energy change       =',dEtot,' W m-2'
-                  dItot = SUM(massarea(:,:)*zdqrain(:,:,igcm_h2o_ice))/totarea*RLVTT/cpp
-                  dItot = dItot + SUM(area(:)*zdqssnow(:))/totarea*RLVTT/cpp
-                  dVtot = SUM(massarea(:,:)*zdqrain(:,:,igcm_h2o_vap))*ptimestep/totarea
-                  dVtot = dVtot + SUM(area(:)*zdqsrain(:))/totarea*RLVTT/cpp
+                  call planetwide_sumval(cpp*massarea(:,:)*zdtrain(:,:)/totarea_planet,dEtot)
+                  if (is_master) print*,'In rain atmospheric T energy change       =',dEtot,' W m-2'
+                  call planetwide_sumval(massarea(:,:)*zdqrain(:,:,igcm_h2o_ice)/totarea_planet*RLVTT/cpp,dItot_tmp)
+                  call planetwide_sumval(area(:)*zdqssnow(:)/totarea_planet*RLVTT/cpp,dItot)
+                  dItot = dItot + dItot_tmp
+                  call planetwide_sumval(massarea(:,:)*zdqrain(:,:,igcm_h2o_vap)*ptimestep/totarea_planet,dVtot_tmp)
+                  call planetwide_sumval(area(:)*zdqsrain(:)/totarea_planet*RLVTT/cpp,dVtot)
+                  dVtot = dVtot + dVtot_tmp
                   dEtot = dItot + dVtot
-                 print*,'In rain dItot =',dItot,' W m-2'
-                 print*,'In rain dVtot =',dVtot,' W m-2'
-                 print*,'In rain atmospheric L energy change       =',dEtot,' W m-2'
+                  if (is_master) then
+                        print*,'In rain dItot =',dItot,' W m-2'
+                        print*,'In rain dVtot =',dVtot,' W m-2'
+                        print*,'In rain atmospheric L energy change       =',dEtot,' W m-2'
+                  endif
 
                ! test water conservation
-                  dWtot = SUM(massarea(:,:)*zdqrain(:,:,igcm_h2o_vap))*ptimestep/totarea
-                  dWtot = dWtot + SUM(massarea(:,:)*zdqrain(:,:,igcm_h2o_ice))*ptimestep/totarea
-                  dWtots =  SUM((zdqsrain(:)+zdqssnow(:))*area(:))*ptimestep/totarea
-                 print*,'In rain atmospheric water change        =',dWtot,' kg m-2'
-                 print*,'In rain surface water change            =',dWtots,' kg m-2'
-                 print*,'In rain non-cons factor                 =',dWtot+dWtots,' kg m-2'
+                  call planetwide_sumval(massarea(:,:)*zdqrain(:,:,igcm_h2o_vap)*ptimestep/totarea_planet+      &
+                        massarea(:,:)*zdqrain(:,:,igcm_h2o_ice)*ptimestep/totarea_planet,dWtot)
+                  call planetwide_sumval((zdqsrain(:)+zdqssnow(:))*area(:)*ptimestep/totarea_planet,dWtots)
+                  if (is_master) then
+                        print*,'In rain atmospheric water change        =',dWtot,' kg m-2'
+                        print*,'In rain surface water change            =',dWtots,' kg m-2'
+                        print*,'In rain non-cons factor                 =',dWtot+dWtots,' kg m-2'
+                  endif
               endif
               !-------------------------
@@ -1289 +1315 @@
            if(watertest)then
               iq=igcm_h2o_ice
-              dWtot = SUM(massarea(:,:)*pq(:,:,iq))*ptimestep/totarea
-              dWtots = SUM(massarea(:,:)*pdq(:,:,iq))*ptimestep/totarea
-              print*,'Before sedim pq  =',dWtot,' kg m-2'
-              print*,'Before sedim pdq =',dWtots,' kg m-2'
+              call planetwide_sumval(massarea(:,:)*pq(:,:,iq)*ptimestep/totarea_planet,dWtot)
+              call planetwide_sumval(massarea(:,:)*pdq(:,:,iq)*ptimestep/totarea_planet,dWtots)
+              if (is_master) then
+                print*,'Before sedim pq  =',dWtot,' kg m-2'
+                print*,'Before sedim pdq =',dWtots,' kg m-2'
+              endif
            endif
            !-------------------------
@@ -1303 +1331 @@
            if(watertest)then
               iq=igcm_h2o_ice
-              dWtot = SUM(massarea(:,:)*pq(:,:,iq))*ptimestep/totarea
-              dWtots = SUM(massarea(:,:)*pdq(:,:,iq))*ptimestep/totarea
-              print*,'After sedim pq  =',dWtot,' kg m-2'
-              print*,'After sedim pdq =',dWtots,' kg m-2'
+              call planetwide_sumval(massarea(:,:)*pq(:,:,iq)*ptimestep/totarea_planet,dWtot)
+              call planetwide_sumval(massarea(:,:)*pdq(:,:,iq)*ptimestep/totarea_planet,dWtots)
+              if (is_master) then
+                print*,'After sedim pq  =',dWtot,' kg m-2'
+                print*,'After sedim pdq =',dWtots,' kg m-2'
+              endif
            endif
            !-------------------------
@@ -1319 +1349 @@
            ! test water conservation
            if(watertest)then
-              dWtot = SUM(massarea(:,:)*(zdqsed(:,:,igcm_h2o_vap)+zdqsed(:,:,igcm_h2o_ice)))*ptimestep/totarea
-              dWtots = SUM((zdqssed(:,igcm_h2o_vap)+zdqssed(:,igcm_h2o_ice))*area(:))*ptimestep/totarea
-              print*,'In sedim atmospheric ice change         =',dWtot,' kg m-2'
-              print*,'In sedim surface ice change             =',dWtots,' kg m-2'
-              print*,'In sedim non-cons factor                =',dWtot+dWtots,' kg m-2'
+              call planetwide_sumval(massarea(:,:)*(zdqsed(:,:,igcm_h2o_vap)+zdqsed(:,:,igcm_h2o_ice))*ptimestep/totarea_planet,dWtot)
+              call planetwide_sumval((zdqssed(:,igcm_h2o_vap)+zdqssed(:,igcm_h2o_ice))*area(:)*ptimestep/totarea_planet,dWtots)
+              if (is_master) then
+                print*,'In sedim atmospheric ice change         =',dWtot,' kg m-2'
+                print*,'In sedim surface ice change             =',dWtots,' kg m-2'
+                print*,'In sedim non-cons factor                =',dWtot+dWtots,' kg m-2'
+              endif
            endif
            !-------------------------
@@ -1389 +1421 @@
             ! test energy conservation
             if(enertest)then
-               dEtots = SUM(area(:)*capcal(:)*zdtsurf_hyd(:))/totarea
-               print*,'In hydrol surface energy change     =',dEtots,' W m-2'
+               call planetwide_sumval(area(:)*capcal(:)*zdtsurf_hyd(:)/totarea_planet,dEtots)
+               if (is_master) print*,'In hydrol surface energy change     =',dEtots,' W m-2'
             endif
             !-------------------------
@@ -1397 +1429 @@
             ! test water conservation
             if(watertest)then
-               dWtots =  SUM(dqs_hyd(:,igcm_h2o_ice)*area(:))*ptimestep/totarea
-               print*,'In hydrol surface ice change            =',dWtots,' kg m-2'
-               dWtots =  SUM(dqs_hyd(:,igcm_h2o_vap)*area(:))*ptimestep/totarea
-               print*,'In hydrol surface water change          =',dWtots,' kg m-2'
-               print*,'---------------------------------------------------------------'
+               call planetwide_sumval(dqs_hyd(:,igcm_h2o_ice)*area(:)*ptimestep/totarea_planet,dWtots)
+               if (is_master) print*,'In hydrol surface ice change            =',dWtots,' kg m-2'
+               call planetwide_sumval(dqs_hyd(:,igcm_h2o_vap)*area(:)*ptimestep/totarea_planet,dWtots)
+               if (is_master) then
+                print*,'In hydrol surface water change          =',dWtots,' kg m-2'
+                print*,'---------------------------------------------------------------'
+               endif
             endif
             !-------------------------
@@ -1439 +1473 @@
 ! test energy conservation
       if(enertest)then
-         dEtots = SUM(area(:)*capcal(:)*zdtsurf(:))/totarea
-         print*,'Surface energy change                 =',dEtots,' W m-2'
+         call planetwide_sumval(area(:)*capcal(:)*zdtsurf(:)/totarea_planet,dEtots)      
+         if (is_master) print*,'Surface energy change                 =',dEtots,' W m-2'
       endif
 !-------------------------
@@ -1487 +1521 @@
 !     ---------------------------------------------------------
 
-      Ts1 = SUM(area(:)*tsurf(:))/totarea
-      Ts2 = MINVAL(tsurf(:))
-      Ts3 = MAXVAL(tsurf(:))
+      call planetwide_sumval(area(:)*tsurf(:)/totarea_planet,Ts1)
+      call planetwide_minval(tsurf(:),Ts2)
+      call planetwide_maxval(tsurf(:),Ts3)
       if(callsoil)then
          TsS = SUM(area(:)*tsoil(:,nsoilmx))/totarea        ! mean temperature at bottom soil layer
@@ -1495 +1529 @@
            print*,Ts1,Ts2,Ts3,TsS
       else
+        if (is_master) then
            print*,'  ave[Tsurf]     min[Tsurf]     max[Tsurf]'
            print*,Ts1,Ts2,Ts3
+        endif
       end if
 
@@ -1505 +1541 @@
       if(corrk)then
 
-         ISR = SUM(area(:)*fluxtop_dn(:))/totarea
-         ASR = SUM(area(:)*fluxabs_sw(:))/totarea
-         OLR = SUM(area(:)*fluxtop_lw(:))/totarea
-         GND = SUM(area(:)*fluxgrd(:))/totarea
-         DYN = SUM(area(:)*fluxdyn(:))/totarea
+         call planetwide_sumval(area(:)*fluxtop_dn(:)/totarea_planet,ISR)
+         call planetwide_sumval(area(:)*fluxabs_sw(:)/totarea_planet,ASR)
+         call planetwide_sumval(area(:)*fluxtop_lw(:)/totarea_planet,OLR)
+         call planetwide_sumval(area(:)*fluxgrd(:)/totarea_planet,GND)
+         call planetwide_sumval(area(:)*fluxdyn(:)/totarea_planet,DYN)
          do ig=1,ngrid
             if(fluxtop_dn(ig).lt.0.0)then
@@ -1525 +1561 @@
             DYN=0.0
          endif
-
-         print*,'  ISR            ASR            OLR            GND            DYN [W m^-2]'
-         print*, ISR,ASR,OLR,GND,DYN
-
-         if(enertest)then
+         
+         if (is_master) then
+                print*,'  ISR            ASR            OLR            GND            DYN [W m^-2]'
+                print*, ISR,ASR,OLR,GND,DYN
+         endif
+
+         if(enertest .and. is_master)then
             print*,'SW flux/heating difference SW++ - ASR = ',dEtotSW+dEtotsSW-ASR,' W m-2'
             print*,'LW flux/heating difference LW++ - OLR = ',dEtotLW+dEtotsLW+OLR,' W m-2'
@@ -1535 +1573 @@
          endif
 
-         if(meanOLR)then
+         if(meanOLR .and. is_master)then
             if((ngrid.gt.1) .or. (mod(icount-1,ecritphy).eq.0))then
                ! to record global radiative balance
@@ -1542 +1580 @@
                close(92)
                open(93,file="tem_bal.out",form='formatted',position='append')
-               write(93,*) zday,Ts1,Ts2,Ts3,TsS
+               if(callsoil)then
+                write(93,*) zday,Ts1,Ts2,Ts3,TsS
+               else
+                write(93,*) zday,Ts1,Ts2,Ts3
+               endif
                close(93)
             endif
@@ -1600 +1642 @@
       if(water)then
 
-         icesrf = SUM(area(:)*qsurf_hist(:,igcm_h2o_ice))/totarea
-         liqsrf = SUM(area(:)*qsurf_hist(:,igcm_h2o_vap))/totarea
-         icecol = SUM(area(:)*qcol(:,igcm_h2o_ice))/totarea
-         vapcol = SUM(area(:)*qcol(:,igcm_h2o_vap))/totarea
+         call planetwide_sumval(area(:)*qsurf_hist(:,igcm_h2o_ice)/totarea_planet,icesrf)
+         call planetwide_sumval(area(:)*qsurf_hist(:,igcm_h2o_vap)/totarea_planet,liqsrf)
+         call planetwide_sumval(area(:)*qcol(:,igcm_h2o_ice)/totarea_planet,icecol)
+         call planetwide_sumval(area(:)*qcol(:,igcm_h2o_vap)/totarea_planet,vapcol)
 
          h2otot = icesrf + liqsrf + icecol + vapcol
-
-         print*,' Total water amount [kg m^-2]: ',h2otot
-         print*,' Surface ice    Surface liq.   Atmos. con.     Atmos. vap. [kg m^-2] '
-         print*, icesrf,liqsrf,icecol,vapcol
-
-         if(meanOLR)then
+         
+         if (is_master) then
+                print*,' Total water amount [kg m^-2]: ',h2otot
+                print*,' Surface ice    Surface liq.   Atmos. con.     Atmos. vap. [kg m^-2] '
+                print*, icesrf,liqsrf,icecol,vapcol
+         endif
+
+         if(meanOLR .and. is_master)then
             if((ngrid.gt.1) .or. (mod(icount-1,ecritphy).eq.0))then
                ! to record global water balance

trunk/LMDZ.GENERIC/libf/phystd/planetwide_mod.F90

@@ -20 +20 @@
 end interface
 
+interface planetwide_sumval ! sum() , over the entire planet
+  module procedure planetwide_sumval_i1, planetwide_sumval_i2, &
+                   planetwide_sumval_r1, planetwide_sumval_r2
+end interface
+
 contains
 
@@ -200 +205 @@
   end subroutine planetwide_minval_r2
 
-
+!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
+
+  subroutine planetwide_sumval_i1(values,values_sum)
+  use dimphy, only: klon
+  use mod_grid_phy_lmdz, only : klon_glo
+  implicit none
+  integer,intent(in) :: values(:) ! local grid (klon)
+  integer,intent(out) :: values_sum
+#ifdef CPP_PARA
+  integer :: values_glo(klon_glo) ! global grid
+  
+  ! gather field on master:
+  call gather(values,values_glo)
+  ! calculate sum value
+  if (is_master) then
+    values_sum=SUM(values_glo(:))
+  endif
+  ! broadcast information to all cores
+  call bcast(values_sum)
+#else
+  values_sum=SUM(values(:))
+#endif
+  end subroutine planetwide_sumval_i1
+
+  subroutine planetwide_sumval_i2(values,values_sum)
+  use dimphy, only: klon, klev
+  use mod_grid_phy_lmdz, only : klon_glo
+  implicit none
+  integer,intent(in) :: values(:,:) ! local grid (klon,klev)
+  integer,intent(out) :: values_sum
+#ifdef CPP_PARA
+  integer :: values_glo(klon_glo,klev) ! global grid
+  
+  ! gather field on master:
+  call gather(values,values_glo)
+  ! calculate sum value
+  if (is_master) then
+    values_sum=SUM(values_glo)
+  endif
+  ! broadcast information to all cores
+  call bcast(values_sum)
+#else
+  values_sum=SUM(values)
+#endif
+  end subroutine planetwide_sumval_i2
+
+  subroutine planetwide_sumval_r1(values,values_sum)
+  use dimphy, only: klon
+  use mod_grid_phy_lmdz, only : klon_glo
+  implicit none
+  real,intent(in) :: values(:) ! local grid (klon)
+  real,intent(out) :: values_sum
+#ifdef CPP_PARA
+  real :: values_glo(klon_glo) ! global grid
+  
+  ! gather field on master:
+  call gather(values,values_glo)
+  ! calculate sum value
+  if (is_master) then
+    values_sum=SUM(values_glo)
+  endif
+  ! broadcast information to all cores
+  call bcast(values_sum)
+#else
+  values_sum=SUM(values)
+#endif
+  end subroutine planetwide_sumval_r1
+
+  subroutine planetwide_sumval_r2(values,values_sum)
+  use dimphy, only: klon, klev
+  use mod_grid_phy_lmdz, only : klon_glo
+  implicit none
+  real,intent(in) :: values(:,:) ! local grid (klon,klev)
+  real,intent(out) :: values_sum
+#ifdef CPP_PARA
+  real :: values_glo(klon_glo,klev) ! global grid
+  
+  ! gather field on master:
+  call gather(values,values_glo)
+  ! calculate sum value
+  if (is_master) then
+    values_sum=SUM(values_glo)
+  endif
+  ! broadcast information to all cores
+  call bcast(values_sum)
+#else
+  values_sum=SUM(values)
+#endif
+  end subroutine planetwide_sumval_r2
+  
+  
 end module planetwide_mod