Changeset 5b5589b in flex_extract.git for Source

Timestamp:

Aug 21, 2019, 8:11:08 PM (5 years ago)

Author:

pesei <petra seibert -@…>

Branches:

master, ctbto, dev

Children:

7e59510, e0005c9

Parents:

dfa7dbd

git-author:

pesei <petra seibert -@…> (08/21/19 20:06:51)

git-committer:

pesei <petra seibert -@…> (08/21/19 20:11:08)

Message:

Update automated Fortran documentation

and move everything to the right place.
Replace tabs with blanks in two source files (so that FORD works!)

Location:

Source/Fortran

Files:

• phgrreal.f90 (modified) (7 diffs)
• rwgrib2.f90 (modified) (2 diffs)

Source/Fortran/phgrreal.f90

@@ -14 +14 @@
 !          CX00,CX01,CX11,CX02,....CXMNAUFMNAUF
 ! FELD   = FELD DER METEOROLOGISCHEN VARIABLEN
-!       WSAVE  = Working Array fuer Fouriertransformation
-! Z         = LEGENDREFUNKTIONSWERTE
+! WSAVE  = Working Array fuer Fouriertransformation
+! Z      = LEGENDREFUNKTIONSWERTE
 !
 ! MNAUF  ANZAHL DER FOURIERKOEFFIZIENTEN
@@ -24 +24 @@
     IMPLICIT NONE
 
-!                       Anzahl der Gitterpunkte auf jedem Breitenkreis
+!   Anzahl der Gitterpunkte auf jedem Breitenkreis
     INTEGER MLAT(MAXB/2)
     INTEGER K,MAXL,MAXB,MLEVEL,MNAUF
@@ -84 +84 @@
         LLPS=LLP
         IF (2*I+1 .LT. MLAT(K)) THEN
-!               Innerste Schleife aufgespalten um if-Abfrage zu sparen
+!         Innerste Schleife aufgespalten um if-Abfrage zu sparen
           DO 18 J=I,MNAUF,2
             SCR=SCR+Z(LLP,K)*CXMN(2*LL,L)
@@ -124 +124 @@
 !! FELDVARIABLEN VOM PHASENRAUM IN  DEN PHYSIKALISCHEN
 !! RAUM AUF KUGELKOORDINATEN DURCH. Es kann ein Teilausschnitt
-!!      Der Erde angegeben werden. Diese Routine ist langsamer als phgrph
+!! Der Erde angegeben werden. Diese Routine ist langsamer als phgrph
 
 ! CXMN   = SPEKTRALKOEFFIZIENTEN IN DER REIHENFOLGE
@@ -168 +168 @@
       DO 11 J=5,5
         IF (ABS(ABS(Z(100,J))-ABS(Z(100,MAXB+1-J))) .GT. 1E-11) SYM=.FALSE.
-!             WRITE(*,*) ABS(Z(100,J)),ABS(Z(100,MAXB+1-J))
+!       WRITE(*,*) ABS(Z(100,J)),ABS(Z(100,MAXB+1-J))
 11    CONTINUE
 !!      WRITE(*,*) 'Symmetrisch: ',SYM
@@ -225 +225 @@
         LLS=LL
         LLPS=LLP
-!             Innerste Schleife aufgespalten um if-Abfrage zu sparen
+!       Innerste Schleife aufgespalten um if-Abfrage zu sparen
         DO 18 K=I,MNAUF,2
           SCR=SCR+Z(LLP,J)*CXMN(2*LL,L)
@@ -398 +398 @@
 ! Die Polynome sind wie im ECMWF indiziert, d.h.
 ! P00,P10,P11,P20,P21,P22,...
-!       Ansonsten ist die Routine analog zu PLGNDN
+! Ansonsten ist die Routine analog zu PLGNDN
 ! X IST DER COSINUS DES ZENITWINKELS ODER
 !       DER SINUS DER GEOGRAFISCHEN BREITE

Source/Fortran/rwgrib2.f90

@@ -131 +131 @@
      do i=1,LLEN
        if (MLEVELIST(i:i) .eq. '/') THEN
-                l=l+1
-                MLINDEX(l)=i
-              end if
-            end do
-            MLINDEX(l+1)=LLEN+1
-            do i=1,l
-              read(MLEVELIST(MLINDEX(i)+1:MLINDEX(i+1)-1),*) ILEVEL(i)
-            end do
+         l=l+1
+         MLINDEX(l)=i
+       end if
+     end do
+     MLINDEX(l+1)=LLEN+1
+     do i=1,l
+       read(MLEVELIST(MLINDEX(i)+1:MLINDEX(i+1)-1),*) ILEVEL(i)
+     end do
    end if 
 
@@ -153 +153 @@
 !           call grib_set(igrib,"units","s**-1")
 !         end if
-              if (l .ne. mlevel) then
+       if (l .ne. mlevel) then
          zsec4(1:maxl*maxb)=RESHAPE(FELD(:,:,ILEVEL(k)),(/maxl*maxb/))
-              else
-                zsec4(1:maxl*maxb)=RESHAPE(FELD(:,:,k),(/maxl*maxb/))
-              end if
+       else
+         zsec4(1:maxl*maxb)=RESHAPE(FELD(:,:,k),(/maxl*maxb/))
+       end if
        call grib_set(igrib,"values",zsec4)