| 1 | PROGRAM sla2fb |
|---|
| 2 | !!--------------------------------------------------------------------- |
|---|
| 3 | !! |
|---|
| 4 | !! ** PROGRAM sla2fb ** |
|---|
| 5 | !! |
|---|
| 6 | !! ** Purpose : Convert AVISO SLA format to feedback format |
|---|
| 7 | !! |
|---|
| 8 | !! ** Method : Use of utilities from obs_fbm. |
|---|
| 9 | !! |
|---|
| 10 | !! ** Action : |
|---|
| 11 | !! |
|---|
| 12 | !! Usage: |
|---|
| 13 | !! sla2fb.exe [-s type] outputfile inputfile1 inputfile2 ... |
|---|
| 14 | !! Option: |
|---|
| 15 | !! -s Select altimeter data_source |
|---|
| 16 | !! |
|---|
| 17 | !! History : |
|---|
| 18 | !! ! 2010 (K. Mogensen) Initial version |
|---|
| 19 | !!---------------------------------------------------------------------- |
|---|
| 20 | USE obs_fbm |
|---|
| 21 | USE obs_sla_io |
|---|
| 22 | USE convmerge |
|---|
| 23 | IMPLICIT NONE |
|---|
| 24 | ! |
|---|
| 25 | ! Command line arguments for output file and input files |
|---|
| 26 | ! |
|---|
| 27 | #ifndef NOIARGCPROTO |
|---|
| 28 | INTEGER,EXTERNAL :: iargc |
|---|
| 29 | #endif |
|---|
| 30 | INTEGER :: nargs |
|---|
| 31 | CHARACTER(len=256) :: cdoutfile |
|---|
| 32 | CHARACTER(len=256),ALLOCATABLE :: cdinfile(:) |
|---|
| 33 | CHARACTER(len=256) :: cdtmp |
|---|
| 34 | CHARACTER(len=5) :: cdsource |
|---|
| 35 | ! |
|---|
| 36 | ! Input data |
|---|
| 37 | ! |
|---|
| 38 | TYPE(obfbdata), POINTER :: slaf(:) |
|---|
| 39 | INTEGER :: ninfiles,ntotobs |
|---|
| 40 | ! |
|---|
| 41 | ! Output data |
|---|
| 42 | ! |
|---|
| 43 | TYPE(obfbdata) :: fbdata |
|---|
| 44 | ! |
|---|
| 45 | ! Loop variables |
|---|
| 46 | ! |
|---|
| 47 | INTEGER :: ip,ia,ji,jk,noff |
|---|
| 48 | ! |
|---|
| 49 | ! Get number of command line arguments |
|---|
| 50 | ! |
|---|
| 51 | nargs=IARGC() |
|---|
| 52 | IF (nargs < 1) THEN |
|---|
| 53 | WRITE(*,'(A)')'Usage:' |
|---|
| 54 | WRITE(*,'(A)')'sla2fb [-s type] outputfile inputfile1 inputfile2 ...' |
|---|
| 55 | CALL abort() |
|---|
| 56 | ENDIF |
|---|
| 57 | cdsource='' |
|---|
| 58 | ! |
|---|
| 59 | ! Get input data |
|---|
| 60 | ! |
|---|
| 61 | noff=1 |
|---|
| 62 | IF ( nargs > 1 ) THEN |
|---|
| 63 | CALL getarg(1,cdtmp) |
|---|
| 64 | IF (TRIM(cdtmp)=='-s') THEN |
|---|
| 65 | IF ( nargs < 3 ) THEN |
|---|
| 66 | WRITE(*,*)'Missing arguments to -s <datasource>' |
|---|
| 67 | CALL abort |
|---|
| 68 | ENDIF |
|---|
| 69 | CALL getarg(2,cdsource) |
|---|
| 70 | noff=3 |
|---|
| 71 | ENDIF |
|---|
| 72 | ENDIF |
|---|
| 73 | CALL getarg(noff,cdoutfile) |
|---|
| 74 | ninfiles = nargs - noff |
|---|
| 75 | ALLOCATE( slaf(MAX(nargs-noff,1)) ) |
|---|
| 76 | ALLOCATE( cdinfile(nargs-noff) ) |
|---|
| 77 | ntotobs = 0 |
|---|
| 78 | DO ia=1,ninfiles |
|---|
| 79 | CALL getarg( ia + noff, cdinfile(ia) ) |
|---|
| 80 | WRITE(*,'(2A)')'File = ',TRIM(cdinfile(ia)) |
|---|
| 81 | CALL read_avisofile( TRIM(cdinfile(ia)), slaf(ia), 6, .TRUE., .FALSE. ) |
|---|
| 82 | WRITE(*,'(A,I9,A)')'has',slaf(ia)%nobs,' observations' |
|---|
| 83 | IF (LEN_TRIM(cdsource)>0) THEN |
|---|
| 84 | DO ji=1,slaf(ia)%nobs |
|---|
| 85 | slaf(ia)%cdwmo(ji)=TRIM(slaf(ia)%cdwmo(ji))//'_'//TRIM(cdsource) |
|---|
| 86 | ENDDO |
|---|
| 87 | ENDIF |
|---|
| 88 | ntotobs = ntotobs + slaf(ia)%nobs |
|---|
| 89 | ENDDO |
|---|
| 90 | IF (ninfiles==0) THEN |
|---|
| 91 | CALL init_obfbdata( slaf(1) ) |
|---|
| 92 | CALL alloc_obfbdata( slaf(1), 1, 0, 1, 0, 0, .FALSE. ) |
|---|
| 93 | slaf(1)%cname(1) = 'SLA' |
|---|
| 94 | slaf(1)%cdjuldref = '19500101000000' |
|---|
| 95 | ENDIF |
|---|
| 96 | WRITE(*,'(A,I8)') 'Total observations : ',ntotobs |
|---|
| 97 | ! |
|---|
| 98 | ! Merge and output the data. |
|---|
| 99 | ! |
|---|
| 100 | CALL conv_fbmerge( TRIM(cdoutfile), ninfiles, slaf ) |
|---|
| 101 | |
|---|
| 102 | END PROGRAM sla2fb |
|---|
