cvsp_p.f

Go to the documentation of this file.
00001 C:\opentelemac\v7p0\sources\sisyphe\cvsp_p.f
00002 !
00078                          SUBROUTINE CVSP_P
00079 !                        *****************
00080 !
00081      &(PATH_PRE,FILE_PRE,JG)
00082 !
00083 !***********************************************************************
00084 ! SISYPHE   V7P0                                   14/03/2013
00085 !***********************************************************************
00086 !
00087 !
00088 !
00089 !
00090 !
00091 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00092 !| JG             |<--| GLOBAL POINT NUMBER
00093 !| PATH_PRE       |<--| WHERE TO SAVE
00094 !| FILE_PRE       |<--| FILENAMETRUNK
00095 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00096 !
00097       USE DECLARATIONS_SISYPHE
00098       USE BIEF
00099       USE BIEF_DEF
00100       USE CVSP_OUTPUTFILES
00101 !
00102       IMPLICIT NONE
00103       INTEGER LNG,LU
00104       COMMON/INFO/LNG,LU
00105 !
00106 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00107 !
00108       INTEGER     , INTENT(IN)  :: JG
00109       CHARACTER(*), INTENT(IN) :: PATH_PRE
00110       CHARACTER(*), INTENT(IN) :: FILE_PRE
00111 !
00112 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00113 !
00114       CHARACTER(LEN=100) DEBUGFILE
00115       CHARACTER(LEN=5) OCSTR
00116       INTEGER I, K, J
00117       DOUBLE PRECISION AT, BSUM
00118 !
00119 !-----------------------------------------------------------------------
00120 !
00121       AT = DT*LT/PERCOU
00122       J = JG
00123       OUTPUTCOUNTER = OUTPUTCOUNTER + 1
00124 !
00125 !     GLOBAL NUMBERS TO GLOBAL NUMBERS
00126 !
00127       IF(NCSIZE.GT.1) THEN
00128         J = GLOBAL_TO_LOCAL_POINT(JG,MESH)
00129       ENDIF
00130 !
00131       WRITE(UNIT=OCSTR, FMT='(I5)') OUTPUTCOUNTER
00132       DO I=1,5
00133         IF(OCSTR(I:I)==' ') OCSTR(I:I)='0'
00134       ENDDO
00135 !
00136       WRITE(UNIT=DEBUGFILE, FMT='(A,A,A,A,I8,A,G15.8,A)')
00137      &     PATH_PRE,OCSTR,'_',FILE_PRE,
00138      &     JG,'_T_',AT,'.VSP.CSV'
00139 
00140       DO I=1,LEN_TRIM(DEBUGFILE)
00141         IF(DEBUGFILE(I:I)==' ') DEBUGFILE(I:I)='_'
00142       ENDDO
00143 !
00144       IF(J > 0) THEN
00145         OPEN(81, FILE=DEBUGFILE, STATUS='UNKNOWN' )
00146         REWIND 81
00147         WRITE(81,*)
00148      &       "J K FD50(I) AT PRO_D(K_I) PRO_F(K_I) X Y D50 ALT T H"
00149 !
00150         DO K=1,PRO_MAX(J)
00151           BSUM = 0.D0
00152           DO I=1,NSICLA
00153             BSUM = FDM(I)*PRO_F(J,PRO_MAX(J)+1-K,I) + BSUM
00154           ENDDO
00155 !
00156           DO I=1,NSICLA
00157             IF(K.EQ.1) THEN
00158 ! FULL OUTPUT WITH COORDINATES ETC. ON SURFACE
00159               WRITE (81,'(I8,1X,I4,1X,10(G20.12,1X))')
00160      &              JG,PRO_MAX(J)+1-K,FDM(I),AT,
00161      &              PRO_D(J,PRO_MAX(J)+1-K,I),
00162      &              PRO_F(J,PRO_MAX(J)+1-K,I),X(J),Y(J),
00163      &              BSUM,ES(J,1),TOB%R(J), Z%R(J)
00164             ELSE
00165 ! FOLLOWING SECTIONS
00166               WRITE (81,'(I8,1X,I4,1X,5(G20.12,1X))')
00167      &              JG,PRO_MAX(J)+1-K,FDM(I),AT,
00168      &              PRO_D(J,PRO_MAX(J)+1-K,I),
00169      &              PRO_F(J,PRO_MAX(J)+1-K,I)
00170             ENDIF
00171           ENDDO
00172         ENDDO
00173 !
00174         BSUM = 0.D0
00175         DO I=1,NSICLA
00176           BSUM = FDM(I)*PRO_F(J,1,I) + BSUM
00177         ENDDO
00178 !
00179         CLOSE(81)
00180 !
00181       ENDIF
00182 !
00183 !-----------------------------------------------------------------------
00184 !
00185       RETURN
00186       END SUBROUTINE CVSP_P

Generated on Fri Aug 31 2013 18:12:58 by S.E.Bourban (HRW) using doxygen 1.7.0