Changeset 658 for trunk/LMDZ.MARS/libf/aeronomars/param_read.F

Timestamp:

May 12, 2012, 8:10:08 AM (13 years ago)

Author:

emillour

Message:

Mars GCM:

• updated high atmosphere photochemistry (jthermcalc.F, param_v4.h, iono.h, paramfoto_compact.F, param_read.F , thermosphere.F).
• minor change in calchim.F90 (to not use maxloc(zycol, dim = 2) function which seems to be a problem for g95) .
• minor bug fix in perosat.F; set tendency on pdqscloud for h2o2 to zero if none is computed.
• in "moldiff.F", changed "tridag" to "tridag_sp", "LUBKSB" to "LUBKSB_SP" and "LUDCMP" to "LUDCMP_SP" to avoid possible conflicts with same routines defined in "moldiff_red.F". Added use of automatic-sized array in "tridag" and "tridag_sp" local array "gam" (to avoid using an a priori oversized local array).

FGG+JYC+EM

File:

• trunk/LMDZ.MARS/libf/aeronomars/param_read.F (modified) (2 diffs)

trunk/LMDZ.MARS/libf/aeronomars/param_read.F

@@ -95 +95 @@
       IF (ierr.NE.0) THEN 
        write(*,*)'cant find directory EUVDAT and content coln.dat'
-       write(*,*)'(in aeronomars/param_read_iono.F)'
+       write(*,*)'(in aeronomars/param_read.F)'
        write(*,*)'It should be in :', datafile,'/'
        write(*,*)'1) You can change this directory address in '
@@ -181 +181 @@
 
       do i=nz2,1,-1
-         read(220,*) (jabsifotsintpar(j,2,i),j=1,16)
-      end do
-      
-      do i=nz2,1,-1
-         read(230,*) (jabsifotsintpar(j,3,i),j=1,16)
-      end do
-
-      do i=nz2,1,-1
-         read(240,*) (jabsifotsintpar(j,1,i),j=1,16)
-      end do
-
-      do i=nz2,1,-1
-         read(250,*) (jabsifotsintpar(j,2,i),j=17,24)
-      end do
-
-
-      do i=nz2,1,-1
-         read(260,*) (jabsifotsintpar(j,2,i),j=25,31)
-      end do
-
-      do i=nz2,1,-1
-         read(270,*) (jabsifotsintpar(j,1,i),j=17,24)
-      end do
-
-      do i=nz2,1,-1
-         read(280,*) (jabsifotsintpar(j,1,i),j=25,31)
-      end do
-
-      do i=nz2,1,-1
-         read(290,*) jabsifotsintpar(32,1,i)
-      end do
-
-      do i=nz2,1,-1
-         read(300,*) (jabsifotsintpar(j,2,i),j=32,34)
-      end do
-
-      do i=nz2,1,-1
-         read(160,*) (jabsifotsintpar(j,5,i),j=1,15)
-      end do
-
-      do i=nz2,1,-1
-         read(150,*) (jabsifotsintpar(j,4,i),j=25,31)
-      end do
-
-      do i=nz2,1,-1
-         read(170,*) (jabsifotsintpar(j,6,i),j=25,35)
-      end do
-
-      do i=nz2,1,-1
-         read(180,*) (jabsifotsintpar(j,7,i),j=34,36)
-      end do
-
-      do i=nz2,1,-1
-         read(390,*) (jabsifotsintpar(j,8,i),j=2,16)
-      enddo
-
-      do i=nz2,1,-1
-         read(400,*) (jabsifotsintpar(j,8,i),j=17,24)
-      enddo
-
-      do i=nz2,1,-1
-         read(410,*) (jabsifotsintpar(j,9,i),j=1,16)
-      enddo
-
-      do i=nz2,1,-1
-         read(420,*) (jabsifotsintpar(j,10,i),j=2,16)
-      enddo
-
-      do i=nz2,1,-1
-         read(430,*) (jabsifotsintpar(j,10,i),j=17,24)
-      enddo
-
-      do i=nz2,1,-1
-         read(440,*) (jabsifotsintpar(j,10,i),j=25,32)
-      enddo
-
-      do i=nz2,1,-1
-         read(450,*) (jabsifotsintpar(j,11,i),j=2,16)
-      enddo
-
-      do i=nz2,1,-1
-         read(460,*) (jabsifotsintpar(j,11,i),j=17,24)
-      enddo
-
-      do i=nz2,1,-1
-         read(470,*) (jabsifotsintpar(j,11,i),j=25,29)
-      enddo
-      
-      do i=nz2,1,-1
-         read(480,*) (jabsifotsintpar(j,12,i),j=2,16)
-      enddo
-
-      do i=nz2,1,-1
-         read(490,*) (jabsifotsintpar(j,13,i),j=2,16)
-      enddo
-      
-      do i=nz2,1,-1
-         read(500,*) (jabsifotsintpar(j,13,i),j=17,24)
-      enddo
-      
-      do i=nz2,1,-1
-         read(510,*) (jabsifotsintpar(j,13,i),j=25,36)
+         read(220,*) (jabsifotsintpar(i,2,j),j=1,16)
+      end do
+      
+      do i=nz2,1,-1
+         read(230,*) (jabsifotsintpar(i,3,j),j=1,16)
+      end do
+
+      do i=nz2,1,-1
+         read(240,*) (jabsifotsintpar(i,1,j),j=1,16)
+      end do
+
+      do i=nz2,1,-1
+         read(250,*) (jabsifotsintpar(i,2,j),j=17,24)
+      end do
+
+
+      do i=nz2,1,-1
+         read(260,*) (jabsifotsintpar(i,2,j),j=25,31)
+      end do
+
+      do i=nz2,1,-1
+         read(270,*) (jabsifotsintpar(i,1,j),j=17,24)
+      end do
+
+      do i=nz2,1,-1
+         read(280,*) (jabsifotsintpar(i,1,j),j=25,31)
+      end do
+
+      do i=nz2,1,-1
+         read(290,*) jabsifotsintpar(i,1,32)
+      end do
+
+      do i=nz2,1,-1
+         read(300,*) (jabsifotsintpar(i,2,j),j=32,34)
+      end do
+
+      do i=nz2,1,-1
+         read(160,*) (jabsifotsintpar(i,5,j),j=1,15)
+      end do
+
+      do i=nz2,1,-1
+         read(150,*) (jabsifotsintpar(i,4,j),j=25,31)
+      end do
+
+      do i=nz2,1,-1
+         read(170,*) (jabsifotsintpar(i,6,j),j=25,35)
+      end do
+
+      do i=nz2,1,-1
+         read(180,*) (jabsifotsintpar(i,7,j),j=34,36)
+      end do
+
+      do i=nz2,1,-1
+         read(390,*) (jabsifotsintpar(i,8,j),j=2,16)
+      enddo
+
+      do i=nz2,1,-1
+         read(400,*) (jabsifotsintpar(i,8,j),j=17,24)
+      enddo
+
+      do i=nz2,1,-1
+         read(410,*) (jabsifotsintpar(i,9,j),j=1,16)
+      enddo
+
+      do i=nz2,1,-1
+         read(420,*) (jabsifotsintpar(i,10,j),j=2,16)
+      enddo
+
+      do i=nz2,1,-1
+         read(430,*) (jabsifotsintpar(i,10,j),j=17,24)
+      enddo
+
+      do i=nz2,1,-1
+         read(440,*) (jabsifotsintpar(i,10,j),j=25,32)
+      enddo
+
+      do i=nz2,1,-1
+         read(450,*) (jabsifotsintpar(i,11,j),j=2,16)
+      enddo
+
+      do i=nz2,1,-1
+         read(460,*) (jabsifotsintpar(i,11,j),j=17,24)
+      enddo
+
+      do i=nz2,1,-1
+         read(470,*) (jabsifotsintpar(i,11,j),j=25,29)
+      enddo
+      
+      do i=nz2,1,-1
+         read(480,*) (jabsifotsintpar(i,12,j),j=2,16)
+      enddo
+
+      do i=nz2,1,-1
+         read(490,*) (jabsifotsintpar(i,13,j),j=2,16)
+      enddo
+      
+      do i=nz2,1,-1
+         read(500,*) (jabsifotsintpar(i,13,j),j=17,24)
+      enddo
+      
+      do i=nz2,1,-1
+         read(510,*) (jabsifotsintpar(i,13,j),j=25,36)
       enddo