/[lmdze]/trunk/Sources/phylmd/readsulfate_preind.f

Diff of /trunk/Sources/phylmd/readsulfate_preind.f

Parent Directory | Revision Log | View Patch Patch

-revision 190 by guez,
Tue Mar 22 16:31:39 2016 UTC
+revision 191 by guez,
Mon May  9 19:56:28 2016 UTC
 Line 4 
 module readsulfate_preind_m
  contains
-   SUBROUTINE readsulfate_preind(dayvrai, time, first, pi_sulfate)
+   SUBROUTINE readsulfate_preind(dayvrai, time, first, sulfate)
-     ! Read in /calculate pre-industrial values of sulfate
+     ! Read and calculate pre-industrial values of sulfate. This
+     ! routine reads monthly mean values of sulfate aerosols and
+     ! returns a linearly interpolated daily-mean field. It does so for
+     ! the preindustriel values of the sulfate, to a large part
+     ! analogous to the routine readsulfate.
+     ! Author:
+     ! Johannes Quaas (quaas@lmd.jussieu.fr)
+     ! June 26th, 2001
-     use dimens_m
+     use dimens_m, only: iim, jjm
-     use dimphy
+     use dimphy, only: klon, klev
      use dynetat0_m, only: annee_ref
-     use SUPHEC_M
-     use chem
      use getso4fromfile_m, only: getso4fromfile
-     ! Content:
-     ! --------
-     ! This routine reads in monthly mean values of sulfate aerosols and
-     ! returns a linearly interpolated daily-mean field.
-     !
-     ! It does so for the preindustriel values of the sulfate, to a large part
-     ! analogous to the routine readsulfate.
-     !
-     ! Only Pb: Variables must be saved and don t have to be overwritten!
-     !
-     ! Author:
-     ! -------
-     ! Johannes Quaas (quaas@lmd.jussieu.fr)
-     ! 26/06/01
-     !
-     ! Input:
-     ! ------
      integer, intent(in):: dayvrai
      ! current day number, based at value 1 on January 1st of annee_ref
      REAL, intent(in):: time ! heure de la journ\'ee en fraction de jour
-     LOGICAL, intent(in):: first                 ! First timestep
+     LOGICAL, intent(in):: first ! First timestep
      ! (and therefore initialization necessary)
-     !
-     ! Output:
+     real, intent(out):: sulfate (klon, klev)
-     ! -------
+     ! number concentration sulfate (monthly mean data, from file)
-     real pi_sulfate (klon, klev)  ! Number conc. sulfate (monthly mean data,
-     !  from file)
+     ! Local:
-     !
-     ! Local Variables:
-     ! ----------------
      INTEGER i, ig, k, it
      INTEGER j, iday, iyr
      INTEGER im, day1, day2, im2
-     double precision pi_so4_1(iim, jjm+1, klev, 12)
+     real so4_1(iim, jjm + 1, klev, 12)
+     double precision, save:: so4(klon, klev, 12) ! SO4 in right dimension
-     double precision pi_so4(klon, klev, 12)  ! SO4 in right dimension
+     double precision, save:: so4_out(klon, klev)
-     SAVE pi_so4
-     double precision pi_so4_out(klon, klev)
-     SAVE pi_so4_out
-     CHARACTER(len=4) cyear
      LOGICAL lnewday
+     !----------------------------------------------------------------
      iday = dayvrai
      ! Get the year of the run
-     iyr  = iday/360
+     iyr = iday/360
      ! Get the day of the actual year:
-     iday = iday-iyr*360
+     iday = iday - iyr*360
      ! 0.02 is about 0.5/24, namly less than half an hour
      lnewday = time < 0.02
-     ! ---------------------------------------------
+     ! All has to be done only, if a new day begins
-     ! All has to be done only, if a new day begins!
-     ! ---------------------------------------------
-     IF (lnewday.OR.first) THEN
+     test_newday: IF (lnewday .OR. first) THEN
-        im = iday/30 +1 ! the actual month
+        im = iday/30 + 1 ! the actual month
         ! annee_ref is the initial year (defined in temps.h)
         iyr = iyr + annee_ref
         IF (first) THEN
-           cyear='.nat'
-           CALL getso4fromfile(cyear,pi_so4_1)
-           ! Transform the horizontal 2D-field into the physics-field
-           ! (Also the levels and the latitudes have to be inversed)
            ! Initialize field
-           DO it=1,12
+           DO it = 1, 12
-              DO k=1,klev
+              DO k = 1, klev
-                 DO i=1,klon
+                 DO i = 1, klon
-                    pi_so4(i,k,it)=0.
+                    so4(i, k, it) = 0.
                  ENDDO
               ENDDO
            ENDDO
-           write (*,*) 'preind: finished reading', FLOAT(iim)
+           CALL getso4fromfile('.nat', so4_1)
-           DO it=1,12
-              DO k=1, klev
+           ! Transform the horizontal 2D-field into the physics-field
+           ! (Also the levels and the latitudes have to be inversed)
+           DO it = 1, 12
+              DO k = 1, klev
                  ! a) at the poles, use the zonal mean:
-                 DO i=1,iim
+                 DO i = 1, iim
                     ! North pole
-                    pi_so4(1,k,it)=pi_so4(1,k,it)+pi_so4_1(i,jjm+1,klev+1-k,it)
+                    so4(1, k, it) = so4(1, k, it) &
+                         + so4_1(i, jjm + 1, klev + 1 - k, it)
                     ! South pole
-                    pi_so4(klon,k,it)=pi_so4(klon,k,it)+pi_so4_1(i,1,klev+1-k,it)
+                    so4(klon, k, it) = so4(klon, k, it) &
+                         + so4_1(i, 1, klev + 1 - k, it)
                  ENDDO
-                 pi_so4(1,k,it)=pi_so4(1,k,it)/FLOAT(iim)
+                 so4(1, k, it) = so4(1, k, it)/REAL(iim)
-                 pi_so4(klon,k,it)=pi_so4(klon,k,it)/FLOAT(iim)
+                 so4(klon, k, it) = so4(klon, k, it)/REAL(iim)
                  ! b) the values between the poles:
-                 ig=1
+                 ig = 1
-                 DO j=2,jjm
+                 DO j = 2, jjm
-                    DO i=1,iim
+                    DO i = 1, iim
-                       ig=ig+1
+                       ig = ig + 1
-                       if (ig.gt.klon) write (*,*) 'shit'
+                       if (ig > klon) stop 1
-                       pi_so4(ig,k,it) = pi_so4_1(i,jjm+1-j,klev+1-k,it)
+                       so4(ig, k, it) = so4_1(i, jjm + 1 - j, klev + 1 - k, it)
                     ENDDO
                  ENDDO
-                 IF (ig.NE.klon-1) STOP 'Error in readsulfate (var conversion)'
+                 IF (ig /= klon - 1) then
+                    print *, 'Error in readsulfate (var conversion)'
+                    STOP 1
+                 end IF
               ENDDO ! Loop over k (vertical)
            ENDDO ! Loop over it (months)
+        ENDIF ! Had to read new data?
-        ENDIF                     ! Had to read new data?
         ! Interpolate to actual day:
-        IF (iday.LT.im*30-15) THEN
+        IF (iday < im*30 - 15) THEN
            ! in the first half of the month use month before and actual month
-           im2=im-1
+           im2 = im - 1
-           day1 = im2*30+15
+           day1 = im2*30 + 15
-           day2 = im2*30-15
+           day2 = im2*30 - 15
-           IF (im2.LE.0) THEN
+           IF (im2 <= 0) THEN
               ! the month is january, thus the month before december
-              im2=12
+              im2 = 12
            ENDIF
-           DO k=1,klev
+           DO k = 1, klev
-              DO i=1,klon
+              DO i = 1, klon
-                 pi_sulfate(i,k) = pi_so4(i,k,im2)   &
+                 sulfate(i, k) = so4(i, k, im2) &
-                      - FLOAT(iday-day2)/FLOAT(day1-day2) &
+                      - REAL(iday - day2)/REAL(day1 - day2) &
-                      * (pi_so4(i,k,im2) - pi_so4(i,k,im))
+                      * (so4(i, k, im2) - so4(i, k, im))
-                 IF (pi_sulfate(i,k).LT.0.) THEN
+                 IF (sulfate(i, k) < 0.) THEN
-                    IF (iday-day2.LT.0.) write(*,*) 'iday-day2',iday-day2
+                    IF (iday - day2 < 0.) write(*, *) 'iday - day2', iday - day2
-                    IF (pi_so4(i,k,im2) - pi_so4(i,k,im).LT.0.) &
+                    IF (so4(i, k, im2) - so4(i, k, im) < 0.) &
-                         write(*,*) 'pi_so4(i,k,im2) - pi_so4(i,k,im)', &
+                         write(*, *) 'so4(i, k, im2) - so4(i, k, im)', &
-                         pi_so4(i,k,im2) - pi_so4(i,k,im)
+                         so4(i, k, im2) - so4(i, k, im)
-                    IF (day1-day2.LT.0.) write(*,*) 'day1-day2',day1-day2
+                    IF (day1 - day2 < 0.) write(*, *) 'day1 - day2', day1 - day2
-                    stop 'pi_sulfate'
+                    stop 1
                  endif
               ENDDO
            ENDDO
         ELSE
            ! the second half of the month
-           im2=im+1
+           im2 = im + 1
-           day1 = im*30+15
+           IF (im2 > 12) THEN
-           IF (im2.GT.12) THEN
               ! the month is december, the following thus january
-              im2=1
+              im2 = 1
            ENDIF
-           day2 = im*30-15
+           day2 = im*30 - 15
+           day1 = im*30 + 15
-           DO k=1,klev
+           DO k = 1, klev
-              DO i=1,klon
+              DO i = 1, klon
-                 pi_sulfate(i,k) = pi_so4(i,k,im2)   &
+                 sulfate(i, k) = so4(i, k, im2) &
-                      - FLOAT(iday-day2)/FLOAT(day1-day2) &
+                      - REAL(iday - day2)/REAL(day1 - day2) &
-                      * (pi_so4(i,k,im2) - pi_so4(i,k,im))
+                      * (so4(i, k, im2) - so4(i, k, im))
-                 IF (pi_sulfate(i,k).LT.0.) THEN
+                 IF (sulfate(i, k) < 0.) THEN
-                    IF (iday-day2.LT.0.) write(*,*) 'iday-day2',iday-day2
+                    IF (iday - day2 < 0.) write(*, *) 'iday - day2', iday - day2
-                    IF (pi_so4(i,k,im2) - pi_so4(i,k,im).LT.0.) &
+                    IF (so4(i, k, im2) - so4(i, k, im) < 0.) &
-                         write(*,*) 'pi_so4(i,k,im2) - pi_so4(i,k,im)', &
+                         write(*, *) 'so4(i, k, im2) - so4(i, k, im)', &
-                         pi_so4(i,k,im2) - pi_so4(i,k,im)
+                         so4(i, k, im2) - so4(i, k, im)
-                    IF (day1-day2.LT.0.) write(*,*) 'day1-day2',day1-day2
+                    IF (day1 - day2 < 0.) write(*, *) 'day1 - day2', day1 - day2
-                    stop 'pi_sulfate'
+                    stop 1
                  endif
               ENDDO
            ENDDO
         ENDIF
+        DO k = 1, klev
-        !JLD      ! The sulfate concentration [molec cm-3] is read in.
+           DO i = 1, klon
-        !JLD      ! Convert it into mass [ug SO4/m3]
+              so4_out(i, k) = sulfate(i, k)
-        !JLD      ! masse_so4 in [g/mol], n_avogadro in [molec/mol]
-        DO k=1,klev
-           DO i=1,klon
-              !JLD            pi_sulfate(i,k) = pi_sulfate(i,k)*masse_so4
-              !JLD     .           /n_avogadro*1.e+12
-              pi_so4_out(i,k) = pi_sulfate(i,k)
            ENDDO
         ENDDO
+     ELSE
-     ELSE ! If no new day, use old data:
+        ! If no new day, use old data:
-        DO k=1,klev
+        DO k = 1, klev
-           DO i=1,klon
+           DO i = 1, klon
-              pi_sulfate(i,k) = pi_so4_out(i,k)
+              sulfate(i, k) = so4_out(i, k)
            ENDDO
         ENDDO
-     ENDIF ! Was this the beginning of a new day?
+     ENDIF test_newday
    END SUBROUTINE readsulfate_preind

 Legend:



Removed from v.190
 


changed lines


 
Added in v.191
 Legend:



Removed from v.190
 


changed lines


 
Added in v.191
-Removed from v.190
+Added in v.191

	ViewVC Help
Powered by ViewVC 1.1.21