Changeset 1301 for trunk/LMDZ.VENUS/libf/phyvenus/load_ksi.F

Timestamp:

Jun 26, 2014, 12:25:31 PM (10 years ago)

Author:

slebonnois

Message:

SL: many bug corrections in phyvenus, some cleaning, and a new ksi matrix format for Venus IR

File:

• trunk/LMDZ.VENUS/libf/phyvenus/load_ksi.F (modified) (4 diffs)

trunk/LMDZ.VENUS/libf/phyvenus/load_ksi.F

@@ -23 +23 @@
 C     MODIFICATIONS.
 C     --------------
-C        version multimatrice (topographie, sommet nuages): 20/12/2006
+C        
+c   New ksi matrix: possibility of different cloud model fct of lat   05/2014
 C     ------------------------------------------------------------------
 C
@@ -30 +31 @@
 c inputs
       real   psurf(klon)           ! Surface pressure
-      real   ztop(klon)            ! Altitude of the top of cloud deck (km)
 c outputs
       real   ksive(0:klev+1,0:klev+1,nnuve,nbmat)  ! ksi matrixes in Vincent's file
 
 c local variables
-      integer i,j,ig,band,pve,nlve
-      integer mat,Nb,m,Nmat,nl_init,mat0
-      parameter(nl_init=8)
+      integer i,j,isza,ips,band,pve,sve,nlve
+      integer lat,Nb,m,mat
       character*9 tmp1
       character*100 file
+      CHARACTER*2 str2
       real   lambda(nnuve)            ! wavelenght in table (mu->m, middle of interval)
       real   lambdamin(nnuve),lambdamax(nnuve) ! in microns
-      real   dlambda                  ! cm-1  
+      real   dlambda                  ! in microns 
 
       nlve = klev
@@ -53 +53 @@
       open(10,file=file)
       
-      do i=1,nl_init-1
-         read(10,*)
-      enddo
-      read(10,*) (tmp1,i=1,4),Nmat
+      read(10,*)
+      read(10,*) nlatve
+      read(10,*)
 
-      if (nbmat.ne.Nmat) then
-         write(*,*) 'This is subroutine load_ksi'
-         print*,'Probleme de dimension entre ksi.txt et le param nbmat'
-         print*,'Nb matrices = ',nbmat,Nmat
-         stop
-      endif
-
-      do mat=1,nbmat
+      write(*,*) 'This is subroutine load_ksi'
+      write(*,*) 'Nb of lat bands:',nlatve
+      
+      do lat=1,nlatve
+        read(10,*) !line for lat range
+        read(10,*) indexve(lat)
+        read(10,*) nbpsve(lat)
         read(10,*)
-        read(10,*) 
+        read(10,*) nbszave(lat)
+        read(10,*)
+        
+        do isza=1,nbszave(lat)
+          do ips=1,nbpsve(lat)
+          
+        read(10,*) (tmp1,j=1,3),mat    !line for matrix number
         read(10,*) (tmp1,j=1,2),pve
-        psurfve(mat) = pve*1.e5  ! pve en bar, psurfve en Pa
-        read(10,*) (tmp1,j=1,7),ztopve(mat)
-        ztopve(mat) = ztopve(mat)*1.e-3 ! passage en km
+        psurfve(ips,lat) = pve*1.e5    ! pve in bar, psurfve in Pa
+        read(10,*) (tmp1,j=1,3),sve 
+        szave(isza,lat) = cos(sve*RPI/180.) ! conversion in mu0
         read(10,*)
         read(10,*) m,Nb
         if (m.ne.nlve) then
          write(*,*) 'This is subroutine load_ksi'
-         print*,'Probleme de dimension entre ksi.txt et le param nlve'
+         print*,'Dimension problem between ksi.txt and nlve'
          print*,'N levels = ',m,nlve
          stop
@@ -82 +86 @@
         if (Nb.ne.nnuve) then
          write(*,*) 'This is subroutine load_ksi'
-         print*,'Probleme de dimension entre ksi.txt et le param nnuve'
+         print*,'Dimension problem between ksi.txt and nnuve'
          print*,'N freq = ',Nb,nnuve
          stop
         endif
 c     Now reading ksi matrix index "mat"
+        write(str2,'(i2.2)') m+2
         do band=1,Nb
          read(10,*) lambdamin(band),lambdamax(band)
          do i=0,m+1
-            read(10,'(100e17.9)') (ksive(i,j,band,mat),j=0,m+1) ! sr/µm/cm¯¹
+            read(10,'('//str2//'e17.9)') (ksive(i,j,band,mat),j=0,m+1) ! no unit
          enddo                  ! i
         enddo                     ! band
 c       print*,"Matrice ",mat," lue"
-c       print*,"   psurf=",psurfve(mat)," bars, Ztop=",ztopve(mat)," km"
-      enddo  ! mat
+c       print*,"   psurf=",psurfve(ips,lat)," bars"
+        if (mat+1.gt.nbmat) then
+         write(*,*) 'This is subroutine load_ksi'
+         print*,'Dimension problem between ksi.txt and nbmat'
+         print*,'(max number of matrixes)'
+         print*,'nbmat (in comcstVE.h) should be raised'
+         stop
+        endif
+
+          enddo    ! ips
+        enddo      ! isza
+      enddo        ! lat
+      
+      write(*,*) 'Total nb of matrixes:',mat
       
       close(10)
 
-c longueur d'onde centrale et largeur de chaque bande
+c central wavelength and wavelength width
       do band=1,nnuve
-         lambda(band)=(lambdamin(band)+lambdamax(band))/2.*1.e-6   ! en m
-         dlambda     =(1./lambdamin(band)-1./lambdamax(band))*1.e4 ! en cm-1
+         lambda(band)=(lambdamin(band)+lambdamax(band))/2.*1.e-6   ! in m
+         dlambda     =(lambdamax(band)-lambdamin(band))            ! in microns
 c        print*,band,lambdamin(band),dlambda,lambdamax(band)
 
-c changement de convention (signe) pour ksi,
-c et prise en compte de la largeur de bande (en cm-1):
+c sign convention for ksi,
+c and taking into account the wavelength width (in microns):
          do mat=1,nbmat
          do i=0,nlve+1
            do j=0,nlve+1
-              ksive(i,j,band,mat) = -ksive(i,j,band,mat)*dlambda
+              ksive(i,j,band,mat) = +ksive(i,j,band,mat)*dlambda    ! in µm
            enddo
          enddo
          enddo
-c calcul des coeff al et bl pour luminance Planck 
+c computing coeff al and bl for Planck luminance 
          al(band) = 2.*RHPLA*RCLUM*RCLUM/(lambda(band))**5.
-c cette luminance doit etre en W/m²/sr/µm pour correspondre au calcul
-c des ksi. Ici, elle est en W/m²/sr/m donc il faut mettre un facteur 1.e-6
+c in W/m²/m
+c We need W/m²/µm :
      .                * 1.e-6
          bl(band) = RHPLA*RCLUM/(RKBOL*lambda(band))