/[MITgcm]/MITgcm_contrib/jscott/igsm/src/r95mit.F
ViewVC logotype

Diff of /MITgcm_contrib/jscott/igsm/src/r95mit.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.1 by jscott, Fri Aug 11 19:35:31 2006 UTC revision 1.3 by jscott, Wed Sep 2 15:30:28 2009 UTC
# Line 1  Line 1 
1    C $Header$
2    C $Name$
3    
4  #include "ctrparam.h"  #include "ctrparam.h"
5    
# Line 199  C Line 201  C
201  #include "B83XX.COM"                                                     180.    #include "B83XX.COM"                                                     180.  
202  #include "chem_para"  #include "chem_para"
203  #include "chem_com"  #include "chem_com"
204    #if ( defined OCEAN_3D )
205          PARAMETER (IM0=01,JM0=N_LAT,LM0=N_LEV,IO0=N_LON0,KAIJ0=75)
206    #include "AGRID.h"
207    #endif
208  #if ( defined CLM )  #if ( defined CLM )
209  #include "CLM.COM"  #include "CLM.h"
210  #endif  #endif
211    
212        DIMENSION ALVISK(11,4),ALNIRK(11,4),FIELDC(11,3),VTMASK(11)        241.          DIMENSION ALVISK(11,4),ALNIRK(11,4),FIELDC(11,3),VTMASK(11)        241.  
# Line 444  C$    XEANIR=XEANIR+XSNNIR*(1.-EXPSNE) Line 450  C$    XEANIR=XEANIR+XSNNIR*(1.-EXPSNE)
450    
451  #if ( defined CLM )  #if ( defined CLM )
452  c     if(ncallclm.ge.1)then  c     if(ncallclm.ge.1)then
453        BEAVIS=0.7*asdirclm(JLAT)+0.3*asdifclm(JLAT)        i=1
454        BEANIR=0.7*aldirclm(JLAT)+0.3*aldifclm(JLAT)        BEAVIS=0.7*asdirclm(i,JLAT)+0.3*asdifclm(i,JLAT)
455          BEANIR=0.7*aldirclm(i,JLAT)+0.3*aldifclm(i,JLAT)
456        XEAVIS=BEAVIS        XEAVIS=BEAVIS
457        XEANIR=BEANIR        XEANIR=BEANIR
458  c     endif  c     endif
# Line 499  C Line 506  C
506        EXPSNO=EXP(-SNOWOI/DMOICE)                                         478.          EXPSNO=EXP(-SNOWOI/DMOICE)                                         478.  
507        BOIVIS=AOIVIS*EXPSNO+BSNVIS*(1.-EXPSNO)                            479.          BOIVIS=AOIVIS*EXPSNO+BSNVIS*(1.-EXPSNO)                            479.  
508        BOINIR=AOINIR*EXPSNO+BSNNIR*(1.-EXPSNO)                            480.          BOINIR=AOINIR*EXPSNO+BSNNIR*(1.-EXPSNO)                            480.  
509    Cnext lines - use albedos from ice model
510    #if ( defined OCEAN_3D )
511          BOIVIS=mmsAlb(JLAT)
512          BOINIR=mmsAlbNIR(JLAT)
513    #endif
514        XOIVIS=BOIVIS                                                      481.          XOIVIS=BOIVIS                                                      481.  
515        XOINIR=BOINIR                                                      482.          XOINIR=BOINIR                                                      482.  
516  C                                                                        483.    C                                                                        483.  
# Line 539  C Line 551  C
551    
552  #if ( defined CLM )  #if ( defined CLM )
553  c     if(ncallclm.ge.1)then  c     if(ncallclm.ge.1)then
554        BLIVIS=0.7*asdirclm(JLAT)+0.3*asdifclm(JLAT)        i=1
555        BLINIR=0.7*aldirclm(JLAT)+0.3*aldifclm(JLAT)        BLIVIS=0.7*asdirclm(i,JLAT)+0.3*asdifclm(i,JLAT)
556          BLINIR=0.7*aldirclm(i,JLAT)+0.3*aldifclm(i,JLAT)
557  c     endif  c     endif
558  #endif  #endif
559    
# Line 698  C Line 711  C
711  C                                                                        685.    C                                                                        685.  
712        IF(LASTVC.GE.0) GO TO 107                                          686.          IF(LASTVC.GE.0) GO TO 107                                          686.  
713  C                                                                        687.    C                                                                        687.  
       print *,' Before DO 105'  
714        DO 105 L=1,NL                                                      688.          DO 105 L=1,NL                                                      688.  
715        P=PLB(L)                                                           689.          P=PLB(L)                                                           689.  
716        DO 101 N=2,8                                                       690.          DO 101 N=2,8                                                       690.  
# Line 1809  C Line 1821  C
1821    
1822  #include "B83XX.COM"  #include "B83XX.COM"
1823  #if ( defined CLM )  #if ( defined CLM )
1824  #include "CLM.COM"  #include "CLM.h"
1825  #endif  #endif
1826        
1827        DATA R6,R24/.1666667,4.166667E-02/                                1577.          DATA R6,R24/.1666667,4.166667E-02/                                1577.  
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22