source: trunk/SOURCES/firstoutput-0.2.f @ 25

!> \file firstoutput-0.2.f
!! TOOOOOOOO DOOOOOOOO
!<

!> SUBROUTINE: FIRSTOUTPUT
!! TOOOOOOOO DOOOOOOOO
!! \author ... 
!! \date ...
!! @note Used modules:
!! @note    - use module3D_phy
!! @note    - use ISO_DECLAR 
!! @note    - use DEFORM_DECLAR
!<     
      subroutine FIRSTOUTPUT()

c     ************************************************
c     * LEGENDE DE LA SORTIE ECRAN + 1ere SORTIE  ****
c     ************************************************

      USE module3D_phy
      USE ISO_DECLAR   ! pour nbed
      USE DEFORM_DECLAR

      implicit none

c     ** for line display **
      write(num_ritz,*) 'RunName  ',runname
      write(num_ritz,fmt='(a8,i1)') 'reprise=', icompteur
      if (icompteur.eq.1) then
         write(num_ritz,*) 'fichier sortie pour reprise=', iout
         write(num_ritz,*) 'fichier reprise=', reprcptr
      endif
      write(num_ritz,*)
!      write(num_ritz,fmt='(a9,i1,2x,a9,i1)')
      write(num_ritz,*) 'ICOUPLE=', ICOUPLE,'    GEOPLACE=',GEOPLACE



      write(num_ritz,*) 'CSI=',CSI,'  Csnow=',Csnow,'  Cice=',Cice
!     write(num_ritz,*) 'deltaacc=',deltaacc
      write(num_ritz,*) 'DZM=',DZM,'  ROM=',ROM,'  CL=',CL
      write(num_ritz,*) 'CM=',CM,'  CPM=',CPM,'  NCOND=',NCOND
      write(num_ritz,*) 'SF=',SF,'  BAT1=',BAT1,'  Q1=',Q1
      write(num_ritz,*) 'BAT2=',BAT2,'  Q2=',Q2
      write(num_ritz,*) 'RGAS=',RGAS,'  RO=',RO,'  ROW=',ROW,'  DA=',DA
      write(num_ritz,*)
      write(num_ritz,*) '*******************************'
      write(num_ritz,*) 'shelf et glissement'
      write(num_ritz,*) '*******************************'
      write(num_ritz,*) ' cf=',cf,'      toblim=',toblim/1.e5,'bar'
      write(num_ritz,*) ' Neffgz=',neffgz/rog,'m au dessus flotaison'
      write(num_ritz,*) ' moteurmax=',moteurmax/1.e5,'    betamax=',
     & betamax/1.e5,'bar'
      write(num_ritz,*) 'fusion ice shelf : voir le module'
      write(num_ritz,*) '*******************************'
      write(num_ritz,*)
      write(num_ritz,*) 'HWATERMAX=',HWATERMAX,'   COEFBMAX=',COEFBMAX
      write(num_ritz,*) 'NICE=',NICE,'  NBED=',NBED
      write(num_ritz,*)

        write(num_ritz,902) '   nt    TIME   TAFOR  sealevel    vol '//
     &  '  NP  <H>  HMAX  '// 
     &  '  <BM>   <TB>   <VS>'//
     &  '   <TBdot>   <VSdot>   <Hdot>   dt   <BMELT>'//
     &  '  <ACC>'

      write(num_ritz,*)
      write(num_ritz,*)
c     (****** ICE SHEET STATISTICS and SHORT OUTPUT *******)
c     if (NTRACE.eq.1) write (6,*)'%%%% ICE SHEET STATISTICS %%%%'
      NP=0
      VOL=0.
      HMAX=0.  
      do I=1,NX
      do J=1,NY 
        if (H(I,J).gt.0.) then 
          NP=NP+1 
          VOL=VOL+H(I,J) 
          if (H(I,J).gt.HMAX)  HMAX=H(I,J) 
        endif 
      end do
      end do
      if (NP.ne.0) then 
        HMEAN=vol/NP 
        VOL=VOL*DX*DY 
        SURF=NP*DX*DY 
      endif

c     **** short display ****
      write(num_ritz,903) TIME,TAFOR,VOL,NP,int(HMEAN)

902   format (A130)
903   format(f9.1,1x,f5.1,1x,e10.5,1x,i6,1x,i4,1x,i5,1X,
     &  f7.3,1x,f6.2,1X,i3,4(1x,e7.1),1X,f6.2,1x,
     &  f6.2,1X,i5,1X,f6.2,1X,i5,1X,f7.3) 

      end