Changeset 644 for LMDZ4/trunk/libf/phylmd/radlwsw.F

Timestamp:

May 25, 2005, 3:10:09 PM (19 years ago)

Author:

Laurent Fairhead

Message:

Synchronisation avec tous les diagnostiques de Ionela IM
Inclusion du slab ocean IM
LF

File:

• LMDZ4/trunk/libf/phylmd/radlwsw.F (modified) (10 diffs)

LMDZ4/trunk/libf/phylmd/radlwsw.F

@@ -89 +89 @@
       real topsw0(klon), toplw0(klon), solsw0(klon), sollw0(klon)
       real sollwdown(klon)
-cccIM 
+cIM output 3D 
       REAL*8 ZFSUP(KDLON,KFLEV+1)
       REAL*8 ZFSDN(KDLON,KFLEV+1)
       REAL*8 ZFSUP0(KDLON,KFLEV+1)
       REAL*8 ZFSDN0(KDLON,KFLEV+1)
-cIM
-cIM
-cIM   real sollwdownclr(klon) !LWdnSFCclr
-cIM   real toplwdown(klon) !LWdnTOA
-cIM   real toplwdownclr(klon) !LWdnTOAclr
+c
       REAL*8 ZFLUP(KDLON,KFLEV+1)
       REAL*8 ZFLDN(KDLON,KFLEV+1)
@@ -134 +130 @@
       REAL*8 ztopsw(kdlon), ztoplw(kdlon)
       REAL*8 zsolsw(kdlon), zsollw(kdlon), zalbpla(kdlon)
-cIM BEG
+cIM
       REAL*8 zsollwdown(kdlon)
-cIM   REAL*8 zsollwdown(kdlon), zsollwdownclr(kdlon)
-cIM   REAL*8 ztoplwdown(kdlon), ztoplwdownclr(kdlon)
-cIM END
+c
       REAL*8 ztopsw0(kdlon), ztoplw0(kdlon)
       REAL*8 zsolsw0(kdlon), zsollw0(kdlon)
       REAL*8 zznormcp
-cIM 080304   REAL swdn(klon,2),swdn0(klon,2),swup(klon,2),swup0(klon,2)
+cIM output 3D : SWup, SWdn, LWup, LWdn
       REAL swdn(klon,kflev+1),swdn0(klon,kflev+1)
       REAL swup(klon,kflev+1),swup0(klon,kflev+1)
-cIM BEG
       REAL lwdn(klon,kflev+1),lwdn0(klon,kflev+1)
       REAL lwup(klon,kflev+1),lwup0(klon,kflev+1)
-cIM END
 c-OB
 cjq the following quantities are needed for the aerosol radiative forcings
@@ -308 +300 @@
      .        ztoplw,zsollw,ztoplw0,zsollw0,
      .        zsollwdown,
-cIM  .        zsollwdown,zsollwdownclr,
-cIM  .        ztoplwdown,ztoplwdownclr)
      .        ZFLUP, ZFLDN, ZFLUP0,ZFLDN0)
 cIM ctes ds clesphys.h   CALL SW(PSCT, RCO2, zrmu0, zfract,
@@ -334 +324 @@
          sollwdown(iof+i) = zsollwdown(i)
 cIM
-cIM      sollwdownclr(iof+i) = zsollwdownclr(i)
-cIM BEG
-cIM      toplwdown(iof+i) = ztoplwdown(i)
-cIM      toplwdownclr(iof+i) = ztoplwdownclr(i)
-cIM END
-cIM 110304 BEG
          DO k = 1, kflev+1
          lwdn0 ( iof+i,k)   = ZFLDN0 ( i,k)
@@ -346 +330 @@
          lwup  ( iof+i,k)   = ZFLUP  ( i,k)
          ENDDO
-cIM 110304 END
+c
          topsw0(iof+i) = ztopsw0(i)
          toplw0(iof+i) = ztoplw0(i)
@@ -352 +336 @@
          sollw0(iof+i) = zsollw0(i)
          albpla(iof+i) = zalbpla(i)
-cIM 080304 BEG
+cIM
          DO k = 1, kflev+1
          swdn0 ( iof+i,k)   = ZFSDN0 ( i,k)
@@ -359 +343 @@
          swup  ( iof+i,k)   = ZFSUP  ( i,k)
          ENDDO !k=1, kflev+1
-cIM 080304 END
-c        swdn0 ( iof+i,1)   = ZFSDN0 ( i,1 )
-c        swdn0 ( iof+i,2)   = ZFSDN0 ( i,kflev + 1 )
-c        swdn  ( iof+i,1)   = ZFSDN  ( i,1 )
-c        swdn  ( iof+i,2)   = ZFSDN  ( i,kflev + 1 )
-c        swup0 ( iof+i,1)   = ZFSUP0 ( i,1 )
-c        swup0 ( iof+i,2)   = ZFSUP0 ( i,kflev + 1 )
-c        swup  ( iof+i,1)   = ZFSUP  ( i,1 )
-c        swup  ( iof+i,2)   = ZFSUP  ( i,kflev + 1 )
       ENDDO
 cjq-transform the aerosol forcings, if they have
@@ -2803 +2778 @@
      .              PTOPLW,PSOLLW,PTOPLW0,PSOLLW0,
      .              psollwdown,
-cIM  .              psollwdown,psollwdownclr,
-cIM  .              ptoplwdown,ptoplwdownclr)
      .              plwup, plwdn, plwup0, plwdn0)
       IMPLICIT none
@@ -2871 +2844 @@
 c Rajout LF
       real*8 psollwdown(kdlon)    ! LONGWAVE downwards flux at surface
-c Rajout IM
-cIM   real*8 psollwdownclr(kdlon) ! LONGWAVE CS downwards flux at surface
-cIM   real*8 ptoplwdown(kdlon)    ! LONGWAVE downwards flux at T.O.A.
-cIM   real*8 ptoplwdownclr(kdlon) ! LONGWAVE CS downwards flux at T.O.A.
 cIM
       REAL*8 plwup(KDLON,KFLEV+1)  ! LW up total sky
@@ -2953 +2922 @@
          PTOPLW0(i) = ZFLUC(i,1,KFLEV+1) + ZFLUC(i,2,KFLEV+1)
          psollwdown(i) = -ZFLUX(i,2,1)
-cIM
-cIM      psollwdownclr(i) = -ZFLUC(i,2,1)
-cIM      ptoplwdown(i) = ZFLUX(i,2,KFLEV+1)
-cIM      ptoplwdownclr(i) = ZFLUC(i,2,KFLEV+1)
-cIM
+c
 cIM attention aux signes !; LWtop >0, LWdn < 0
          DO k = 1, KFLEV+1