differences important for GECKO-A
(ie NOT simply differences between versions)
=================================
4,5c4,5
< *=    Version 4.5                                                            =*
< *=    Sep 2007                                                               =*
---
> *=    Version 4.2                                                            =*
> *=    May 2003                                                               =*
111a111,119
> *cba following added for MCM or saprc99 scheme
>       INTEGER nj4,nj6,nj8
>       REAL sj4(kj,kz,kw),valj4(kj,kz)
>       REAL sj6(kj,kz,kw), valj6(kj,kz),sj8(kj,kz,kw), valj8(kj,kz)
>       CHARACTER*40 jlabel3(kj), jlabel4(kj), jlabel5(kj)
>       CHARACTER*40 jlabel6(kj), jlabel7(kj),jlabel8(kj), jlabel9(kj)
>       REAL resu4(10,kj), resu6(10,kj),resu8(10,kj)
>       REAL toto1, toto2
> 
403a406,412
> *cba change the zenith angle : 90=86, 80=>78
>       DO i=1,nt
>         IF (sza(i).eq.80.) sza(i)=78.
>         IF (sza(i).eq.90.) sza(i)=86.
>       ENDDO
> 
> 
536c503,520
> 
> * cba pchem3,4,IN added for jvalues in SGMCM
>       CALL pchem3(nw,wl,nz,tlev,airden,
>      $     nj4,sj4,jlabel4,jlabel5)
> 
> csl pchem4 added for jvalues for MCM2001
>       CALL pchem4(nw,wl,nz,tlev,airden,
>      $     nj6,sj6,jlabel6,jlabel7)
> 
> csl pchemIN added for jvalues for MCM2001
>       CALL pchemIN(nw,wl,nz,tlev,airden,
>      $     nj8,sj8,jlabel8,jlabel9)
> 
> cjlt pchem08 added for jvalues for difunctional nitrates
>       CALL pchem08(nw,wl,nz,tlev,airden,
>      $     nj8,sj8,jlabel8,jlabel9)
> 
677a647,651
> * cba following lines added for saprc99 and SGMCM 
>          CALL zero2(valj4,kj,kz)
>          CALL zero2(valj6,kj,kz)
>          CALL zero2(valj8,kj,kz)
> 
733a707,722
> * cba - loop for SGMCM
>                DO ij = 1, nj4
>                  djdw = saflux(iz,iw) * sj4(ij,iz,iw)
>                  valj4(ij,iz) = valj4(ij,iz) + djdw * (wu(iw) - wl(iw))
>                ENDDO
> 
>                DO ij = 1, nj6
>                  djdw = saflux(iz,iw) * sj6(ij,iz,iw)
>                  valj6(ij,iz) = valj6(ij,iz) + djdw * (wu(iw) - wl(iw))
>                ENDDO
> 
>                DO ij = 1, nj8
>                  djdw = saflux(iz,iw) * sj8(ij,iz,iw)
>                  valj8(ij,iz) = valj8(ij,iz) + djdw * (wu(iw) - wl(iw))
>                ENDDO
> 
753a742,757
> * cba - save the j-value data for izout
>          DO ij=1,nj4
>           resu4(it,ij)=valj4(ij,izout)
>          ENDDO
>          DO ij=1,nj6
>            resu6(it,ij)=valj6(ij,izout)
>          ENDDO
>          DO ij=1,nj8
>            resu8(it,ij)=valj8(ij,izout)
>          ENDDO
> 
> 
> 
>  20   CONTINUE
>         write(6,*) 'izout=', izout
>         write(6,*) 'z(izout)', z(izout)
780a785,821
> 
>       toto1 = 90.0
>       toto2 = 0.
> 
> * cba - write results for MCM2001 
> * jmlt - and write results for difunctional nitrates
>       DO ij=1,nj8
>         WRITE(90,'(a2)') '/ '
>         WRITE(90,'(a2,a40)') '/ ', jlabel8(ij)
>         WRITE(90,'(a40)') jlabel9(ij)
>         DO it=1,nt
>          WRITE(90,'(f6.1, 1pe12.4)') sza(it),resu8(it,ij)
>         ENDDO
>         WRITE(90,'(f6.1, 1pe12.4)') toto1, toto2
>       ENDDO
> 
>       DO ij=1,nj4
>         WRITE(90,'(a2)') '/ '
>         WRITE(90,'(a2,a40)') '/ ', jlabel4(ij)
>         WRITE(90,'(a40)') jlabel5(ij)
>         DO it=1,nt
>          WRITE(90,'(f6.1, 1pe12.4)') sza(it),resu4(it,ij)
>         ENDDO
>         WRITE(90,'(f6.1, 1pe12.4)') toto1, toto2
>       ENDDO
> 
>       DO ij=1,nj6
>         WRITE(90,'(a2)') '/ '
>         WRITE(90,'(a2,a40)') '/ ', jlabel6(ij)
>         WRITE(90,'(a40)') jlabel7(ij)
>         DO it=1,nt
>          WRITE(90,'(f6.1, 1pe12.4)') sza(it),resu6(it,ij)
>         ENDDO
>         WRITE(90,'(f6.1, 1pe12.4)') toto1, toto2
>       ENDDO
> 
>       WRITE(90,'(a4)') 'END '
