Changeset 37 for branches/petra/src/getfields.f90

Timestamp:

Apr 22, 2015, 3:42:35 PM (9 years ago)

Author:

pesei

Message:

Wet dep quick fix and other small changes. Wet depo quick fix not final yet.

File:

• branches/petra/src/getfields.f90 (modified) (5 diffs)

branches/petra/src/getfields.f90

@@ -41 +41 @@
   !   Variables tth,qvh,tthn,qvhn (on eta coordinates) in common block.
   !   Function of nstop extended.
+  !
+  !  3/2015, PS: put in verbosity output
+  !
   !*****************************************************************************
   !                                                                            *
@@ -83 +86 @@
   integer :: indmin = 1
 
+  call verboutput(verbosity,1,'  getfields')
 
   ! Check, if wind fields are available for the current time step
@@ -125 +129 @@
     do indj=indmin,numbwf-1
        if (ldirect*wftime(indj+1).gt.ldirect*itime) then
+           call verboutput(verbosity,1,'    call readwind')
           call readwind(indj+1,memind(2),uuh,vvh,wwh)
+           call verboutput(verbosity,1,'    call readwind_nests')
           call readwind_nests(indj+1,memind(2),uuhn,vvhn,wwhn)
+           call verboutput(verbosity,1,'    call calcpar')
           call calcpar(memind(2),uuh,vvh,pvh)
+           call verboutput(verbosity,1,'    call calcpar_nests')
           call calcpar_nests(memind(2),uuhn,vvhn,pvhn)
+           call verboutput(verbosity,1,'    call verttransform')
           call verttransform(memind(2),uuh,vvh,wwh,pvh)
+           call verboutput(verbosity,1,'    call verttransform_nests')
           call verttransform_nests(memind(2),uuhn,vvhn,wwhn,pvhn)
           memtime(2)=wftime(indj+1)
@@ -148 +158 @@
              (ldirect*wftime(indj+1).gt.ldirect*itime)) then
            memind(1)=1
+            call verboutput(verbosity,1,'    call readwind')
            call readwind(indj,memind(1),uuh,vvh,wwh)
+            call verboutput(verbosity,1,'    call readwind_nests')
            call readwind_nests(indj,memind(1),uuhn,vvhn,wwhn)
+            call verboutput(verbosity,1,'    call calcpar')
            call calcpar(memind(1),uuh,vvh,pvh)
+            call verboutput(verbosity,1,'    call calcpar_nests')
            call calcpar_nests(memind(1),uuhn,vvhn,pvhn)
+            call verboutput(verbosity,1,'    call verttransform')
            call verttransform(memind(1),uuh,vvh,wwh,pvh)
+            call verboutput(verbosity,1,'    call verttransform_nests')
            call verttransform_nests(memind(1),uuhn,vvhn,wwhn,pvhn)
            memtime(1)=wftime(indj)
            memind(2)=2
+            call verboutput(verbosity,1,'    call readwind')
            call readwind(indj+1,memind(2),uuh,vvh,wwh)
+            call verboutput(verbosity,1,'    call readwind_nests')
            call readwind_nests(indj+1,memind(2),uuhn,vvhn,wwhn)
+            call verboutput(verbosity,1,'    call calcpar')
            call calcpar(memind(2),uuh,vvh,pvh)
+            call verboutput(verbosity,1,'    call calcpar_nests')
            call calcpar_nests(memind(2),uuhn,vvhn,pvhn)
+            call verboutput(verbosity,1,'    call verttransform')
            call verttransform(memind(2),uuh,vvh,wwh,pvh)
+            call verboutput(verbosity,1,'    call verttransform_nests')
            call verttransform_nests(memind(2),uuhn,vvhn,wwhn,pvhn)
            memtime(2)=wftime(indj+1)
@@ -175 +197 @@
   if (lwindinterv.gt.idiffmax) nstop=3
 
+  call verboutput(verbosity,1,'    leaving getfields')
+
 end subroutine getfields