Changeset 304

Timestamp:

Sep 22, 2011, 1:52:53 PM (13 years ago)

Author:

rwordsworth

Message:

Option added for FZERO but not implemented. May be
useful for optimization in future.

File:

• trunk/LMDZ.GENERIC/libf/phystd/sugas_corrk.F90 (modified) (4 diffs)

trunk/LMDZ.GENERIC/libf/phystd/sugas_corrk.F90

@@ -47 +47 @@
       integer ngas, igas
 
-      ! temporary special for addh2
-      double precision testH2
+      double precision testcont ! for continuum absorption initialisation
 
 !=======================================================================
@@ -120 +119 @@
             print*,'Name of a gas in radiative transfer data (',gastype(n),') does not', &
                  'match that in gases.def (',gnom(n),'), exiting.'
+            call abort
          endif
       enddo
@@ -294 +294 @@
          do nw=1,L_NSPECTI
             fzeroi(nw) = 0.
-         end do
+!            do nt=1,L_NTREF
+!               do np=1,L_NPREF
+!                  do nh=1,L_REFVAR
+!                     do ng = 1,L_NGAUSS
+!                        if(gasi8(nt,np,nh,nw,ng).lt.1.0e-25)then
+!                           fzeroi(nw)=fzeroi(nw)+1.
+!                        endif
+!                     end do
+!                  end do
+!               end do
+!            end do
+!            fzeroi(nw)=fzeroi(nw)/dble(L_NTREF*L_NPREF*L_REFVAR*L_NGAUSS)
+         end do
+
+         do nw=1,L_NSPECTV
+            fzerov(nw) = 0.
+!            do nt=1,L_NTREF
+!               do np=1,L_NPREF
+!                  do nh=1,L_REFVAR
+!                     do ng = 1,L_NGAUSS
+!                        if(gasv8(nt,np,nh,nw,ng).lt.1.0e-25)then
+!                           fzerov(nw)=fzerov(nw)+1.
+!                        endif
+!                     end do
+!                  end do
+!               end do
+!            end do
+!            fzerov(nw)=fzerov(nw)/dble(L_NTREF*L_NPREF*L_REFVAR*L_NGAUSS)
+         end do
+
+
          do nw=1,L_NSPECTV
             fzerov(nw) = 0.
@@ -473 +503 @@
       end do
 
+
+
       do igas=1,ngasmx
          if(gnom(igas).eq.'H2_')then 
-            call interpolateH2H2(500.D+0,250.D+0,17500.D+0,testH2,.true.)
+            call interpolateH2H2(500.D+0,250.D+0,17500.D+0,testcont,.true.)
+         elseif(gnom(igas).eq.'H2O')then 
+            call interpolateH2Ocont(990.D+0,296.D+0,683.2D+0*2,0.D+0,testcont,.true.)
          endif
       enddo
 
-
       return
     end subroutine sugas_corrk